{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# Data Source: https://www.kaggle.com/worldbank/world-development-indicators\n",
    "# Folder: 'world-development-indicators'"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<br><p style=\"font-family: Arial; font-size:3.75em;color:purple; font-style:bold\">\n",
    "Matplotlib: Exploring <br> <br> <br>Data Visualization</p><br><br>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<br><br><center><h1 style=\"font-size:2em;color:#2467C0\">World Development Indicators</h1></center>\n",
    "<br>\n",
    "<table>\n",
    "<col width=\"550\">\n",
    "<col width=\"450\">\n",
    "<tr>\n",
    "<td><img src=\"https://upload.wikimedia.org/wikipedia/commons/4/46/North_South_divide.svg\" align=\"middle\" style=\"width:550px;height:360px;\"/></td>\n",
    "<td>\n",
    "This week, we will be using an open dataset from <a href=\"https://www.kaggle.com\">Kaggle</a>. It is  <a href=\"https://www.kaggle.com/worldbank/world-development-indicators\">The World Development Indicators</a> dataset obtained from the World Bank containing over a thousand annual indicators of economic development from hundreds of countries around the world.\n",
    "<br>\n",
    "<br>\n",
    "This is a slightly modified version of the original dataset from <a href=\"http://data.worldbank.org/data-catalog/world-development-indicators\">The World Bank</a>\n",
    "<br>\n",
    "<br>\n",
    "List of the <a href=\"https://www.kaggle.com/benhamner/d/worldbank/world-development-indicators/indicators-in-data\">available indicators</a> and a <a href=\"https://www.kaggle.com/benhamner/d/worldbank/world-development-indicators/countries-in-the-wdi-data\">list of the available countries</a>.\n",
    "</td>\n",
    "</tr>\n",
    "</table>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Step 1: Initial exploration of the Dataset"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import numpy as np\n",
    "import random\n",
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(5656458, 6)"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data = pd.read_csv('./world-development-indicators/Indicators.csv')\n",
    "data.shape"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This is a really large dataset, at least in terms of the number of rows.  But with 6 columns, what does this hold?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>CountryName</th>\n",
       "      <th>CountryCode</th>\n",
       "      <th>IndicatorName</th>\n",
       "      <th>IndicatorCode</th>\n",
       "      <th>Year</th>\n",
       "      <th>Value</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Arab World</td>\n",
       "      <td>ARB</td>\n",
       "      <td>Adolescent fertility rate (births per 1,000 wo...</td>\n",
       "      <td>SP.ADO.TFRT</td>\n",
       "      <td>1960</td>\n",
       "      <td>1.335609e+02</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>Arab World</td>\n",
       "      <td>ARB</td>\n",
       "      <td>Age dependency ratio (% of working-age populat...</td>\n",
       "      <td>SP.POP.DPND</td>\n",
       "      <td>1960</td>\n",
       "      <td>8.779760e+01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Arab World</td>\n",
       "      <td>ARB</td>\n",
       "      <td>Age dependency ratio, old (% of working-age po...</td>\n",
       "      <td>SP.POP.DPND.OL</td>\n",
       "      <td>1960</td>\n",
       "      <td>6.634579e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>Arab World</td>\n",
       "      <td>ARB</td>\n",
       "      <td>Age dependency ratio, young (% of working-age ...</td>\n",
       "      <td>SP.POP.DPND.YG</td>\n",
       "      <td>1960</td>\n",
       "      <td>8.102333e+01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Arab World</td>\n",
       "      <td>ARB</td>\n",
       "      <td>Arms exports (SIPRI trend indicator values)</td>\n",
       "      <td>MS.MIL.XPRT.KD</td>\n",
       "      <td>1960</td>\n",
       "      <td>3.000000e+06</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>Arab World</td>\n",
       "      <td>ARB</td>\n",
       "      <td>Arms imports (SIPRI trend indicator values)</td>\n",
       "      <td>MS.MIL.MPRT.KD</td>\n",
       "      <td>1960</td>\n",
       "      <td>5.380000e+08</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>Arab World</td>\n",
       "      <td>ARB</td>\n",
       "      <td>Birth rate, crude (per 1,000 people)</td>\n",
       "      <td>SP.DYN.CBRT.IN</td>\n",
       "      <td>1960</td>\n",
       "      <td>4.769789e+01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>Arab World</td>\n",
       "      <td>ARB</td>\n",
       "      <td>CO2 emissions (kt)</td>\n",
       "      <td>EN.ATM.CO2E.KT</td>\n",
       "      <td>1960</td>\n",
       "      <td>5.956399e+04</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>Arab World</td>\n",
       "      <td>ARB</td>\n",
       "      <td>CO2 emissions (metric tons per capita)</td>\n",
       "      <td>EN.ATM.CO2E.PC</td>\n",
       "      <td>1960</td>\n",
       "      <td>6.439635e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>Arab World</td>\n",
       "      <td>ARB</td>\n",
       "      <td>CO2 emissions from gaseous fuel consumption (%...</td>\n",
       "      <td>EN.ATM.CO2E.GF.ZS</td>\n",
       "      <td>1960</td>\n",
       "      <td>5.041292e+00</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "  CountryName CountryCode                                      IndicatorName  \\\n",
       "0  Arab World         ARB  Adolescent fertility rate (births per 1,000 wo...   \n",
       "1  Arab World         ARB  Age dependency ratio (% of working-age populat...   \n",
       "2  Arab World         ARB  Age dependency ratio, old (% of working-age po...   \n",
       "3  Arab World         ARB  Age dependency ratio, young (% of working-age ...   \n",
       "4  Arab World         ARB        Arms exports (SIPRI trend indicator values)   \n",
       "5  Arab World         ARB        Arms imports (SIPRI trend indicator values)   \n",
       "6  Arab World         ARB               Birth rate, crude (per 1,000 people)   \n",
       "7  Arab World         ARB                                 CO2 emissions (kt)   \n",
       "8  Arab World         ARB             CO2 emissions (metric tons per capita)   \n",
       "9  Arab World         ARB  CO2 emissions from gaseous fuel consumption (%...   \n",
       "\n",
       "       IndicatorCode  Year         Value  \n",
       "0        SP.ADO.TFRT  1960  1.335609e+02  \n",
       "1        SP.POP.DPND  1960  8.779760e+01  \n",
       "2     SP.POP.DPND.OL  1960  6.634579e+00  \n",
       "3     SP.POP.DPND.YG  1960  8.102333e+01  \n",
       "4     MS.MIL.XPRT.KD  1960  3.000000e+06  \n",
       "5     MS.MIL.MPRT.KD  1960  5.380000e+08  \n",
       "6     SP.DYN.CBRT.IN  1960  4.769789e+01  \n",
       "7     EN.ATM.CO2E.KT  1960  5.956399e+04  \n",
       "8     EN.ATM.CO2E.PC  1960  6.439635e-01  \n",
       "9  EN.ATM.CO2E.GF.ZS  1960  5.041292e+00  "
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "data.head(10)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Looks like it has different indicators for different countries with the year and value of the indicator. "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### How many UNIQUE country names are there ?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "247"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "countries = data['CountryName'].unique().tolist()\n",
    "len(countries)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Are there same number of country codes ?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "247"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# How many unique country codes are there ? (should be the same #)\n",
    "countryCodes = data['CountryCode'].unique().tolist()\n",
    "len(countryCodes)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Are there many indicators or few ?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "1344"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# How many unique indicators are there ? (should be the same #)\n",
    "indicators = data['IndicatorName'].unique().tolist()\n",
    "len(indicators)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### How many years of data do we have ?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "56"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# How many years of data do we have ?\n",
    "years = data['Year'].unique().tolist()\n",
    "len(years)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### What's the range of years?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "1960  to  2015\n"
     ]
    }
   ],
   "source": [
    "print(min(years),\" to \",max(years))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<p style=\"font-family: Arial; font-size:2.5em;color:blue; font-style:bold\">\n",
    "Matplotlib: Basic Plotting, Part 1</p><br>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Lets pick a country and an indicator to explore: CO2 Emissions per capita and the USA"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# select CO2 emissions for the United States\n",
    "hist_indicator = 'CO2 emissions \\(metric'\n",
    "hist_country = 'USA'\n",
    "\n",
    "mask1 = data['IndicatorName'].str.contains(hist_indicator) \n",
    "mask2 = data['CountryCode'].str.contains(hist_country)\n",
    "\n",
    "# stage is just those indicators matching the USA for country code and CO2 emissions over time.\n",
    "stage = data[mask1 & mask2]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>CountryName</th>\n",
       "      <th>CountryCode</th>\n",
       "      <th>IndicatorName</th>\n",
       "      <th>IndicatorCode</th>\n",
       "      <th>Year</th>\n",
       "      <th>Value</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>22232</th>\n",
       "      <td>United States</td>\n",
       "      <td>USA</td>\n",
       "      <td>CO2 emissions (metric tons per capita)</td>\n",
       "      <td>EN.ATM.CO2E.PC</td>\n",
       "      <td>1960</td>\n",
       "      <td>15.999779</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>48708</th>\n",
       "      <td>United States</td>\n",
       "      <td>USA</td>\n",
       "      <td>CO2 emissions (metric tons per capita)</td>\n",
       "      <td>EN.ATM.CO2E.PC</td>\n",
       "      <td>1961</td>\n",
       "      <td>15.681256</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>77087</th>\n",
       "      <td>United States</td>\n",
       "      <td>USA</td>\n",
       "      <td>CO2 emissions (metric tons per capita)</td>\n",
       "      <td>EN.ATM.CO2E.PC</td>\n",
       "      <td>1962</td>\n",
       "      <td>16.013937</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>105704</th>\n",
       "      <td>United States</td>\n",
       "      <td>USA</td>\n",
       "      <td>CO2 emissions (metric tons per capita)</td>\n",
       "      <td>EN.ATM.CO2E.PC</td>\n",
       "      <td>1963</td>\n",
       "      <td>16.482762</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>134742</th>\n",
       "      <td>United States</td>\n",
       "      <td>USA</td>\n",
       "      <td>CO2 emissions (metric tons per capita)</td>\n",
       "      <td>EN.ATM.CO2E.PC</td>\n",
       "      <td>1964</td>\n",
       "      <td>16.968119</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "          CountryName CountryCode                           IndicatorName  \\\n",
       "22232   United States         USA  CO2 emissions (metric tons per capita)   \n",
       "48708   United States         USA  CO2 emissions (metric tons per capita)   \n",
       "77087   United States         USA  CO2 emissions (metric tons per capita)   \n",
       "105704  United States         USA  CO2 emissions (metric tons per capita)   \n",
       "134742  United States         USA  CO2 emissions (metric tons per capita)   \n",
       "\n",
       "         IndicatorCode  Year      Value  \n",
       "22232   EN.ATM.CO2E.PC  1960  15.999779  \n",
       "48708   EN.ATM.CO2E.PC  1961  15.681256  \n",
       "77087   EN.ATM.CO2E.PC  1962  16.013937  \n",
       "105704  EN.ATM.CO2E.PC  1963  16.482762  \n",
       "134742  EN.ATM.CO2E.PC  1964  16.968119  "
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "stage.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Let's see how emissions have changed over time using MatplotLib"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAD8CAYAAABn919SAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAADddJREFUeJzt3X+s3fVdx/HnGwrRwRaL3LKK1MsMQ5s4GLsiZourEhhj\nid1i0LG5NYKpfzCzmU2tPxJNjAkuOn9Eo6mhW4mT+WNbqBGH2DgbHUPaBUpZ2QqkYqXhFotCNBEK\nb/84n47jtefee37ce8553+cjuTnf8znf8z2f9z2nr/s5n++PRmYiSZp+Z427A5Kk0TDQJakIA12S\nijDQJakIA12SijDQJakIA12SijDQJakIA12Sili3mi924YUX5uzs7Gq+pCRNvQMHDjybmTNLrbeq\ngT47O8v+/ftX8yUlaepFxL8sZz2nXCSpCANdkoow0CWpCANdkoow0CWpCANdkoow0CWpCANdkoow\n0CWpiFU9U1Sra3bHX/+f+0dvf9eYeiJpNThCl6QiDHRJKsJAl6QiDHRJKsKdohqphTtiwZ2x0mox\n0LUqDHpp5RnoBXh4oiRwDl2SynCEvsaMaurDKRRp8jhCl6QiDHRJKsIpFy3KqRVpehjoE8gQlTQI\nA32KGPSSFuMcuiQV4QhdY+VJUdLoOEKXpCIMdEkqYskpl4i4BLgTeD3wCrAzM383Ii4A/gyYBY4C\nP5qZz61cV7WS3OEqTb/ljNBPAR/NzO8GrgFui4jNwA5gb2ZeBuxt9yVJY7LkCD0zjwPH2/ILEXEY\nuBjYCmxpq+0Gvgj8/Ir0UtLU8Nve+PR1lEtEzAJvBh4ALmphT2Yej4gNI++dNIU8ckfjsuxAj4jz\ngc8CH8nM5yNiuc/bDmwH2LRp0yB9lCaSwa1Js6xAj4hz6IT5pzPzc635mYjY2EbnG4H5Mz03M3cC\nOwHm5uZyBH0uw6+mkkZpOUe5BHAHcDgzP9H10B5gG3B7u717RXqoNcnR73j083t3QDJ5ljNCfyvw\nAeCRiHiotf0inSD/84i4FXgKuGlluihpnFb6P0Xpd/u9/ug4CFjeUS7/CPSaML92tN2RtBqmIfz8\nBtA/zxSVpCK8OJe0hFGMFMc1beEod20x0CWtOVX/0Bnoq6Dqh0fDc2S9stba79FAV1lr7R+zZKBr\nzTHoVZVHuUhSEY7QpcaRu6adgT5CBsLKcgeitDgDXZKaaR8cOIcuSUUY6JJUhFMuA5j2r2WSanKE\nLklFGOiSVISBLklFGOiSVISBLklFeJTLEqbhv+qSNB6TdsSbI3RJKsJAl6QinHKRpCVM2tRKLwa6\nJI3YuP4AGOiNOz8lTTvn0CWpCANdkoow0CWpiDU3h+5cuaSqHKFLUhEGuiQVYaBLUhFl59Cn5cwu\nSRqVqQl0A1qSFueUiyQVYaBLUhFTM+XSi1MxktThCF2Silgy0CNiV0TMR8ShrrZfjYh/i4iH2s+N\nK9tNSdJSljNC/xRwwxnafzszr2w/94y2W5Kkfi0Z6Jm5Dzi5Cn2RJA1hmDn0D0XEwTYls35kPZIk\nDWTQQP9D4DuBK4HjwG/1WjEitkfE/ojYf+LEiQFfTpK0lIECPTOfycyXM/MV4I+BqxdZd2dmzmXm\n3MzMzKD9lCQtYaBAj4iNXXffAxzqta4kaXUseWJRRNwFbAEujIhjwK8AWyLiSiCBo8BPrWAfJUnL\nsGSgZ+bNZ2i+YwX6IkkagmeKSlIRBrokFWGgS1IRBrokFWGgS1IRBrokFWGgS1IRBrokFWGgS1IR\nBrokFWGgS1IRBrokFWGgS1IRBrokFWGgS1IRBrokFWGgS1IRBrokFWGgS1IRBrokFWGgS1IRBrok\nFWGgS1IRBrokFWGgS1IRBrokFWGgS1IRBrokFWGgS1IRBrokFWGgS1IRBrokFWGgS1IRBrokFWGg\nS1IRBrokFWGgS1IRSwZ6ROyKiPmIONTVdkFE3BcRR9rt+pXtpiRpKcsZoX8KuGFB2w5gb2ZeBuxt\n9yVJY7RkoGfmPuDkguatwO62vBt494j7JUnq06Bz6Bdl5nGAdrthdF2SJA1ixXeKRsT2iNgfEftP\nnDix0i8nSWvWoIH+TERsBGi3871WzMydmTmXmXMzMzMDvpwkaSmDBvoeYFtb3gbcPZruSJIGtZzD\nFu8C7gcuj4hjEXErcDtwXUQcAa5r9yVJY7RuqRUy8+YeD1074r5IkobgmaKSVISBLklFGOiSVISB\nLklFGOiSVISBLklFGOiSVISBLklFGOiSVISBLklFGOiSVISBLklFGOiSVISBLklFGOiSVISBLklF\nGOiSVISBLklFGOiSVISBLklFGOiSVISBLklFGOiSVISBLklFGOiSVISBLklFGOiSVISBLklFGOiS\nVISBLklFGOiSVISBLklFGOiSVISBLklFGOiSVISBLklFGOiSVMS6YZ4cEUeBF4CXgVOZOTeKTkmS\n+jdUoDc/mJnPjmA7kqQhOOUiSUUMG+gJ/G1EHIiI7aPokCRpMMNOubw1M5+OiA3AfRHxWGbu616h\nBf12gE2bNg35cpKkXoYaoWfm0+12Hvg8cPUZ1tmZmXOZOTczMzPMy0mSFjFwoEfEeRHx2tPLwPXA\noVF1TJLUn2GmXC4CPh8Rp7fzp5n5hZH0SpLUt4EDPTOfBK4YYV8kSUPwsEVJKsJAl6QiDHRJKsJA\nl6QiDHRJKsJAl6QiDHRJKsJAl6QiDHRJKsJAl6QiDHRJKsJAl6QiDHRJKsJAl6QiDHRJKsJAl6Qi\nDHRJKsJAl6QiDHRJKsJAl6QiDHRJKsJAl6QiDHRJKsJAl6QiDHRJKsJAl6QiDHRJKsJAl6QiDHRJ\nKsJAl6QiDHRJKsJAl6QiDHRJKsJAl6QiDHRJKsJAl6QiDHRJKmKoQI+IGyLiaxHxeETsGFWnJEn9\nGzjQI+Js4A+AdwKbgZsjYvOoOiZJ6s8wI/Srgccz88nMfBH4DLB1NN2SJPVrmEC/GPjXrvvHWpsk\naQwiMwd7YsRNwDsy8yfb/Q8AV2fmTy9Ybzuwvd29HPja4N39hguBZ0ewnUm3VuoEa61qrdS60nV+\nR2bOLLXSuiFe4BhwSdf9bweeXrhSZu4Edg7xOv9PROzPzLlRbnMSrZU6wVqrWiu1Tkqdw0y5PAhc\nFhGXRsS5wHuBPaPpliSpXwOP0DPzVER8CLgXOBvYlZmPjqxnkqS+DDPlQmbeA9wzor70Y6RTOBNs\nrdQJ1lrVWql1IuoceKeoJGmyeOq/JBUxEYEeEbsiYj4iDnW1XRER90fEIxHxVxHxuq7H3tQee7Q9\n/k2t/S3t/uMR8XsREeOoZzH91BoR74+Ih7p+XomIK9tj1Wo9JyJ2t/bDEfELXc+Z6EtM9FnnuRHx\nydb+cERs6XrONLynl0TE37f36NGI+HBrvyAi7ouII+12fWuPVsvjEXEwIq7q2ta2tv6RiNg2rprO\nZIA6v6u93/8TER9bsK3V+/xm5th/gB8ArgIOdbU9CLy9Ld8C/FpbXgccBK5o978VOLst/zPw/UAA\nfwO8c9y1DVPrgud9D/Bk1/1StQLvAz7Tll8DHAVm6exwfwJ4A3Au8DCwedy1DVHnbcAn2/IG4ABw\n1hS9pxuBq9rya4Gv07n0x8eBHa19B/AbbfnGVksA1wAPtPYLgCfb7fq2vH7c9Q1R5wbge4FfBz7W\ntZ1V/fxOxAg9M/cBJxc0Xw7sa8v3AT/Slq8HDmbmw+25/56ZL0fERuB1mXl/dn6TdwLvXvne96fP\nWrvdDNwFULTWBM6LiHXANwMvAs8zBZeY6LPOzcDe9rx54D+AuSl6T49n5lfa8gvAYTpniG8FdrfV\ndvNq37cCd2bHl4FvabW+A7gvM09m5nN0fkc3rGIpi+q3zsycz8wHgZcWbGpVP78TEeg9HAJ+uC3f\nxKsnMb0RyIi4NyK+EhE/19ovpnOy02nTdCmCXrV2+zFaoFOz1r8E/gs4DjwF/GZmnmR6LzHRq86H\nga0RsS4iLgXe0h6buvc0ImaBNwMPABdl5nHohCGdESv0fv+m5n1dZp29rGqdkxzotwC3RcQBOl95\nXmzt64C3Ae9vt++JiGvpfKVbaFoO4elVKwAR8X3Af2fm6TnairVeDbwMfBtwKfDRiHgD01trrzp3\n0flHvR/4HeBLwCmmrM6IOB/4LPCRzHx+sVXP0JaLtE+UPursuYkztK1YnUMdh76SMvMxOtMrRMQb\ngXe1h44B/5CZz7bH7qEzf/kndC4/cNoZL0UwiRap9bT38uroHDq/g2q1vg/4Qma+BMxHxD8Bc3RG\nN0teYmLS9KozM08BP3N6vYj4EnAEeI4peU8j4hw6IffpzPxca34mIjZm5vE2pTLf2ntdIuQYsGVB\n+xdXst/96rPOXpZ1iZRRmdgRekRsaLdnAb8M/FF76F7gTRHxmjbf+nbgq+3rzwsRcU07OuCDwN1j\n6HrfFqn1dNtNdObegG981atW61PAD7WjIs6jswPtMab0EhO96myf2/Pa8nXAqcycms9v69sdwOHM\n/ETXQ3uA00eqbOPVvu8BPtje12uA/2y13gtcHxHr25Ei17e2iTBAnb2s7ud33HuT257gu+jMnb5E\n5y/arcCH6exZ/jpwO+0kqLb+jwOP0pmn/HhX+1xrewL4/e7nTMrPALVuAb58hu2UqhU4H/iL9r5+\nFfjZru3c2NZ/Avilcdc1ZJ2zdK44ehj4OzpX0Zum9/RtdKYMDgIPtZ8b6RxttpfOt429wAVt/aDz\nH+E8ATwCzHVt6xbg8fbzE+Oubcg6X9/e++fp7Og+Rmcn96p+fj1TVJKKmNgpF0lSfwx0SSrCQJek\nIgx0SSrCQJekIgx0SSrCQJekIgx0SSrifwEFy1A8cPlM6wAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f1160e43940>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# get the years\n",
    "years = stage['Year'].values\n",
    "# get the values \n",
    "co2 = stage['Value'].values\n",
    "\n",
    "# create\n",
    "plt.bar(years,co2)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Turns out emissions per capita have dropped a bit over time, but let's make this graphic a bit more appealing before we continue to explore it."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEWCAYAAACJ0YulAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd8VfX9x/HXJwl7h733ECmCMhU3zlqtWnFr68Bat1Zr\nW2utra3a1lq7fmLdigNHUcGJCxFQQPaQPcNIGAkBMj+/P85JjDHjZtx7k/B+Ph73ce/9nvU5l3A/\n93zPd5i7IyIiApAQ7wBERKTmUFIQEZFCSgoiIlJISUFERAopKYiISCElBRERKaSkIActM7vYzN6r\nwvZHm9mK6owpgmO+bWaXx/KYcnBRUpAqMbOLzGyOme01s5TwS2tMkeUDzewNM9tjZhlm9pGZHVlk\neT8zm2xmO8xsp5m9a2b9yzjeU2aWHR6v4LGgMrG7+/PufnJltg23n+7upcYaDe5+mrs/XZltzczN\nrE+xsnvM7Lki739lZmvDz3WTmb1Uwn6eMrNcM+tUmTikZlNSkEozs1uBh4E/Au2BbsC/gbPC5b2B\nGcAioCfQCXgdeM/MRoe7aQm8AfQP9/EFMLmcQz/o7k2LPA6r1hM7SIVXIJcCY929KTAMmFZsnSbA\nucAe4OKYBylRp6QglWJmLYB7gevc/TV3z3T3HHd/091vD1e7B5jp7r92953unuHujwDPAg8AuPsX\n7v54uDwH+BvQ38xaVyKmHuGv4Z+Y2UYz22VmPzWz4Wa20Mx2m9k/i6z/YzP7LHxtZvY3M9seXtUs\nNLNB4bLTzWxpeKWz2cx+HpYfZ2abiuzvEDP7ODzOEjM7s8iyp8zsX2Y2JdzP7DBplnnsEs7xYzO7\nqmj8ZvaX8FzXmtlpFf3cihgOvOvuqwHcfau7Tyi2zrnAboJ/e1Vj1UFKClJZo4GGBL/8S3MSMKmE\n8peBo8yscQnLjgG2untaFWIbCfQFzie4kvk1MBY4FBhnZseWsM3J4bH7EVy9nA8UxPA4cI27NwMG\nAR8W39jM6gFvAu8B7YAbgOeLVYVdCPwOaAWsAu6L4NiRnOsKoA3wIPC4mVmE2xY3C7jMzG43s2Fm\nlljCOpcDLwAvAgPM7PBKHktqKCUFqazWQKq755axThsgpYTyFIK/vVZFC82sC/Av4NZyjv3z8Nd4\nwaN4Hfvv3f2Au78HZAIvuPt2d98MTAeGlrDPHKAZMAAwd1/m7ilFlg00s+buvsvd55Ww/SigKXC/\nu2e7+4fAWwSJoMBr4ZVRLvA8MCSCY5dnvbs/5u55wNNAR4JquApz9+cIktkpwCfAdjO7s2C5mXUD\njgcmuvs2gqolXS3UMUoKUllpQBszSypjnVSCL6niOgL5wK6CAjNrS/Ar+9/u/kI5x/6Lu7cs8ij+\nxbStyOv9JbxvWnyH4Zf4PwmS0jYzm2BmzcPF5wKnA+vN7JMi90OK6gRsdPf8ImXrgc5F3m8t8npf\nQRzlHLs8hft0933hy++cXygPqFesrB5BUirYx/PuPpbgiuWnwL1mdkq4+FJgmbvPD98/D1wUXiVJ\nHaGkIJU1EzgA/LCMdT4AziuhfBzBvYZ9AGbWiiAhvOHu95Wwfky4+yPufgRBNVM/4Paw/Et3P4ug\nWuh/BNVfxW0BuppZ0f9T3YDNVTl2NdsA9ChW1pMgeRWPJ8fdJwELCarMAC4DepnZVjPbCjxEcDVY\nlfsYUsMoKUiluPse4G7gX2b2QzNrbGb1zOw0M3swXO13wJFmdp+ZJZtZMzO7geDL5RcA4S/id4EZ\n7n5nSceKhfBm9MjwV28mQcLLM7P6FvRnaBHeCE8n+MVd3OxwuzvCz+E44AcEde+VOnb1nNm3vATc\nZWZdzCzBzMaGMb4SxvFjM/t++O+UEN60PhSYHV4d9QZGEFR7DSFIFhNRFVKdoqQglebuDxHU/98F\n7AA2AtcT/JrG3VcCY4DDgHUE9xLOBU5x9xnhbs4maPXyE/t234NuZRz6jmLrplbD6TQHHiOo0lpP\nUD32l3DZpcA6M0snqFK5pPjG7p4NnEnwqzmVoGnuZe6+vIrHrk73Ap8Dn4XHehC42N0Xh8vTgV8R\nXFHsDpdf6+6fEXzxT3b3RWGrpK3uvhX4O3CGmSVHIV6JA9MkOyIiUkBXCiIiUihqScHMulowpMGy\nsCPPTWH5PWEHoPnh4/RoxSAiIhUTteojM+sIdHT3eWbWDJhL0FJlHLDX3aNRZyoiIlVQVhvzKgk7\n36SErzPMbBnfbrMtIiI1TExuNJtZD+BTgiZstwI/JmjpMAe4zd13lbDNeGA8QJMmTY4YMGBA1OMU\nEalL5s6dm+rubSuyTdSTgpk1Jegyf5+7v2Zm7Qma7Dnwe4IqpivK2sewYcN8zpw5UY1TRKSuMbO5\n7j6sIttEtfVR2BnnVeB5d38NwN23uXteOBzAYwSdYUREpAaIZusjIxhdclnYyamgvOhYOGcDi4tv\nKyIi8RG1G83AUQQ9QReZWcEAWr8CLjSzIQTVR+uAa6IYg4iIVEA0Wx99BpQ0rvvUaB1TRESqRj2a\nRUSkkJKCiIgUUlIQEZFCSgoiIlJISUFERAopKYiISCElBRERKRRxUjCzJmaWGM1gREQkvkpNCuHE\n3ReZ2RQz2w4sB1LCCXP+bGZ9YxemiIjEQllXCh8BvYFfAh3cvau7twOOBmYB95vZdyYwFxGR2qus\nYS7GuntO8UJ330kw8umr4SioIiJSR5SaFIonBDNrBzQssnxDSUlDRERqr3JvNJvZmWa2ElhLMFnO\nOuDtKMclIiJxEEnro98Do4Cv3b0ncCIwI6pRiYhIXESSFHLcPQ1IMLMEd/8IGBLluEREJA4imU9h\ndzjP8qfA82Hz1NzohiWxNnH2BnZkZHHjiX0IJs0TkYNRJEnhLGA/cAtwMdAC+F00g5LY+u/0Nfxh\nyjIAOrVsyHnDusY5IhGJl0iqj+5293x3z3X3p939EeAX0Q5MYuPZmev4w5RlnDaoA6N6JXP35CWs\n2r433mGJSJxEkhROKqHstOoORGLvpS838JvJSxh7SDv+fsFQ/n7BUBrVT+T6ifM4kJMX7/BEJA7K\nGubiWjNbBPQ3s4VFHmuBhbELUaLhtXmbuPO1RRzTry3/uvhw6icl0L55Q/563mEs35rBH6cui3lM\nefke82OKyLeVdU9hIkF/hD8BdxYpzwh7NUst9eaCLfx80gJG92rNhEuPoEHSN+McHj+gHVcf3ZPH\npq/lyN6tOXVQx6jHs3rHXu55YwkzV6cxuEsLjurThtG9W3N4t1Y0rKcxGEViydxL/nVmZs3dPd3M\nkktaHsvEMGzYMJ8zZ06sDlenvbN4K9dNnMcR3Vrx1BXDaVz/u78LsnPz+dH/fc661Eym3nQ0XVo1\njkosmVm5/OPDVTz+2Roa1kvkrCGdWLw5nYWbdpPv0CApgeE9kjmyT2suHtGdFo01qopIRZjZXHcf\nVqFtykgKb7n7GWF1kQNF2ym6u/eqfKgVo6RQPT5avp3xz85hUOcWPHvlSJo2KP1CcX1aJt9/5DP6\ntW/KS9eMpl5i9U294e68uTCFP05Zxtb0A5x3RBfuOHUAbZs1ACD9QA5frNnJjNWpfL4qjRXbMhjV\nK5nnrxpFYoKay4pEqlqTQk2ipFB1M1al8pOnvqRf+6Y8f9UoWjQq/1f3Gwu2cOMLX3Hd8b25/ZQB\n1RLHiq0Z/PaNxcxas5NBnZtz71mDOLxbqzK3mTRnI7e/spCbTuzLLSf1q5Y4RA4GlUkKkfRTwMzO\nAcYQXDFMd/f/VSI+iZMv1+3kqqfn0LN1E569YmRECQHgzMM6MWNlKv/+eDVj+rRldO/WlY5h9Y69\nPDJtJW8s2EKLRvW47+xBXDC8W0S//M8b1pWZa9J45MOVjOyZzJF92lQ6DhEpW7lXCmb2b6AP8EJY\ndD6w2t2vi3JshXSlUHkLNu7m4v/Opl3zBrw0fnRhFU2k9mfncfoj08nJy+edm48ps8qpJOtSM3lk\n2kr+N38zDZISufzIHlxzTC9aNalfof1kZuVy5j8/I/1ALlNvPLrC5yFyMIpK9ZGZLQEGebiimSUA\ni9z90EpHWkFKCpWzdEs6Fz42ixaN6vHyNaPp0KJh+RuVYM66nZz36EwuGtGN+87+XkTbbEjbxyMf\nruT1rzZTL9G4bHQPxh/TizZNK/9lvnxrOmf9cwbDeyTzzBUjSND9BZEyRav6aAXQDVgfvu+K+inU\nCM/MXMeBnDwGdGjOgA7NaNusQeG4RSu3ZXDJ47NpUj+R568aWemEADCsRzJXjQmaqZ46qANH921b\n5vovf7mRX72+iMQE4/LRPfjpcb1o16zyxy8woENz7jnzUH752iL+/fEqrj9BM8KKVLdIkkJrYJmZ\nfRG+Hw7MNLM3ANz9zGgFJ6X7aPl27p685FtlrRrXY0CH5vTv0Iypi1JITDCev3oUXZOr3qT0tpP7\nM235dn7xykLeveUYmjUs+b7Eq3M38YvXFjKmTxv+ct5htG9e9WRQ1AXDuzJzdRoPvf81w3skM7JX\n5e9ziMh3RVJ9dGxZy939k2qNqASqPvq2rNw8TvnbpySY8eL4UazekcmKreks35rB8q0ZrNiaQfNG\nSTx35Uj6tm9Wbcf9asMuzv3P54wb1pX7zx38neWT52/mlpfmM7p3ax6/fHjUOp7tzcrljEemsz8n\nj6k3Hk3rKlRJidRlUak+isWXvlTM45+tZV3aPp6+YgTtmjekXfOG32oZlJ/vOFR7m/6h3VpxzbG9\n+c/Hqzl1UAeO69+ucNmUhSnc8tJ8hvdI5r+XRS8hADRtkMQ/Lzqcc/79ObdNWsATlw/X/QWRahLJ\ndJyjzOxLM9trZtlmlmdm6bEITr4rZc9+/vnhKk4a2J5j+5Vct5+QYFHr5HXz2L70a9+UO19dxJ79\nwRTd7yzeyo0vfsUR3VvxxI+H06h+9IemGNS5BXedcQgfr9jBEzPWRv14IgeLSLqp/hO4EFgJNAKu\nCsskDv44dTm5+c7dZwyMy/EbJCXyl/MOY8feLO59cykfLN3GDS/MY3CXFjz5kxE0qWCT1aq4dFR3\nThrYngfeWc7izXtidlyRuiyisQvcfRWQ6O557v4kcFxUo5ISzVqTxpsLtvDTY3tXy83jyhrcpSU/\nO643r87bxE+fm8vAjs15+ooRFe7DUFVmxoPnDqZ1kwbc+MJXZGZpQkCRqookKewzs/rAfDN70Mxu\nAZpEOS4pJjcvn3veWELnlo249tje8Q6HG07oy2FdWnBop+Y8c8VImpfSGinaWjWpz8MXDGFtWib3\nvLGk/A1EpEyRJIVLw/WuBzIJ+imcW95GZtbVzD4ys2VmtsTMbgrLk83sfTNbGT6XPfCNAPD87A0s\n35rBXd8/JCZ19uWpn5TAK9ceyes/Oyruo5eO6tWaG47vw6S5m5g8f3NcY5HYST+Qw5+mLuO/09eQ\nr7k4qk0k1/upQLa7HwB+Z2aJQCRtAHOB29x9npk1A+aa2fvAj4Fp7n6/md1JMFeDpvcsQ9reLP76\n3grG9GnDqYM6xDucQtU5cmpV3XhiX2asTuOu1xcztGsrurWOX/VaZbk7e/bnkLo3ix0Z2aRlZpGa\nkUXq3mz6tm/KWUM6xzvEGsHdmbIohXvfXMr2jCwApq9M5aFxh6l5cjWIJClMA8YCBRP3NgLeA44s\nayN3TwFSwtcZZrYM6AycxTf3JJ4GPkZJoUx/fncF+7LzuOfMgYU9luXbkhIT+PsFQzjt79O58cWv\nmPTT6h3uO9q+WLuTnz0/l9S92aWuk7LnAD+tAVWH8bQhbR+/mbyYT77ewaDOzfnv5cNYuGkP9761\nlNMfmc7fLxjKKHVorJJIkkJDdy+cyd3d95pZhX6GmVkPYCgwG2gfJgzcPcXM2pWyzXhgPEC3bt0q\ncrg6Zc66nbw0ZyNXHtWTPu2qryNaXdSlVWPuP2cw102cx9/e/5o7Tq2e4b6jbeqiFG5+aT5dWzXi\n2uP60KZpfdo0bUCbpg1o3bQ+zRvW47ZJC7j/7eUkJRhXHR2zqUxqjOzcfB6bvoZHpq0kKcG4+4yB\nXDa6O0mJCQzu0pKh3Vpy/cSvuOixWdw8th/XHd9Hc29UUiRJIdPMDnf3eQBmdgSwP9IDmFlT4FXg\n5nAmt4i2c/cJwAQIejRHery6ZOmWdK58eg5dWjXiprEa5ycS3x/ckekru/KfT1ZzTL+2Nf5X41Mz\n1vK7t5ZyeLdWPH75MFo2Lnn02L+NO4z8fOcPU5aRYMYVY3rGONKqc3dWbd/L/I27adGoHt1bN6Fr\ncqMSZ//Lys1j1fa9LE/JYPnWdD5cvp3VOzI59dAO/PbMgXRs0ehb6x/aqQVv3jCGX7++iIfe/5rZ\na9P42/lDqmXMrdy8fLZlZNG4XiItG9er81frkQxzMRx4EdgSFnUEznf3ueXu3Kwe8Bbwrrs/FJat\nAI4LrxI6Ah+7e/+y9nMwDnOxclsG50+YRYOkBF6+ZnRcm6DWNvuyczn14ekkJRhTbzq6Rs7znJ/v\nPPDuch79ZA0nD2zPIxcOLTfOnLx8rp84j3eXbOPesw7lstE9YhNsJbk7q3dkMnNNGrPWpDF7TVqJ\n1WNtmjage+vGdEtuTF6+s3xrOqt3ZJIX3jyun5TAIR2aceOJfTnxkPblHvPlORv57RtLMIyWjetR\nLzGBeokWPgevWzSqR3KTBrRpWp/kJvVpHV6VJSUY69L2sS41k3WpmaxNy2Tjzn3k5IWxJCbQtlkD\n2jVvQLtmDWjfvCFDurbkh0M618he9VGbeS38cu9PMCXncnfPiWAbI7hnsNPdby5S/mcgrciN5mR3\nv6OsfR1sSWFdaibjHp2JAy+NH0Wvtk3jHVKtM33lDi59/AtuPLEvt9aw2dqyc/O545UF/G/+Fi4d\n1Z17zjw04qqO7Nx8fvb8PD5Yto0//HAQl4zqHuVoK2Z/dh6frtzBe0u28cnXO0jdG9wI7hAOxTK6\nV2sO796KzKxcNuzcFzzS9hW+BjikY7PCgR0P6diMHq2bkFTB+0MrtmYwcfZ69ufkkZPnZOflk5Ob\nT26+k52bz579OaTtzSI1M5vs3PzvbN+wXgI9WjehR+sm9GzbhG7Jjdmfncf2jCy2px8InjMOsHXP\nAdIP5DK6V2se/NHgGvfjrUZNx2lmY4DpwCKg4FP/FcF9hZcJhuPeAJzn7jvL2tfBlBQ27drHuP+b\nyYHcfF4cP4p+1Tig3cHm5he/YsqiFN6+6egacz8m40AO1z43j89WpXL7Kf352XG9K1wdkZWbx7XP\nzePD5du5/5zvccGIqt9z25+dx6y1aexIzyI338nNzyc3L3zOdxLNaN+8IR1aNKRji4a0b96w8Mpm\nV2Y205Zv570lW/l05Q4O5OTTvGESx/Vvx5G9WzOqV2u6t25cI6td3J3M7LwgQezNJicvn+6tG9O+\nWcOIfvm7Oy9+uZH7piwj351fnn4IF4/oVmOuGmpUUqhOB0tS2LrnAOMencnufdlMvHoUgzq3iHdI\ntVrq3izGPvQJ/do148Xxo+L+H3Xz7v1c+dSXrNq+l/vPHcyPjuhS6X1l5eZxzbNz+XjFDsb0acNZ\nQzpxyqAOFepEuGnXPj5avp0Pl2/n89VpZJXwi7ksyU2Cqpe1qUFVT4fmDTn50PacPLADI3sl16rW\nX1W1efd+7nx1IdNXpnJk79Y8cG7JVw2ZWbms2ZFJl1aNKjz7YGVUe1IIq4C6uPvGqgZXFQdDUtiR\nkcX5E2ayPT2LZ68cwdByJrOXyLz85UbueHUhfzrne1xYDb+oK2vBxt1c+fQcsnLz+PfFh5c7UVEk\nDuTk8e+PV/O/rzazYec+6iclMPaQdpx5WGeO69+28Jf8/uw8UvbsJ2XPAbbs3s+qHXv5ePkOVmzL\nAKBbcmNOGNCOEwa0o1fbJiQlJJCUaCQlGEmJCSQlGLn5zrb0oLokZc8Btob725aeRf8OTTl5YAcG\nd2lRI68GYsXdeeGLjdw3ZSkAt57cn2YNkli5PYOV2/eyctteNu8O2ui0aVqfF64eVa1D25ckWtNx\nznX3I6oUWRXV9aQwd/0ufj5pASl79vPMFSMZ0TM53iHVGe7OBRNmsSwlnQ9uO7ZaWqNU1NuLUrjl\n5fm0bdaAJ388vNqrstydrzbu5o35W3hr4RZS92bTrGESXVo1JmXPfnbv+/YtwKQEY3iPZE4Y0I7j\nB7Sjd9smB/WXeXXbtGsfd766iM9WpQLBjfLebZvSt13w6NyqEX96eznu8OL4UfRpF717htFKCv8C\nnnL3L6sSXFXU1aRwICePv763gv9+tpaOzRvy8AVDlRCiYPWOvZz28HROGdSBf1w4NOLtCnrOPvTe\n1+zYm0XXVo3p0qoRXZMb07VVI7q0akyPNo3p1aZpiVVT7s5/PlnNg++s4PBuLZlw2bAqzVEdidy8\nfGasDgZO3JmZTccWDenUshEdWzSkY4tGdGoZ3BdokFTzWmTVJe7O/I27adW4Pl2TG3+nIcGq7Xu5\nYMIszILE0DtKjUmilRSWErQ8Wkcw9pEB7u7fnXorSupiUpizbid3vLKQNamZXDSyG788bUCpU1xK\n1T38wdc8/MFKnvrJ8G9NDlSaWWvS+NPUZSzYtIcBHZoxvEcym3btY9Ou/WzctY8DOd/UvzdrkMRh\nXYMOVEO7tWRI11Y0bZDEr19fxKS5mzjzsE48+KPBNbJprMTPym0ZXPjYLBLMeOma0fRsU/3jjEYr\nKZTY5s3d11fkQFVRl5LC/uw8/vzuCp78fC2dWjTigXMHM6Zvm3iHVedl5eZx+t+nk5Wbz3u3HFNi\nhykI/qM+8M5yPli2nQ7NG3Lbyf045/Au3/ql5+6k7s1m0659rN6RyfyNu/hqw26Wb80obFvfolE9\n9uzP4aYT+3Lz2L6qnpESrdgaJIb6iQm8OH4UPao5MUSzn8IYoK+7P2lmbYGm7h6z6a7qSlKYtSaN\nO19dyLq0fVw6qju/OG1AzOcgOJjNXpPG+RNmMaJHMj3bNKF+UtCZqX5SAvUTjc27D/D6V5toUj+J\na4/vzRVH9azQr/t92bks2rSHrzbuZumWdE4+tD1nDO4UxTOSumD51nQunDCLhvUSeXH8KLq3rr7E\nEK0rhd8Cw4D+7t7PzDoBk9z9qMqHWjG1PSkEQ/wu54UvNtA1Obg6OLK3rg7i4ZFpK5k0dyPZuflB\np6bc/OCRl0/9xAQuHtWNG07oS3IMmguKFFi6JZ2L/juLxvUSefOGMdU22mu0ksJ8gsHs5rn70LBs\noe4pROa9JVv5zeTF7MjI4oqjenLryf1KrbqQ+HF38h0NoiZxs2jTHn747xlcMLwr9539vWrZZ2WS\nQiTfTtnu7mbm4UE061oEdmRkcc8bS5iyKIUBHZox4dJhHNa1ZbzDklKYGYnKBxJH3+vSgktHdeeZ\nmeu4dHR3BnRoHpc4Iuly+LKZPQq0NLOrgQ+Ax6IbVu321sItjH3oE95fuo3bTurHG9ePUUIQkXLd\nPLYvzRrW4/dvLSVeo02UmxTc/S/AKwTDX/cD7nb3f0Q7sNpqwqeruX7iV/Rq24SpN43hhhP7Uj/p\n4OnuLyKV17JxfW4Z25cZq9KYtmx7XGKI9NtqEcHgdp+Gr6WY/HznvilL+ePU5Xx/cMewp2LNGIRN\nRGqPi0d1p3fbJtw3dVmJI7hGW7lJwcyuAr4AzgF+BMwysyuiHVhtkpOXz22TFvDY9LVcPro7/7hg\nqHqMikil1EtM4K4zBrI2NZNnZq6L+fEjudF8OzDU3dMAzKw18DnwRDQDqy0ys3K59vl5fPr1jkoP\nhSwiUtTx/dtxbL+2/H3aSs4e2rnamqhGIpLqo01ARpH3GUBcR02tKdL2ZnHRY7OYsSqVB88dzHXH\n91FCEJFqcdf3D2Ffdh5/++DrmB43kiuFzcBsM5sMOHAW8IWZ3QpQMM3mwWbx5j3c8MJXbNm9n0cv\nOYKxA8ueJlBEpCL6tm/GJSO78eys9Vw6qgf9O3xzj3J/dh5TFqXw4hcbWJeWyY0n9uXikd2rpZ9N\nJFcKq4H/ESQEgMlACtAsfBxUDuTk8cA7yznrXzPIzMpl4tUjlRBEJCpuHtvvW01Ul25J5+7Jixnx\nxw/4+aQFpGVm06tNU+6evIRz//M5y1LSq3xMzbxWAXPW7eSOVxeyZkcm44Z14denD6RFY41sKiLR\n8+SMtfzuzaX0bdeUldv3Uj8pgdMHdeCCEd0YGQ61P3n+Fn7/1lJ278/hqqN7cvOJ/WhUPzFqPZoP\neplZufz53RU8PXMdnVo04tkrR1TLzFkiIuW5ZFR3/vfVZrJy8/ntDwZy9tDOtGz87bG5fji0M8f2\na8uf3l7Go5+sYcrCFP7ww0GVOp6uFMrx2cpU7nxtIZt37+fy0T24/ZT+NNHIpiJSQ81ak8avXl/E\nmh2ZrH/gjApfKairbSnSD+Rw56sLueTx2dRPTODla0Zzz5mHKiGISI02qldr3r7paG4e27dS20fS\nee1BM2tuZvXMbJqZpZrZJZU6Wi0xbdk2Tn7oU16es5Frju3F1JuOZngPTZMpIrVDg6REbh7br1Lb\nRvKz92R3v8PMzibos3Ae8BHwXKWOWIPtyszmd28u4X/zt9C/fTMevfQIDWQnIgeVSJJCQfOa04EX\n3H1nXeug5e68vXgrd09ezO59wRSK1x3fRwPZichBJ5Kk8KaZLQf2Az8Lp+M8EN2wYmfL7v3cPXkJ\nHyzbxqDOzXn2ypEc0jE+45iLiMRbuUnB3e80sweAdHfPM7NMgl7NtVpevvPcrPU8+M5y8tz51ekD\nuOKoniQl6upARA5ekTalOQToYWZF138mCvHExPKt6dz56iLmb9zN0X3b8Mezv0fX5MbxDktEJO7K\nTQpm9izQG5gP5IXFTpyTQlZuHlv3HKB768hnB92fncc/PlzJhE/X0LxRPR4+fwhnDemkQexEREKR\nXCkMAwZ6Devl9uaCFG5/ZQEn9G/Hj4/qwZg+bUr9ct+ecYBnPl/Pc7PXs3tfDj86ogu/Pv0QWjWp\nX+L6IiIHq0iSwmKgA8EgeDXGMX3bcMMJfZk4ez2XPv4Ffdo15fIje3DO0M6FHcxWbM3gv9PXMHn+\nFnLy8zkcMcBYAAAS/UlEQVTpkPZcc2wvjuiuPgciIiUpd5gLM/sIGEIw+1pWQbm7nxnd0L5R1jAX\nWbl5TFmYwpMz1rFo8x6aNUzi3MO7sCY1k0+/3kHDegmcd0RXrhjTk55tIq9qEhGp7aI1IN49lQsn\nNhokJXLO4V04e2hn5m3YxZMz1vHsrPUkN6nPz0/ux8Uju6uaSEQkQpE0Sf3EzNoDw8OiL9x9e3TD\nqjgz44juyRzRPZn0Azk0TEpU5zMRkQqKZOyjcQRVR+cB4whmYftRtAOriuYN6ykhiIhUQiTfnL8G\nhrv75e5+GTAC+E15G5nZE2a23cwWFym7x8w2m9n88HF65UMXEZHqFklSSChWXZQW4XZPAaeWUP43\ndx8SPqZGsB8REYmRSG40v2Nm7wIvhO/PB94ubyN3/9TMelQ+NBERibVyf/G7++3Ao8Bg4DBggrvf\nUYVjXm9mC8PqpValrWRm481sjpnN2bFjRxUOJyIikYrkRvMD7v6au9/q7re4++vhAHmV8R+CITOG\nEHSG+2tpK7r7BHcf5u7D2rbVfMgiIrEQyb2Bk0ooO60yB3P3be6e5+75wGMEN61FRKSGKPWegpld\nC/wM6GVmC4ssagbMqMzBzKyjuxcMl3E2wRAaIiJSQ5R1o3kiwQ3lPwF3FinPcPed5e3YzF4AjgPa\nmNkm4LfAcWY2hGCU1XXANZULW0REoqHUpODue4A9wIWV2bG7l7Td45XZl4iIxIa6/YqISCElBRER\nKRRJk9SeZtawyPtG6pQmIlI3RXKlMAnIL/I+LywTEZE6JpKkkOTu2QVvwteaoEBEpA6KJCnsMLPC\nWdbM7CwgNXohiYhIvEQyIN5PgefN7J+AARuBy6IalYiIxEUkM6+tBkaZWVOCOZ0zoh+WiIjEQ1nD\nXFzi7s+Z2a3FygFw94eiHJuIiMRYWVcKTcLnZrEIRERE4q+sYS4eNbNEIN3d/xbDmEREJE7KbH3k\n7nnAmWWtIyIidUckrY8+D1sevQRkFhS6+7yoRSUiInERSVI4Mny+t0iZAydUfzgiIhJPkSSFK919\nTdECM+sVpXhERCSOIunR/EoJZRr7SESkDiqrn8IA4FCghZmdU2RRc6BhyVuJiEhtVlb1UX/gDKAl\n8IMi5RnA1dEMSkRE4qOsfgqTgclmNtrdZ8YwJhERiZNI7imkmdk0M1sMYGaDzeyuKMclIiJxEElS\neAz4JZAD4O4LgQuiGZSIiMRHJEmhsbt/UawsNxrBiIhIfEWSFFLNrDdBhzXM7EdASlSjEhGRuIik\n89p1wARggJltBtYCl0Q1KhERiYtIJtlZA4w1syZAgibZERGpu8pNCmbWkmD6zR5AUpFJdm6MamQi\nIhJzkVQfTQVmAYuA/OiGIyIi8RRJUmjo7reWv5qIiNR2kbQ+etbMrjazjmaWXPCIemQiIhJzkVwp\nZAN/Bn5N2Cw1fNbw2SIidUwkSeFWoI+7p0Y7GBERia9Iqo+WAPuiHYiIiMRfJFcKecB8M/sIyCoo\nVJNUEZG6J5Kk8L/wISIidVwkPZqfjkUgIiISf6XeUzCzN83sB2ZWr4RlvczsXjO7ooztnzCz7QXz\nMIRlyWb2vpmtDJ9bVf0URESkupR1o/lq4GhguZl9aWZTzexDM1sLPArMdfcnytj+KeDUYmV3AtPc\nvS8wLXwvIiI1RFnTcW4F7gDuMLMeQEdgP/C1u5fbGsndPw23K+os4Ljw9dPAx8AvKhayiIhESyQ3\nmnH3dcC6ajhee3dPCfeZYmbtSlvRzMYD4wG6detWDYcWEZHyRNJPIS7cfYK7D3P3YW3bto13OCIi\nB4VYJ4VtZtYRIHzeHuPji4hIGSqUFMyslZkNrsLx3gAuD19fDkyuwr5ERKSalZsUzOxjM2sejoy6\nAHjSzB6KYLsXgJlAfzPbZGZXAvcDJ5nZSuCk8L2IiNQQkdxobuHu6WZ2FfCku//WzBaWt5G7X1jK\nohMrFKGIiMRMJNVHSWH9/zjgrSjHIyIicRRJUrgXeBdY5e5fmlkvYGV0wxIRkXiIZOyjScCkIu/X\nAOdGMygREYmPcpOCmbUlGPKiR9H13b3UcY9ERKR2iuRG82RgOvABwdwKIiJSR0WSFBq7u8YnEhE5\nCERyo/ktMzs96pGIiEjcRZIUbiJIDAfMLCN8pEc7MBERib1IWh81i0UgIiISfxENnW1mZwLHhG8/\ndnd1YhMRqYMiGfvofoIqpKXh46awTERE6phIrhROB4a4ez6AmT0NfIWm0hQRqXMiHTq7ZZHXLaIR\niIiIxF8kVwp/Ar4ys48AI7i38MuoRiUiInERSeujF8zsY2A4QVL4hbtvjXZgIiISe6VWH5nZgPD5\ncKAjsAnYCHQKy0REpI4p60rhVmA88NcSljlwQlQiEhGRuCk1Kbj7+PD5+NiFIyIi8RRJP4XzzKxZ\n+PouM3vNzIZGPzQREYm1SJqk/sbdM8xsDHAK8DTwf9ENS0RE4iGSpFAwh8L3gf+4+2SgfvRCEhGR\neIkkKWw2s0eBccBUM2sQ4XYiIlLLRPLlPg54FzjV3XcDycDtUY1KRETiIpIezR2BKe6eZWbHAYOB\nZ6IalYiIxEUkVwqvAnlm1gd4HOgJTIxqVCIiEheRJIV8d88FzgEedvdbCK4eRESkjokkKeSY2YXA\nZUDB5Dr1oheSiIjESyRJ4SfAaOA+d19rZj2B56IbloiIxEMko6QuBW4s8n4toJnXRETqoFKTgpm9\n7O7jzGwRwQB4hYsAd/fBUY9ORERiqqwrhZvC5zNiEYiIiMRfWaOkpoTP6wHMrHlZ64uISO1X7pe8\nmV0D3Avs55tqJAd6RTEuERGJg0h++f8cONTdU6MdjIiIxFckSWE1sK86D2pm64AMghFYc919WHXu\nX0REKieSpPBL4HMzmw1kFRS6+42lbxKR43X1ISJSs0SSFB4FPgQWAfnRDUdEROIpkqSQ6+63VvNx\nHXjPzBx41N0nFF/BzMYD4wG6detWzYcXEZGSRDLMxUdmNt7MOppZcsGjisc9yt0PB04DrjOzY4qv\n4O4T3H2Yuw9r27ZtFQ8nIiKRiORK4aLw+ZdFyqrUJNXdt4TP283sdWAE8Gll9yciItUjkrGPelbn\nAc2sCZDg7hnh65MJ+kGIiEicRdJ5rTFwK9DN3cebWV+gv7u/Vc6mpWkPvG5mBcef6O7vVHJfIiJS\njSKpPnoSmAscGb7fBEzim7kVKsTd1wCHVWZbERGJrkhuNPd29weBHAB3308wUqqIiNQxkSSFbDNr\nRDjukZn1pkgnNhERqTsiqT76LfAO0NXMngeOAn4czaBERCQ+Iml99L6ZzQNGEVQb3aThKURE6qaI\n5kdw9zRgSpRjERGROIvknoKIiBwklBRERKSQkoKIiBQqNSmY2ffMbJaZbTSzCWbWqsiyL2ITnoiI\nxFJZVwr/Ae4Bvgd8DXwW9lEAqBfluEREJA7Kan3UtMiYRH8xs7nAO2Z2KWFHNhERqVvKSgpmZi3c\nfQ+Au39kZucCrwJVnU9BRERqoLKqjx4ADila4O4LgROB16IZlIiIxEepVwruPrHgtZk1DYo80903\nAFfHIjgREYmtMpukmtm1ZrYBWA9sNLP1Zvaz2IQmIiKxVlaT1LuAHwDHuXtrd08GjgdOC5eJiEgd\nU9aVwqXAOeGkOEDhBDnjgMuiHZiIiMRemdVH7n6ghLL9QH7UIhIRkbgpKylsMrMTixea2QlASvRC\nEhGReCmrn8KNwGQz+4xgjmYHhhNMsnNWDGITEZEYK/VKwd2XAIOAT4EeQK/w9aBwmYiI1DGlXimY\nWR+gvbs/Uaz8aDPb4u6rox6diIjEVFn3FB4GMkoo3x8uExGROqaspNAjHNbiW9x9DkF1koiI1DFl\nJYWGZSxrVN2BiIhI/JWVFL40s++McWRmVxK0RhIRkTqmrCapNwOvm9nFfJMEhgH1gbOjHZiIiMRe\nWaOkbgOONLPjCZqmAkxx9w9jEpmIiMRcWVcKQDC5DvBRDGIREZE4K3PsIxERObgoKYiISCElBRER\nKaSkICIihZQURESkkJKCiIgUiktSMLNTzWyFma0yszvjEYOIiHxXzJOCmSUC/wJOAwYCF5rZwFjH\nISIi3xWPK4URwCp3X+Pu2cCLaCY3EZEaodwezVHQGdhY5P0mYGTxlcxsPDA+fLvXzFaUsK82QGq1\nR1i7HOyfwcF+/qDPQOdf+vl3r+jO4pEUrIQy/06B+wRgQpk7Mpvj7sOqK7Da6GD/DA728wd9Bjr/\n6j3/eFQfbQK6FnnfBdgShzhERKSYeCSFL4G+ZtbTzOoDFwBvxCEOEREpJubVR+6ea2bXA+8CicAT\n7r6kkrsrs3rpIHGwfwYH+/mDPgOdfzUy9+9U54uIyEFKPZpFRKSQkoKIiBSqcUnBzJ4ws+1mtrhI\n2WFmNtPMFpnZm2bWvMiyweGyJeHyhmH5EeH7VWb2iJmV1BS2xqnI+ZvZxWY2v8gj38yGhMtq5flD\nhT+Demb2dFi+zMx+WWSbWjmcSgXPv76ZPRmWLzCz44psUyv/Bsysq5l9FP57LjGzm8LyZDN738xW\nhs+twnILz2+VmS00s8OL7OvycP2VZnZ5vM6pIipx/gPCv40sM/t5sX1V/P+Au9eoB3AMcDiwuEjZ\nl8Cx4esrgN+Hr5OAhcBh4fvWQGL4+gtgNEG/iLeB0+J9btV9/sW2+x6wpsj7Wnn+lfgbuAh4MXzd\nGFgH9CBoxLAa6AXUBxYAA+N9blE4/+uAJ8PX7YC5QEJt/hsAOgKHh6+bAV8TDInzIHBnWH4n8ED4\n+vTw/AwYBcwOy5OBNeFzq/B1q3ifXxTOvx0wHLgP+HmR/VTq/0CNu1Jw90+BncWK+wOfhq/fB84N\nX58MLHT3BeG2ae6eZ2YdgebuPtODT+cZ4IfRj77qKnj+RV0IvABQm88fKvwZONDEzJKARkA2kE4t\nHk6lguc/EJgWbrcd2A0Mq81/A+6e4u7zwtcZwDKCkRDOAp4OV3uab87nLOAZD8wCWobnfwrwvrvv\ndPddBJ/bqTE8lUqp6Pm7+3Z3/xLIKbarSv0fqHFJoRSLgTPD1+fxTee3foCb2btmNs/M7gjLOxN0\nkiuwKSyrrUo7/6LOJ0wK1L3zh9I/g1eATCAF2AD8xd13UvJwKrX5Myjt/BcAZ5lZkpn1BI4Il9WJ\nvwEz6wEMBWYD7d09BYIvToJfyFD6v3Wt/xuI8PxLU6nzry1J4QrgOjObS3A5lR2WJwFjgIvD57PN\n7EQiHEqjFint/AEws5HAPncvqIOua+cPpX8GI4A8oBPQE7jNzHpR9z6D0s7/CYL/7HOAh4HPgVzq\nwPmbWVPgVeBmd08va9USyryM8lqhAudf6i5KKCv3/OMx9lGFuftygqoizKwf8P1w0SbgE3dPDZdN\nJaiLfY5g+IwCtXoojTLOv8AFfHOVAMHnUmfOH8r8DC4C3nH3HGC7mc0AhhH8Qqozw6mUdv7ungvc\nUrCemX0OrAR2UYv/BsysHsEX4vPu/lpYvM3MOrp7Slg9tD0sL23onE3AccXKP45m3NWlgudfmkoN\nKVQrrhTMrF34nADcBfxfuOhdYLCZNQ7rlI8FloaXVhlmNipscXEZMDkOoVeLMs6/oOw8gvpCoPDS\nss6cP5T5GWwATghboDQhuNG4nDo2nEpp5x/+7TcJX58E5Lp7rf4/EMb7OLDM3R8qsugNoKAF0eV8\ncz5vAJeFfwOjgD3h+b8LnGxmrcKWOieHZTVaJc6/NJX7PxDvO+0l3Hl/gaB+OIcg010J3ERwB/5r\n4H7Cntjh+pcASwjqXB8sUj4sLFsN/LPoNjX5UYnzPw6YVcJ+auX5V/QzAJoCk8K/gaXA7UX2c3q4\n/mrg1/E+ryidfw9gBcHNyA+A7rX9b4CgKtgJWhbODx+nE7QunEZwJTQNSA7XN4KJu1YDi4BhRfZ1\nBbAqfPwk3ucWpfPvEP6dpBM0NNhE0MigUv8HNMyFiIgUqhXVRyIiEhtKCiIiUkhJQURECikpiIhI\nISUFEREppKQgUkzY3v0zMzutSNk4M3snnnGJxIKapIqUwMwGEfR/GEow2uR84FR3X12FfSZ50ANZ\npMZSUhAphZk9SDDYXhMgw91/H47Jfx3BUMSfA9e7e76ZTSAYYqUR8JK73xvuYxPwKMHonA+7+6Q4\nnIpIxGrF2EcicfI7YB7B4HPDwquHs4Ej3T03TAQXABMJxrnfGQ638pGZveLuS8P9ZLr7UfE4AZGK\nUlIQKYW7Z5rZS8Bed88ys7EEk5nMCYanoRHfDE18oZldSfB/qhPBPAcFSeGl2EYuUnlKCiJlyw8f\nEIyx84S7/6boCmbWl2BsohHuvtvMngMaFlklMyaRilQDtT4SidwHwDgzawNgZq3NrBvQHMgA0ovM\n+CVSK+lKQSRC7r7IzH4HfBAOYZ0D/JRggpulBCOSrgFmxC9KkapR6yMRESmk6iMRESmkpCAiIoWU\nFEREpJCSgoiIFFJSEBGRQkoKIiJSSElBREQK/T/IJdbEsNfpxAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f115fe0d5c0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# switch to a line plot\n",
    "plt.plot(stage['Year'].values, stage['Value'].values)\n",
    "\n",
    "# Label the axes\n",
    "plt.xlabel('Year')\n",
    "plt.ylabel(stage['IndicatorName'].iloc[0])\n",
    "\n",
    "#label the figure\n",
    "plt.title('CO2 Emissions in USA')\n",
    "\n",
    "# to make more honest, start they y axis at 0\n",
    "plt.axis([1959, 2011,0,25])\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Using Histograms to explore the distribution of values\n",
    "We could also visualize this data as a histogram to better explore the ranges of values in CO2 production per year. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# If you want to just include those within one standard deviation fo the mean, you could do the following\n",
    "# lower = stage['Value'].mean() - stage['Value'].std()\n",
    "# upper = stage['Value'].mean() + stage['Value'].std()\n",
    "# hist_data = [x for x in stage[:10000]['Value'] if x>lower and x<upper ]\n",
    "\n",
    "# Otherwise, let's look at all the data\n",
    "hist_data = stage['Value'].values"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "52\n"
     ]
    }
   ],
   "source": [
    "print(len(hist_data))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYwAAAEWCAYAAAB1xKBvAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XucnGV99/HPl0SCsOFkYA0EjSBgkUcou4CKh11RGvKg\nQYuVFClUbCrWitbUc4HUarVGrRYtRkG0YhZFo4hUyCMuASVAFgOEcj6VmAjKeTkn/J4/7mtlmMzM\nXpvdmXsm+b5fr3ntfZr7+s69u/Ob+3SNIgIzM7PRbFF2ADMz6wwuGGZmlsUFw8zMsrhgmJlZFhcM\nMzPL4oJhZmZZXDCs7Um6XlJf2TkMJN0p6Y1l57ByuGBYqWq9AUk6XtJlI+MR8fKIGBxlPTMlhaTJ\nTYraVOk1r5c0XPXYpexsZiM68p/LrNUkTY6IdU1u5vKIeE2T2zDbaN7DsLZXuRci6SBJKyQ9LOke\nSV9Miy1LPx9Mn8xfJWkLSZ+UdJekeyV9R9J2Fev9qzTvPkn/VNXOqZLOlfRdSQ8Dx6e2L5f0oKS1\nkk6TtGXF+kLSeyXdIukRSZ+StEd6zsOSvl+5/Bhe/x6S7pd0QBrfRdIfRg7TSfprSTekNm+X9LcV\nz+2TtFrSh9M2WCvpSEmzJd2c1vvxiuVHXvc5aX1XS9qvTq4tJH1U0m1pG35f0o5jfX3WOVwwrNN8\nGfhyRGwL7AF8P01/Xfq5fUR0RcTlwPHp0Q/sDnQBpwFI2gf4GnAMMB3YDti1qq05wLnA9sDZwHrg\ng8A04FXAocB7q54zC+gBXgl8GFiU2tgN2BeYO9YXHBG3AR8Bzpa0NfAt4KyKw3T3AkcA2wJ/DXxp\npLgkLwS2Sq/vZOAbwDtTztcCJ0vavep1/wDYEfge8GNJz6sR7f3AkcDrgV2AB4CvjvX1WQeJCD/8\nKO0B3AkMAw9WPB4DLqta5o1peBmwAJhWtZ6ZQACTK6b9AnhvxfjewNMUh2JPBhZXzNsaeKqinVOB\nZaNk/wCwpGI8gEMqxoeAj1SMfwH49zrrOh5YV7Udbqta5jzgOuBaYEqDXD8GTkrDfcDjwKQ0PjXl\nPLgq55EVr3t5xbwtgLXAa2v8Lm4ADq1YdvrI9i3778qP5jy8h2Ht4MiI2H7kwYaf2iudAOwF3Cjp\nKklHNFh2F+CuivG7KIpFd5p398iMiHgMuK/q+XdXjkjaS9L5kn6XDlN9hmJvo9I9FcOP1xjvapB3\neeV2iIg9quZ/g2Iv5T8i4smKXIdLWp4OLz0IzK7KdV9ErK/IUCtnZa7K7fIMsJpie1V7MbAkHaJ7\nkKKArKfYvrYJcsGwjhIRt0TEXGBn4HPAuZK2ofjUXG0NxZvaiBdRfIq/h+JT84yRGZKeD7ygurmq\n8f8EbgT2jOKQ2McBbfyrySepC/h34Azg1JFzBZKmAD8EFgLdqeBeMM5cu1W0uwXFdlpTY7m7gcOr\nitxWEfHbcbRtbcwFwzqKpHdK2il98n0wTV4P/B54huJcxYjFwAclvSS94X4GOCeKq53OBd4s6dXp\nRPQCRn+TnQo8DAxLehlw4oS9sNF9GRiKiHcDPwNOT9O3BKZQvP51kg4HDhtnWz2S3pYuUf4A8CSw\nvMZypwOflvRiAEk7SZozzratjblgWKeZBVwvaZjiTfToiHgiHVL6NPCrdIjklcCZwH9RnPe4A3gC\n+HuAiLg+DQ9Q7G08QnHy+Enqmw/8ZVr2G8A5E/zaXlXjPowD05vwLOA9abl/AA6QdExEPEJx8vn7\nFCed/5LiXMd4/AR4R1rfscDbIuLpGst9ObV1kaRHKIrKweNs29qYIvwFSmZpD+RBisNNd5SdpyyS\nTgVeGhHvLDuLtR/vYdhmS9KbJW2dzoEspLgC6c5yU5m1LxcM25zNoTiZuwbYk+Lwlne5zerwISkz\nM8viPQwzM8uySXU+OG3atJg5c2ZL23z00UfZZpttWtrmeHVa5k7LC52XudPyQudlbte8Q0NDf4iI\nnXKW3aQKxsyZM1mxYkVL2xwcHKSvr6+lbY5Xp2XutLzQeZk7LS90XuZ2zSvprtGXKviQlJmZZXHB\nMDOzLC4YZmaWxQXDzMyyuGCYmVkWFwwzM8vigmFmZllcMMzMLIsLhpmZZdmk7vQ2a1da0JJvct1A\nnOLORW3ieA/DzMyyNG0PQ9KZwBHAvRGxb5p2DrB3WmR74MGI2L/Gc++k+BrM9cC6iOhtVk4zM8vT\nzENSZwGnAd8ZmRAR7xgZlvQF4KEGz++PiD80LZ2ZmY1J0wpGRCyTNLPWPEkC/gJ4Q7PaNzOzidXU\nb9xLBeP8kUNSFdNfB3yx3qEmSXcADwABfD0iFjVoYx4wD6C7u7tnYGBgYsJnGh4epqurq6Vtjlen\nZe60vLBh5qG1Q6Xk6Jnek7XcprCN21275u3v7x/KPexf1lVSc4HFDeYfEhFrJO0MLJV0Y0Qsq7Vg\nKiaLAHp7e6PV/c23ax/3jXRa5k7LCxtm7l/QX0qOmJv3gXBT2MbtrtPy1tLyq6QkTQbeBpxTb5mI\nWJN+3gssAQ5qTTozM6unjMtq3wjcGBGra82UtI2kqSPDwGHAqhbmMzOzGppWMCQtBi4H9pa0WtIJ\nadbRVB2OkrSLpAvSaDdwmaRrgCuBn0XEz5uV08zM8jTzKqm5daYfX2PaGmB2Gr4d2K9ZuczMbOP4\nTm8zM8vigmFmZllcMMzMLIsLhpmZZXHBMDOzLC4YZmaWxQXDzMyyuGCYmVkWFwwzM8vigmFmZllc\nMMzMLIsLhpmZZXHBMDOzLC4YZmaWxQXDzMyyuGCYmVkWFwwzM8vigmFmZllcMMzMLIsLhpmZZWla\nwZB0pqR7Ja2qmHaqpN9KWpkes+s8d5akmyTdKumjzcpoZmb5mrmHcRYwq8b0L0XE/ulxQfVMSZOA\nrwKHA/sAcyXt08ScZmaWoWkFIyKWAfdvxFMPAm6NiNsj4ilgAJgzoeHMzGzMFBHNW7k0Ezg/IvZN\n46cCxwMPAyuAD0XEA1XPOQqYFRHvTuPHAgdHxPvqtDEPmAfQ3d3dMzAw0IyXUtfw8DBdXV0tbXO8\nOi1zp+WFDTMPrR0qJUfP9J6s5TaFbdzu2jVvf3//UET05iw7udlhqvwn8Ckg0s8vAO+qWkY1nle3\nqkXEImARQG9vb/T19U1I0FyDg4O0us3x6rTMnZYXNszcv6C/lBwxN+8D4aawjdtdp+WtpaVXSUXE\nPRGxPiKeAb5Bcfip2mpgt4rxGcCaVuQzM7P6WlowJE2vGH0rsKrGYlcBe0p6iaQtgaOB81qRz8zM\n6mvaISlJi4E+YJqk1cApQJ+k/SkOMd0J/G1adhfgmxExOyLWSXofcCEwCTgzIq5vVk4zM8vTtIIR\nEXNrTD6jzrJrgNkV4xcAG1xya2Zm5fGd3mZmlsUFw8zMsrhgmJlZFhcMMzPL4oJhZmZZXDDMzCyL\nC4aZmWVxwTAzsywuGGZmlsUFw8zMsrhgmJlZFhcMMzPL4oJhZmZZXDDMzCyLC4aZmWVxwTAzsywu\nGGZmlsUFw8zMsrhgmJlZFhcMMzPL0rSCIelMSfdKWlUx7fOSbpR0raQlkrav89w7JV0naaWkFc3K\naGZm+Zq5h3EWMKtq2lJg34h4BXAz8LEGz++PiP0jordJ+czMbAyaVjAiYhlwf9W0iyJiXRpdDsxo\nVvtmZjaxFBHNW7k0Ezg/IvatMe+nwDkR8d0a8+4AHgAC+HpELGrQxjxgHkB3d3fPwMDAxITPNDw8\nTFdXV0vbHK9Oy9xpeWHDzENrh0rJ0TO9J2u5TWEbt7t2zdvf3z+UeyRncrPD1CLpE8A64Ow6ixwS\nEWsk7QwslXRj2mPZQComiwB6e3ujr6+vGZHrGhwcpNVtjlenZe60vLBh5v4F/aXkiLl5Hwg3hW3c\n7jotby0tv0pK0nHAEcAxUWf3JiLWpJ/3AkuAg1qX0MzMamlpwZA0C/gI8JaIeKzOMttImjoyDBwG\nrKq1rJmZtU4zL6tdDFwO7C1ptaQTgNOAqRSHmVZKOj0tu4ukC9JTu4HLJF0DXAn8LCJ+3qycZmaW\np2nnMCJibo3JZ9RZdg0wOw3fDuzXrFxmZrZxfKe3mZllccEwM7MsLhhmZpbFBcPMzLK4YJiZWZZR\nC4ak50tSGt5D0mxJpdwhbmZm5cnZw7gUeL6k6cAlwInAmU1NZWZmbSenYGyR7sr+c+C0iHgz8Irm\nxjIzs3aTVTAkHQj8JXB+mjapeZHMzKwd5RSMDwILKLroWCVpd4rDVGZmthlpePJa0iTgzyJi9si0\n1HXHe5sdzMzM2kvDPYyIWI+7FjczM/I6H7xa0o+AHwCPjkyMiPOalsrMzNpOTsHopigUsyumBeCC\nYWa2GRm1YETEsa0IYmZm7W3UgiFpCnA88HJgq5HpETGvebHMzKzd5FxW+x1gJsX3cF8B7AE80cRM\nZmbWhnIKxl4R8TFgOCLOAGYB+zY3lpmZtZucgvF0+vmgpD+h+E7uFzcvkpmZtaOcq6TOkLQDcApw\nIbB1GjYzs83IqHsYEfH1iHggIn4ZES+KiGkR8dWclUs6U9K9klZVTNtR0lJJt6SfO9R57nFpmVsk\nHZf/kszMrBlyvg9jJ0lfl3R+Gt9H0vGZ6z+L4pxHpY8Cv4iIPYFfpPHqNnek2Is5mOJO81PqFRYz\nM2uNnHMYZ1F8D8ZuafwW4EM5K4+IZcD9VZPnAN9Ow98Gjqzx1D8DlkbE/RHxALCUDQuPmZm1kCKi\n8QLSVRFxoKTfRMSfpmkrI2L/rAakmcD5EbFvGn8wIravmP9AROxQ9Zz5wFYR8S9p/J+AxyNiYY31\nzwPmAXR3d/cMDAzkxJoww8PDdHV1tbTN8eq0zJ2WFzbMPLR2qJQcPdN7spbbFLZxu2vXvP39/UMR\n0ZuzbM5J70fTIaIASN+N8cg48uVQjWk1K1tELAIWAfT29kZfX18TY21ocHCQVrc5Xp2WudPywoaZ\n+xf0l5Ij5jb+QDhiU9jG7a7T8taSc0hqPvBTYHdJlwCLgb8fR5v3pK97Jf28t8Yyq3n2EBjADGDN\nONo0M7NxqlswJG0LEBErgH7g9cBJwD4RsXIcbZ4HjFz1dBzwkxrLXAgcJmmHdLL7sDTNzMxK0mgP\n4zeSjgaIiKci4pqIWBkRT+WuXNJi4HJgb0mrJZ0AfBZ4k6RbgDelcST1Svpmau9+4FPAVenxz2ma\nmZmVpNE5jDcA/57e5E+MiFvHuvKImFtn1qE1ll0BvLti/EzgzLG2aWZmzVG3YETEXcBbJc0CfiXp\nKuCZivlvaUE+MzNrE6N9p/fewIeBS4GvUlEwzMxs81K3YEj6LPAW4EMR8d+ti2RmZu2o0R7GeuCA\niPB3X5iZWcNzGJ9oZRAzM2tvOXd6m20ytKBWJwITb+FeC0u7u9usWRrduHdI+jmldXHMzKxdNbpx\n7yvp5+WtCGJmZu2t0SGppyV9C9hV0leqZ0bE+5sXy8zM2k2jgnEE8EaKO77L6ZvZzMzaRqOrpP4A\nDEi6ISKuaWEmMzNrQzndm98naUn6bu57JP1Q0oymJzMzs7aSUzC+RdEl+S7ArhTfjfGtZoYyM7P2\nk1Mwdo6Ib0XEuvQ4C9ipybnMzKzN5BSM30t6p6RJ6fFO4L5mBzMzs/aSUzDeBfwF8DtgLXBUmmZm\nZpuRUbsGiYj/pei11szMNmM5exhmZmYuGGZmlqflBUPS3pJWVjwelvSBqmX6JD1UsczJrc5pZmbP\nNeo5DEmfjIh/ScNTIuLJ8TQYETcB+6f1TQJ+CyypseilEXHEeNoyM7OJ06h78w9LehXFVVEjJrrn\n2kOB2yLirgler5mZTbBGh6RuAt4O7C7pUkmLgBdI2nsC2z8aWFxn3qskXSPpvyW9fALbNDOzjaCI\nqD1Deh1wJfBr4EDgT4CfARcDe0fEq8fVsLQlsAZ4eUTcUzVvW+CZiBiWNBv4ckTsWWc984B5AN3d\n3T0DAwPjiTVmw8PDdHV1tbTN8eq0zBOZd2htazpenjFlBqufXN2Sthrpmd6TtVyn/U1A52Vu17z9\n/f1DEdGbs2yjgvEZ4GCgFzgLuAaYHxH7TERISXOAv4uIwzKWvRPoTT3o1tXb2xsrVqyYiHjZBgcH\n6evra2mb49VpmScybyu/onX+zfNb0lYjcUrt/+9qnfY3AZ2XuV3zSsouGHUPSUXExyPiUOBO4LsU\nJ8h3knSZpJ9OQM651DkcJemFkpSGD0o53R2JmVmJRr1KCrgwIq4CrpJ0YkS8RtK08TQqaWvgTcDf\nVkx7D0BEnE5xov1ESeuAx4Gjo96ukJmZtURO1yAfrhg9Pk1reGgoY52PAS+omnZ6xfBpwGnjacPM\nzCbWmG7c8zfvmZltvtw1iJmZZXHBMDOzLC4YZmaWxQXDzMyyuGCYmVkWFwwzM8vigmFmZllcMMzM\nLIsLhpmZZXHBMDOzLC4YZmaWxQXDzMyyuGCYmVkWFwwzM8vigmFmZllcMMzMLIsLhpmZZXHBMDOz\nLC4YZmaWpbSCIelOSddJWilpRY35kvQVSbdKulbSAWXkNDOzwuSS2++PiD/UmXc4sGd6HAz8Z/pp\nZmYlaOdDUnOA70RhObC9pOllhzIz21wpIsppWLoDeAAI4OsRsahq/vnAZyPisjT+C+AjEbGiarl5\nwDyA7u7unoGBgVbE/6Ph4WG6urpa2uZ4dVrmicw7tHZoQtYzmhlTZrD6ydUtaWsiNCNvz/SeCV1f\ntXp/F636HdfS6DW36/9df3//UET05ixb5iGpQyJijaSdgaWSboyIZRXzVeM5G1S3VGgWAfT29kZf\nX19TwtYzODhIq9scr07LPJF5+xf0T8h6RrNwr4XMv3l+S9qaCM3IG3Ob+2G03t9Fq37HtTR6zZ32\nf1dLaYekImJN+nkvsAQ4qGqR1cBuFeMzgDWtSWdmZtVKKRiStpE0dWQYOAxYVbXYecBfpaulXgk8\nFBFrWxzVzMySsg5JdQNLJI1k+F5E/FzSewAi4nTgAmA2cCvwGPDXJWU1MzNKKhgRcTuwX43pp1cM\nB/B3rcxlZmb1tfNltWZm1kZcMMzMLIsLhpmZZXHBMDOzLC4YZmaWxQXDzMyyuGCYmVkWFwwzM8vi\ngmFmZllcMMzMLIsLhpmZZXHBMDOzLC4YZmaWxQXDzMyyuGCYmVkWFwwzM8vigmFmZlnK+orWtqMF\n2qjnLdxrIf0L+ic4TXOVnTlOidLattbZ2P+pXGX/HW+OvIdhZmZZXDDMzCxLywuGpN0k/VLSDZKu\nl3RSjWX6JD0kaWV6nNzqnGZm9lxlnMNYB3woIq6WNBUYkrQ0Iv6narlLI+KIEvKZmVkNLd/DiIi1\nEXF1Gn4EuAHYtdU5zMxsbBRR3hUrkmYCy4B9I+Lhiul9wA+B1cAaYH5EXF9nHfOAeQDd3d09AwMD\nG5VlaO3QRj1vxpQZrH5y9UY9tyxlZ+6Z3jOm5YeHh+nq6pqQtjf29zxWZW/jseq0vNCemRv9bU/k\n3/FE6u/vH4qI3pxlSysYkrqAS4BPR8SPquZtCzwTEcOSZgNfjog9R1tnb29vrFixYuPyjOOy2vk3\nz9+o55al7Mxjvax2cHCQvr6+CWm72Zd6jih7G49Vp+WF9szc6G97Iv+OJ5Kk7IJRylVSkp5HsQdx\ndnWxAIiIhyNiOA1fADxP0rQWxzQzswplXCUl4Azghoj4Yp1lXpiWQ9JBFDnva11KMzOrVsZVUocA\nxwLXSVqZpn0ceBFARJwOHAWcKGkd8DhwdJR5ssXMzFpfMCLiMqDhgeSIOA04rTWJzMwsh/uSMjOb\nII0uqmhm31et6p/NXYOYmVkWFwwzM8vigmFmZllcMMzMLIsLhpmZZXHBMDOzLC4YZmaWxQXDzMyy\nuGCYmVkWFwwzM8virkGs5cb6nRTN7FLBzPJ5D8PMzLK4YJiZWRYXDDMzy+KCYWZmWVwwzMwsiwuG\nmZllccEwM7MsLhhmZpallIIhaZakmyTdKumjNeZPkXROmn+FpJmtT2lmZpVaXjAkTQK+ChwO7APM\nlbRP1WInAA9ExEuBLwGfa21KMzOrVsYexkHArRFxe0Q8BQwAc6qWmQN8Ow2fCxwqaWz9SZiZ2YRS\nRLS2QekoYFZEvDuNHwscHBHvq1hmVVpmdRq/LS3zhxrrmwfMS6N7Azc1+SVUmwZskKvNdVrmTssL\nnZe50/JC52Vu17wvjoidchYso/PBWnsK1VUrZ5liYsQiYNF4Q20sSSsiores9jdGp2XutLzQeZk7\nLS90XuZOy1tLGYekVgO7VYzPANbUW0bSZGA74P6WpDMzs5rKKBhXAXtKeomkLYGjgfOqljkPOC4N\nHwVcHK0+dmZmZs/R8kNSEbFO0vuAC4FJwJkRcb2kfwZWRMR5wBnAf0m6lWLP4uhW5xyD0g6HjUOn\nZe60vNB5mTstL3Re5k7Lu4GWn/Q2M7PO5Du9zcwsiwuGmZllccEYA0lnSro33SdSOf3vU1cn10v6\nt7LyVauVN3W5sjI97pS0ssyM1epk3l/S8pR5haSDysxYqU7e/SRdLuk6ST+VtG2ZGatJ2k3SLyXd\nkP5mT0rTd5S0VNIt6ecOZWeFhnnfnsafkdRWl6s2yPx5STdKulbSEknbl511TCLCj8wH8DrgAGBV\nxbR+4P8BU9L4zmXnbJS3av4XgJPLzpmxjS8CDk/Ds4HBsnOOkvcq4PVp+F3Ap8rOWZV5OnBAGp4K\n3EzRTc+/AR9N0z8KfK7srKPk/ROKm3UHgd6yc2ZmPgyYnKZ/rl22ce7DexhjEBHL2PB+kBOBz0bE\nk2mZe1serI46eQFIXa38BbC4paFGUSdzACOf0rdjw/t2SlMn797AsjS8FPjzloYaRUSsjYir0/Aj\nwA3Arjy3S55vA0eWk/C56uWNiBsiotU9O2RpkPmiiFiXFltOcR9ax3DBGL+9gNemXnUvkXRg2YEy\nvRa4JyJuKTtIhg8An5d0N7AQ+FjJeUazCnhLGn47z71Rta2knqD/FLgC6I6ItVC84QE7l5estqq8\nHaFB5ncB/93qPOPhgjF+k4EdgFcC/wh8v0M6SpxLm+1dNHAi8MGI2A34IMV9Ou3sXcDfSRqiOBzx\nVMl5apLUBfwQ+EBEPFx2ntF0Wl6on1nSJ4B1wNllZdsYZfQltalZDfwoioOSV0p6hqKTsd+XG6u+\n1N3K24CesrNkOg44KQ3/APhmiVlGFRE3UhyrRtJewP8tN9GGJD2P4o3s7Ij4UZp8j6TpEbFW0nSg\nbQ6v1snb1upllnQccARwaHrf6Bjewxi/HwNvgD++OWxJe/ZIWemNwI2RegPuAGuA16fhNwBtfRhN\n0s7p5xbAJ4HTy030XGkP+Azghoj4YsWsyi55jgN+0upstTTI27bqZZY0C/gI8JaIeKysfBut7LPu\nnfSgOISzFniaYs/iBIoC8V2K49ZXA28oO2ejvGn6WcB7ys43hm38GmAIuIbiOHBP2TlHyXsSxVUx\nNwOfJfWo0C6PtD0DuBZYmR6zgRcAv6AoyL8Adiw76yh535q2+ZPAPcCFZWfNyHwrcHfFtNPLzjqW\nh7sGMTOzLD4kZWZmWVwwzMwsiwuGmZllccEwM7MsLhhmZpbFBWMzIemFkgYk3SbpfyRdkO4bQdLL\nJV0s6ebUU+k/jdytLumY1LPmtZJ+LWm/JmTbRdK5G/G8b0raZ6LzVKx/uqTzJ2hdRzbKKuk9kv4q\nc13bS3rvROTqNOnvdvvcbSBpS0nL0s2qNk6+rHYzkN78fw18OyJOT9P2p+i2YgXFPSQnRsRFkram\nuDv1/Ij4qqRXU9x89ICkw4FTI+Lgcl5Ja0n6PHBZRIz7BjZJZ1Fs0w0Ko6TJ8WyHdDnrmpnWte94\nczVD+ntTRDzTxDZmkrkNJJ0C3BoRHdUNR1sq+0YQP5r/oLg7elmdeScA36matgdwd41ldwB+W2c9\nhwGXU9y8+AOgK02/E/hMmreCoivwC4HbSDcPAjNJ3YMDLweupLip6VpgT2Ab4GcUN+6tAt6Rlh0k\ndWtN0TfWdWn+5ypyDQOfTs9dTtHBHhSdAq5K0+ttm9t5ttv64ynu6v8pcAfwPuAfgN+k9e5Yse1+\nTnGj4aXAy4BXU/Roe0d6XXuk7J8BLgE+BJwKzE/reClFl/nXpO25R1WuAeDxtK7PA0o/V6VtMLJ9\n+lI75wI3UvRbNPIh8bPA/6RtvLDGaz8V+C/gYoob+f6mYt4/UnThfi2woOJ3eAPwtbRNXly1vgMp\nPrRck36/U9NzLk2v8Wrg1RW5lwFLUsbTgS0q/p6m1dgGXRQ3G16dtsGcirb3Ay4o+/9wU3iUHsCP\nFvyS4f3Al+rM+yJwUo3pDwDbVk2bD3yzxrLT0j/4Nmn8I6Tv2Uj/4Cem4S+lN5mpwE7AvWn6TJ4t\nGP8BHJOGtwSeT9E9+Dcq2tsu/RwEeoFdgP9N65yc3uSOTMsE8OY0/G/AJ9PwdRTdTQNsX+M1vQQY\nqhg/nuIu3ZHsD/FswfsSRedypDetPdPwwcDFafgs4KiK9Q0CX6sYP5VnC8YVwFvT8FbA1lXZ/ri9\n0vifU3SjPgnoTttiOsUb70MUXWhvQVG0XwPsCNzEs8Wj1us/leLN/fnp93t32s6HAYsoitQWwPkU\n3wkyE3gGeGWNdW1JUXwPTOPbpt/T1sBWadqewIo03Ac8AeyeXtPSkW3HswWjehtMJv29pvm3Vry+\nScDvy/4/3BQePq5nonhTreWP0yX182w3HdVeSfHlML9Kpz62pHhzGnFe+nkdxZ7HI8Ajkp6o8Y1j\nlwOfkDSDolPHWyRdByyU9DmKwxCXVj3nQIovVfp9yno2xZvYjyl6ih05DzEEvCkN/wo4S9L3gVqd\n2U1nww4kf1mR/SGKvY2R1/WK1DPpq4EfVHRYPKXGukecUz1B0lSKQrYEICKeaPD8Ea8BFkfEeooO\nBC+h2CYPA1dG6jNMxbcrzqTYI3oC+Kakn/Hs9qn2k4h4HHhc0i+Bg1Jbh1HsRUDxyX5PiiJ1V0Qs\nr7GevYG1EXFVek0PpzzbAKelw6PrKb4qYMSVEXF7Wm5xarfReS4Bn5H0OorCtStF8fxdRKyX9JSk\nqen3ZxsJCIL0AAADTklEQVTJBWPzcD1wVIN5r6ucIGl3YHjkn0vSKyh6iD08Iu6rsQ4BSyNibp02\nnkw/n6kYHhl/zt9gRHxP0hUUPbxeKOndEXGxpB6Kvnj+VdJFEfHPVe3X83Skj5kUb0qTUzvvkXRw\namelpP2rXtvjFJ/ua72O6tcy8jq2AB6MiP0b5Kn0aI1pG9M1fqPnVGZeT/Ftb+tUfM3tocDRFIfX\n3lDjudUfJCK19a8R8fXnBCjOKdR6PSP5an0o+SBFH1D7UWy7yuJYq+1GjqHY8+uJiKcl3clzf39T\nqtZvG8FXSW0eLgamSPqbkQmSDpT0eorj2q+R9MY0/fnAVygO3yDpRRSfwI+NiJvrrH85cIikl6bn\nbD1yBdZYpWJ1e0R8hWLP5BWSdgEei4jvUnyB0gFVT7sCeL2kaZImUZzPuGSUdvaIiCsi4mSK3oWr\nv+ToZopP49nSJ+c7JL09taGKq8oeoTiclbOO1ZKOTOuYki5EqFS9rmXAOyRNkrQTxQeAK+u1kfaE\ntouICyi+nKpegZsjaStJL6A4THQVxfmnd6V1IGnXkd55G7gR2EXpy8UkTU1XLW1HsefxDHAsxaGj\nEQdJeomKHn/fAVw2yjbYjuIQ59Npb/jFFa/3BRSHpJ4eJaeNwgVjM5A+Yb8VeFO6rPZ6imPUa9Ih\nhznAJyXdRHF45SrgtPT0kyl6Mf2apJWSVtRY/+8pjvEvlnQtRQF52UbGfQewKh0+eRnwHeD/UHzX\nyErgE8C/VLW/luJb+H5JOlEco1/Z9HlJ10laRfGGe03VOh8FbhspgmNwDHCCpGso9t7mpOkDwD9K\n+o2kPUZZx7HA+9O2/DXwwqps91Ec/luVruRaQnFu6BqKDwcfjojfNVj/VOD8tP5LKD7p13IlxcUG\nyym+l3xNRFwEfA+4PB0qPJdRCmFEPEXxe/2PtF2WUnz6/xpwnKTlFIejKvdQLqc4Mb+K4mKBJaNs\ng7OB3vT3eQxFkRrRD1zQKKPl8WW1ZnVIeivFIY5Plp2l1SSdSnFYcmEJbfdRXABwxASt70fAx6JN\nv/+7k/gchlkdEbEkHc6wDiVpS+DHLhYTw3sYZmaWxecwzMwsiwuGmZllccEwM7MsLhhmZpbFBcPM\nzLL8f4TRWK7x6j9BAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f115fd78a20>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# the histogram of the data\n",
    "plt.hist(hist_data, 10, normed=False, facecolor='green')\n",
    "\n",
    "plt.xlabel(stage['IndicatorName'].iloc[0])\n",
    "plt.ylabel('# of Years')\n",
    "plt.title('Histogram Example')\n",
    "\n",
    "plt.grid(True)\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "So the USA has many years where it produced between 19-20 metric tons per capita with outliers on either side."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### But how do the USA's numbers relate to those of other countries?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>CountryName</th>\n",
       "      <th>CountryCode</th>\n",
       "      <th>IndicatorName</th>\n",
       "      <th>IndicatorCode</th>\n",
       "      <th>Year</th>\n",
       "      <th>Value</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>5026275</th>\n",
       "      <td>Arab World</td>\n",
       "      <td>ARB</td>\n",
       "      <td>CO2 emissions (metric tons per capita)</td>\n",
       "      <td>EN.ATM.CO2E.PC</td>\n",
       "      <td>2011</td>\n",
       "      <td>4.724500</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5026788</th>\n",
       "      <td>Caribbean small states</td>\n",
       "      <td>CSS</td>\n",
       "      <td>CO2 emissions (metric tons per capita)</td>\n",
       "      <td>EN.ATM.CO2E.PC</td>\n",
       "      <td>2011</td>\n",
       "      <td>9.692960</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5027295</th>\n",
       "      <td>Central Europe and the Baltics</td>\n",
       "      <td>CEB</td>\n",
       "      <td>CO2 emissions (metric tons per capita)</td>\n",
       "      <td>EN.ATM.CO2E.PC</td>\n",
       "      <td>2011</td>\n",
       "      <td>6.911131</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5027870</th>\n",
       "      <td>East Asia &amp; Pacific (all income levels)</td>\n",
       "      <td>EAS</td>\n",
       "      <td>CO2 emissions (metric tons per capita)</td>\n",
       "      <td>EN.ATM.CO2E.PC</td>\n",
       "      <td>2011</td>\n",
       "      <td>5.859548</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5028456</th>\n",
       "      <td>East Asia &amp; Pacific (developing only)</td>\n",
       "      <td>EAP</td>\n",
       "      <td>CO2 emissions (metric tons per capita)</td>\n",
       "      <td>EN.ATM.CO2E.PC</td>\n",
       "      <td>2011</td>\n",
       "      <td>5.302499</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                     CountryName CountryCode  \\\n",
       "5026275                               Arab World         ARB   \n",
       "5026788                   Caribbean small states         CSS   \n",
       "5027295           Central Europe and the Baltics         CEB   \n",
       "5027870  East Asia & Pacific (all income levels)         EAS   \n",
       "5028456    East Asia & Pacific (developing only)         EAP   \n",
       "\n",
       "                                  IndicatorName   IndicatorCode  Year  \\\n",
       "5026275  CO2 emissions (metric tons per capita)  EN.ATM.CO2E.PC  2011   \n",
       "5026788  CO2 emissions (metric tons per capita)  EN.ATM.CO2E.PC  2011   \n",
       "5027295  CO2 emissions (metric tons per capita)  EN.ATM.CO2E.PC  2011   \n",
       "5027870  CO2 emissions (metric tons per capita)  EN.ATM.CO2E.PC  2011   \n",
       "5028456  CO2 emissions (metric tons per capita)  EN.ATM.CO2E.PC  2011   \n",
       "\n",
       "            Value  \n",
       "5026275  4.724500  \n",
       "5026788  9.692960  \n",
       "5027295  6.911131  \n",
       "5027870  5.859548  \n",
       "5028456  5.302499  "
      ]
     },
     "execution_count": 26,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# select CO2 emissions for all countries in 2011\n",
    "hist_indicator = 'CO2 emissions \\(metric'\n",
    "hist_year = 2011\n",
    "\n",
    "mask1 = data['IndicatorName'].str.contains(hist_indicator) \n",
    "mask2 = data['Year'].isin([hist_year])\n",
    "\n",
    "# apply our mask\n",
    "co2_2011 = data[mask1 & mask2]\n",
    "co2_2011.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "For how many countries do we have CO2 per capita emissions data in 2011"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "232\n"
     ]
    }
   ],
   "source": [
    "print(len(co2_2011))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAEWCAYAAAB8LwAVAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XucVXW9//HXW/AOiooRCkcI0TJNkvGCWg3eDnlUtFMp\nmWHp4dftaBcqNU/o6Z5Y2elKWmoSeL+kZRo6oiUqkHdTUFFJBAUV8ZaXz++P73dkM7NmZjMze/Zm\n9vv5eOzHrL1u38/6zt77s9Z3rfVdigjMzMxaWq/aAZiZWW1ygjAzs0JOEGZmVsgJwszMCjlBmJlZ\nIScIMzMr5ASxjpB0n6TGasdRTZKOkPSEpFWS3lvteLqbpFMknd2F5Y+WdF13xtSb5c/RO6odRy1z\ngqgBkhZJOqDFuGMl3dL8PiLeHRFNHaxnmKSQ1LdCoVbbVODzEdEvIv7ecqKSEyTdK+lFSYslXSxp\nl5J59pZ0g6QXJD0v6Q+SdiqZvpek6yWtkPR0Xn5wWwFJapL0Sv6xaX79oTMbFxHfiYjjO7NsXn56\nRBzU2eU7I392X87bvVTSbyX168b175D/B8/k/9fdkr4kqU9X150/R4/kcs6V9K2uR9y7OEFY2Wog\n8WwH3NfO9LOAE4ETgC2BHYArgP8AkDQGuA64EtgGGA7cBfy1ZE9yC2AaMCyX9wLw2w7iak5aza9D\n137T1mmHRkQ/YDdgd+DUtV1B0WdL0gjgNuAJYJeI2Bz4CNAA9O9SxFaeiPCryi9gEXBAi3HHArcU\nzQPsAcwFVgJLgR/m8Y8DAazKrzGknYBTgceAZcD5wOYl6/1EnrYc+J8W5ZwGXAJckMs6Ppd9K/Ac\nsAT4KbBByfoC+CywgPTj+k1gRF5mJXBR6fwttrkwVmDDvD0BvAg8XLDsSOANYI926vlm4OcF4/8E\nnN/GMrsBL7Szzibg+DamNQKLga/m7VkCHA4cDDwErABOKZn/NOCCPLxRrvflua7vAAaVfDYeyfX7\nKHB0G5+ZvfNyz+e/e7eI+5vAX/N6rgMGdlR2R59d4Azg6jy8OXBO3u5/At8C+pTE+lfgR7kevlWw\n7guAazr47lwMPJW3cTbw7pJp5wK/BK7P23gTsF2Lz+r2wCTgNeBf+XP2hzz9JODhvOz9wBHV/q3o\n6ZePINZNZwFnRcRmpB/fi/L49+e/AyLtyd5K+iIeC4wF3gH0I/2ok5tWfg4cDQwmfaG3bVHWeFKS\nGABMJ/0IfxEYSEpA+5MSQqlxwGhgL9KP47RcxlBgZ2BCG9tVGGtEvBppDxVg14gYUbDs/sDiiLi9\naMWSNiH9YF5cMPki4MA2Yno/7R+1dOTtpB/cbYFvAL8GPk6qn/cB32ijHXwi6f8xFNgK+DTwsqRN\ngZ8AH4yI/qRturPlwpK2BK7J824F/BC4RtJWJbN9DPgk8DZgA2Bye2V3tKGShpKSX3Pz33nA66Qf\n4fcCB5F2MprtSUp0bwO+XbDKA0ifvfb8ibRz8DZgPukzWupoUiIcSKqnltOJiGl5/A9izSPAh0n/\no82B04EL2mtu7JWqnaH8emsvbBVpb6359RJtH0HMJn1gB7ZYzzDSXlHfknGzgM+WvN+RtLfUl/SD\nNaNk2iakvajSI4jZHcT+BeDykvcB7FPyfh7wtZL3ZwI/bmNdbcZasu7t21j268CcduIckpd/Z8G0\nccBrBePfQ9q7fV87623K/6vS/90387RG0g9r815z/xzDni3q5/CS+m4+gvgU8DfgPS3K2zSX8Z/A\nxi2mHdv8mQGOAW5vMf1W4NiSuE8tmfZZ4Nr2yi7js/sYaYdjY2AQ8GppjKQdgxtLYn28g3W/Boxb\ni+/RgFy/m+f35wIzS6b3I+3gDG35ecrztjqKabH+O4Hx5cbTG14+gqgdh0fEgOYXrffKSx1Hal//\nh6Q7JB3SzrzbkL64zR4jJYdBedoTzRMi4iVSs0KpJ0rf5JOGV0t6StJK4DukvbNSS0uGXy5439ZJ\nzPZi7chy0lFQW54F3mxjnsHAM6UjJG1P2js9MSJu7qDsE0r/dxHxP6VxRcQbebh5L7yc+vgd8Gdg\npqQnJf1A0voR8SJwJGmvfomkayS9s2D5lnVJfl96hPhUyfBLJXEUll245UnzZ3e7iPhsRLxMOn+z\nfo7xOUnPAb8i7ek3e6JoZSXa/Z9K6iPpe5Iezp/FRXlS6eex9PO9ipTwt+mg3Ob1f0LSnSXx70zr\nz3qv5gSxDoqIBRExgfRl+z5wSW56KOqa90nSl7XZv5EO+5eS2oaHNE+QtDGpSWGN4lq8/wXwD2Bk\npCauUwB1fmvKjrUjs4AhkhqKJuYf1ltJJzlb+mheHgBJ2wF/IR0J/K680LtXRLwWEadHxE6kZqRD\nSOeLiIg/R8SBpB/Pf5CarVpqWZeQ6vOfXSl7LTxBOoIYWJI4N4uId5cW1cE6/kI6UmrLx0hNoAeQ\nmoGG5fGln8ehzQP56qotSXXT0hqx5M/Ar4HPA1vlnbZ76b7P+jrBCWIdJOnjkraOiDdJh/aQDp2f\nJu0ll7ZpzwC+KGl4/oJ8B7gwIl4nte8emi/93IDUbNXRF6A/6WTzqrzn+plu27D2Y21XRCwgNW/M\nkNQoaQNJG0k6StJJebaTgIn5Utj+krbIlzaOIW07krYFbgB+FhG/7MZtWyuSxkraJV/OuZLU3PKG\npEGSDss7BK+SmnfeKFjFH4EdJH1MUl9JRwI7AVd3tuy1iT8ilpBOfJ8paTNJ60kaIekDa7GaKcDe\nks6Q9PYc2/aSLpA0gPRZfJV0pLEJ6fPS0sGS9s2f728Ct0VE0ZHLUtb83jTvcD2dy/0k6QiirjhB\nrJvGAfdJWkU6YX1URLySm4i+Tbps8zlJewG/ITUZzCZd8fIK8N8AEXFfHp5JOpp4gXS1zavtlD2Z\ntOf2AmkP68Ju3K42Yy3TCaQT8D8jJc6HgSOAPwBExC3AvwMfIm3vY6STp/vmBAPpJOo7gCml9zZ0\nUO5PW9wHMW8tYm7L20kJfCXwAOkKnAtI39kvk/aCVwAfoKA5MiKWk/b8v0z6Af0qcEhEPNNy3rUo\ne219gnTy+35SE98ltN8MuIaIeJiUvIeRPu/PA5eSruB7gXSV22Oko6L7gTkFq/k9KdGsIF0YcHQb\nxZ0D7JS/N1dExP2k82W3kpLHLqSrruqK8skXs+ZD8OdIzUePVjses66QdC7pyra1vi/DEh9B1DlJ\nh0raJDdZTAXuYfXJPjOrY04QNp7UXPEk6Xryo8KHlWaGm5jMzKwNPoIwM7NC1e58rUsGDhwYw4YN\n69SyL774Iptuumn3BrSOc50Uc7205jppbV2qk3nz5j0TEVt3NF/FEoSk35Aus1sWETu3mDaZ1KnX\n1hHxjCSRLtc8mHRH57ERMb+jMoYNG8bcuXM7FV9TUxONjY2dWra3cp0Uc7205jppbV2qE0kt77Iv\nVMkmpnNJ1+uvIXfodSCp59FmHySdIB1J6lnxFxWMy8zMylCxBBERs0k3p7T0I9JNO6Vnx8eTuluO\niJgDDKi7XhPNzGpMj56klnQY8M+IuKvFpG1Zs+OuxbTudtrMzHpQj52kzv3xf53UJ3yryQXjCq+/\nlTSJ1AzFoEGDaGpq6lQ8q1at6vSyvZXrpJjrpTXXSWu9sU568iqmEeRHPKZz0gwB5kvag3TEMLRk\n3iEU97hIpId7TANoaGiIzp4UWpdOKPUU10kx10trrpPWemOd9FgTU0TcExFvi4hhETGMlBR2i4in\ngKuATyjZC3g+9wZpZmZVUrEEIWkGqSfEHSUtlnRcO7P/kfTowYWkHkLbe1iOmZn1gIo1MeUH2rQ3\nfVjJcACfq1QsZma29tzVhpmZFVqnu9roinlL5jH29LFVKTumuINEM6t9PoIwM7NCThBmZlbICcLM\nzAo5QZiZWSEnCDMzK+QEYWZmhZwgzMyskBOEmZkVcoIwM7NCThBmZlbICcLMzAo5QZiZWSEnCDMz\nK+QEYWZmhZwgzMyskBOEmZkVcoIwM7NCThBmZlbICcLMzAo5QZiZWaGKJQhJv5G0TNK9JePOkPQP\nSXdLulzSgJJpJ0taKOlBSf9eqbjMzKw8lTyCOBcY12Lc9cDOEfEe4CHgZABJOwFHAe/Oy/xcUp8K\nxmZmZh2oWIKIiNnAihbjrouI1/PbOcCQPDwemBkRr0bEo8BCYI9KxWZmZh3rW8WyPwVcmIe3JSWM\nZovzuFYkTQImAQwaNIimpqZOFT5kwyFM3WFqp5btqs7GXGmrVq2q2diqyfXSmuuktd5YJ1VJEJK+\nDrwOTG8eVTBbFC0bEdOAaQANDQ3R2NjYqRjOnHEmkx+a3KlluyomFG5a1TU1NdHZ+uzNXC+tuU5a\n64110uMJQtJE4BBg/4ho/qVcDAwtmW0I8GRPx2ZmZqv16GWuksYBXwMOi4iXSiZdBRwlaUNJw4GR\nwO09GZuZma2pYkcQkmYAjcBASYuBKaSrljYErpcEMCciPh0R90m6CLif1PT0uYh4o1KxmZlZxyqW\nICJiQsHoc9qZ/9vAtysVj5mZrR3fSW1mZoWcIMzMrJAThJmZFXKCMDOzQk4QZmZWyAnCzMwKOUGY\nmVkhJwgzMyvkBGFmZoWcIMzMrJAThJmZFXKCMDOzQk4QZmZWyAnCzMwKOUGYmVkhJwgzMyvkBGFm\nZoWcIMzMrJAThJmZFXKCMDOzQk4QZmZWqGIJQtJvJC2TdG/JuC0lXS9pQf67RR4vST+RtFDS3ZJ2\nq1RcZmZWnkoeQZwLjGsx7iRgVkSMBGbl9wAfBEbm1yTgFxWMy8zMylCxBBERs4EVLUaPB87Lw+cB\nh5eMPz+SOcAASYMrFZuZmXWsp89BDIqIJQD579vy+G2BJ0rmW5zHmZlZlfStdgCZCsZF4YzSJFIz\nFIMGDaKpqalTBQ7ZcAhTd5jaqWW7qrMxV9qqVatqNrZqcr205jpprTfWSU8niKWSBkfEktyEtCyP\nXwwMLZlvCPBk0QoiYhowDaChoSEaGxs7FciZM85k8kOTO7VsV8WEwtxXdU1NTXS2Pnsz10trrpPW\nemOd9HQT01XAxDw8EbiyZPwn8tVMewHPNzdFmZlZdVTsCELSDKARGChpMTAF+B5wkaTjgMeBj+TZ\n/wgcDCwEXgI+Wam4zMysPBVLEBExoY1J+xfMG8DnKhWLmZmtPd9JbWZmhZwgzMyskBOEmZkVcoIw\nM7NCHSYISftI2jQPf1zSDyVtV/nQzMysmso5gvgF8JKkXYGvAo8B51c0KjMzq7pyEsTr+TLU8cBZ\nEXEW0L+yYZmZWbWVcx/EC5JOBo4B3iepD7B+ZcMyM7NqK+cI4kjgVeBTEfEUqZfVMyoalZmZVV2H\nCSInhUuBDfOoZ4DLKxmUmZlVXzlXMf0XcAnwqzxqW+CKSgZlZmbVV04T0+eAfYCVABGxgNUP+jEz\ns16qnATxakT8q/mNpL608TAfMzPrPcpJEDdJOgXYWNKBwMXAHyoblpmZVVs5CeIk4GngHuD/kZ7d\ncGolgzIzs+rr8D6IiHgT+HV+mZlZnWgzQUi6KCI+KukeCs45RMR7KhqZmZlVVXtHECfmv4f0RCBm\nZlZb2kwQEbEkd6txTkQc0IMxmZlZDWj3JHVEvEHqyXXzHorHzMxqRDmd9b0C3CPpeuDF5pERcULF\nojIzs6orJ0Fck1+lfKOcmVkvV06CGJCfAfEWSSe2NXM5JH0ROJ6UaO4BPgkMBmYCWwLzgWNK7+A2\nM7OeVc6NchMLxh3b2QIlbQucADRExM5AH+Ao4PvAjyJiJPAscFxnyzAzs65r7z6ICcDHgOGSriqZ\n1B9Y3g3lbizpNWATYAmwXy4P4DzgNNLjTs3MrAqUniZaMEHaDhgOfJfU3UazF4C7I+L1Theamqi+\nDbwMXEe652JORGyfpw8F/pSPMFouOwmYBDBo0KDRM2fO7FQMS1csZfGrizu3AV00evDoqpTbkVWr\nVtGvX79qh1FzXC+tuU5aW5fqZOzYsfMioqGj+dq7D+Ix4DFgTHcGJmkL0vOthwPPkTr/+2BRCG3E\nNQ2YBtDQ0BCNjY2diuPMGWcy+aHJnVq2q2JCbZ7jb2pqorP12Zu5XlpznbTWG+uknAcGfUjSAknP\nS1op6QVJK7tQ5gHAoxHxdES8BlwG7A0MyF2JAwwBnuxCGWZm1kXlnKT+AXBYRGweEZtFRP+I2KwL\nZT4O7CVpE0kC9gfuB24EPpznmQhc2YUyzMysi8pJEEsj4oHuKjAibiM9wnQ+6RLX9UhNRl8DviRp\nIbAVcE53lWlmZmuvnPsg5kq6kPQc6lebR0bEZZ0tNCKmAFNajH4E2KOz6zQzs+5VToLYDHgJOKhk\nXJDOHZiZWS9VzgODPtkTgZiZWW3pMEFI+i3FDwz6VEUiMjOzmlBOE9PVJcMbAUfgS1DNzHq9cpqY\nLi19L2kG8JeKRWRmZjWhnMtcWxoJ/Ft3B2JmZrWlnHMQL5DOQSj/fYp0z4KZmfVi5TQx9e+JQMzM\nrLaUc5IaSYcB789vmyLi6vbmNzOzdV85nfV9j9Qd9/35daKk71Y6MDMzq65yjiAOBkZFxJsAks4D\n/g6cXMnAzMysusq9imlAyfDmlQjEzMxqSzlHEN8F/i7pRtKVTO/HRw9mZr1eOVcxzZDUBOxOShBf\ni4inKh2YmZlVV5sJQtK/A/0j4pKIWAJclccfLWlZRFzfU0GamVnPa+8cxOnATQXjZwH/W5lwzMys\nVrSXIDaJiKdbjszNS5tWLiQzM6sF7SWIjSS1aoKStD6wceVCMjOzWtBegrgM+LWkt44W8vAv8dPk\nzMx6vfYSxKnAUuAxSfMkzQMWAU/naWZm1ou1eRVTRLwOnCTpdGD7PHphRLzcI5GZmVlVdXgndUS8\nHBH35Fe3JAdJAyRdIukfkh6QNEbSlpKul7Qg/92iO8oyM7PO6cwDg7rDWcC1EfFOYFfgAeAkYFZE\njCRdSntSlWIzMzPaSRCS9sl/N+zOAiVtRuqu4xyAiPhXRDwHjAfOy7OdBxzeneWamdnaUUQUT5Dm\nRcRoSfMjYrduK1AaBUwjdR2+KzCP1J34PyNiQMl8z0ZEq2YmSZOASQCDBg0aPXPmzE7FsXTFUha/\nurhTy3bV6MGjq1JuR1atWkW/fv2qHUbNcb205jppbV2qk7Fjx86LiIaO5msvQcwhNf0cDFzYcnpE\nnNCZwCQ1AHOAfSLiNklnASuB/y4nQZRqaGiIuXPndiYMzpxxJpMfmtypZbsqphTXebU1NTXR2NhY\n7TBqjuulNddJa+tSneQDgA4TRHud9R0CHADsR9rL7y6LgcURcVt+fwnpfMNSSYMjYomkwcCybizT\nzMzWUnuXuT4DzJT0QETc1V0FRsRTkp6QtGNEPAjsz+qn1U0Evpf/XtldZZqZ2dor53kQyyVdDuwD\nBHALcGJEdKUB/7+B6ZI2AB4BPkk6YX6RpOOAx4GPdGH9ZmbWReUkiN8Cv2f1D/bH87gDO1toRNwJ\nFLV/7d/Zda5LdLqqUm6tnvsws9pUzn0Qb4uI30bE6/l1LrB1heMyM7MqKydBPC3p45L65NfHgeWV\nDszMzKqrnATxKeCjwFPAEuDDeZyZmfVi5TyT+nHgsB6IxczMaki1+mIyM7Ma5wRhZmaFnCDMzKxQ\nhwlC0qklw93as6uZmdWu9rr7/qqkMaSrlprdWvmQzMysFrR3FdODpLun3yHpZlLPrluV9KFkZma9\nWHtNTM8CpwALgUbgJ3n8SZL+VuG4zMysyto7ghgHTAFGAD8E7gJejIhP9kRgZmZWXW0eQUTEKRGx\nP7AIuICUTLaWdIukP/RQfGZmViXl9Ob654i4A7hD0mciYl9JAysdmJmZVVeHl7lGxFdL3h6bxz1T\nqYDMzKw2rNWNct35ZDkzM6ttvpPazMwKOUGYmVkhJwgzMyvkBGFmZoWcIMzMrJAThJmZFapagpDU\nR9LfJV2d3w+XdJukBZIulLRBtWIzM7PqHkGcSOohttn3gR9FxEhSR4HHVSUqMzMDqpQgJA0B/gM4\nO78XsB9wSZ7lPODwasRmZmaJIqLnC5UuAb4L9Acmk7rwmBMR2+fpQ4E/RcTOBctOAiYBDBo0aPTM\nmTM7FcPSFUtZ/OriTi27rho9eHS701etWkW/fv16KJp1h+ulNddJa+tSnYwdO3ZeRDR0NF85nfV1\nK0mHAMsiYp6kxubRBbMWZq6ImAZMA2hoaIjGxsai2Tp05owzmfzQ5E4tu66KCe3vDDQ1NdHZ+uzN\nXC+tuU5a64110uMJAtgHOEzSwcBGwGbAj4EBkvpGxOvAEODJKsRmZmZZj5+DiIiTI2JIRAwDjgJu\niIijgRtZ/fzricCVPR2bmZmtVkv3QXwN+JKkhcBWwDlVjsfMrK5Vo4npLRHRBDTl4UeAPaoZj5mZ\nrVZLRxBmZlZDnCDMzKyQE4TVlUWLFrHzzmveXnPaaacxdepU5syZw5577smoUaN417vexWmnnbbG\nfOPHj2fMmDE9GK1ZdVX1HIRZLZk4cSIXXXQRu+66K2+88QYPPvjgW9Oee+455s+fT79+/Xj00UcZ\nPnx4FSM16xk+gjDLli1bxuDBgwHo06cPO+2001vTLr30Ug499FCOOuooOnv3vtm6xgnCLPviF7/I\njjvuyBFHHMGvfvUrXnnllbemzZgxgwkTJjBhwgRmzJhRxSjNeo4ThNWV1C9k8fhvfOMbzJ07l4MO\nOojf//73jBs3DoAVK1awcOFC9t13X3bYYQf69u3Lvffe25Nhm1WFE4TVla222opnn312jXErVqxg\n4MCBAIwYMYLPfOYzzJo1i7vuuovly5dz44038uyzzzJ8+HCGDRvGokWL3MxkdcEJwupKv379GDx4\nMLNmzQJScrj22mvZd999ueaaa2ju3XjBggX06dOHAQMGcMMNN3DttdeyaNEiFi1axLx585wgrC44\nQVjdOf/88/nWt77FqFGj2G+//ZgyZQojRozgd7/7HTvuuCOjRo3imGOOYfr06TzxxBMsW7aMvfba\n663lhw8fzmabbcZtt91Wxa0wqzxf5mp1Z6edduLGG29sNb6to4KLL7641bmL+fPnVyQ2s1riIwgz\nMyvkBGFmZoWcIMzMrJAThBlwxx138Nxzz1U7DLOa4gRhBhx//PEsWLCg2mGY1RQnCKt7y5cvZ9Gi\nRbz3ve+tdihmNcUJwure7Nmz2Xvvvenb11d9m5VygrC619TURGNjY7XDMKs5ThBW95wgzIo5QVhd\nW758OY8++ii77bZbtUMxqzk93ugqaShwPvB24E1gWkScJWlL4EJgGLAI+GhEPNvWemzt6fTirq6b\nTd1hKmNPH1uRsmNKVGS9XTV79mz22Wcf1l9//WqHYlZzqnEE8Trw5Yh4F7AX8DlJOwEnAbMiYiQw\nK783qyg3L5m1rccTREQsiYj5efgF4AFgW2A8cF6e7Tzg8J6OzerPTTfd5ARh1gY1939flcKlYcBs\nYGfg8YgYUDLt2YjYomCZScAkgEGDBo3ubL/8S1csZfGrizu1bG81ZMMhFauT0YNHV2S9XbFy5Uo+\n9rGPccUVV7R7ieuqVavo169fD0ZW+1wnra1LdTJ27Nh5EdHQ0XxVu/BbUj/gUuALEbGyrUdBthQR\n04BpAA0NDdHZvb8zZ5zJ5Icmd2rZ3mrqDlMrVicxofbOQVxxxRW8733v44ADDmh3PjdDteY6aa03\n1klVrmKStD4pOUyPiMvy6KWSBufpg4Fl1YjN6kdv/EKbdaceTxBKhwrnAA9ExA9LJl0FTMzDE4Er\nezo2qy9NTU184AMfqHYYZjWrGk1M+wDHAPdIujOPOwX4HnCRpOOAx4GPVCE2qxMrVqzgkUceYfTo\n2js3YlYrejxBRMQtQFsnHPbvyVisfjX3v+T7H8za5juprS65ecmsY04QVpeampoYO7Yyd42b9RZO\nEFZ3fP7BrDxOEFZ3fP7BrDxOEFZ3fP+DWXmcIKzuuP8ls/I4QVhdWbFiBQ8//LDPP5iVwQnC6srN\nN9/MmDFjfP7BrAxOEFZXfP7BrHxOEFYXpkyZwssvv+wEYbYWnCCsLlx77bXMmjWLhx9+mG222Yav\nfOUr1Q7JrOY5QVhdGDNmDBdffDG77747H/rQhxgwYEDHC5nVOScIqwtjxozhlltuYcmSJeyyyy6c\ncsop1Q7JrOY5QVhdGDNmDI888ggbb7wx06ZNo9wnGJrVs6o9ctTqi06vzg9yTEmPOh06dCj7778/\n06dPZ4MNNqhKLGbrGicIqwuS+Mtf/lLtMMzWKW5iMjOzQj6CsF6tO5q2pu4wlbGnr/2zI5qbt8zW\nVT6CMDOzQk4QZmZWyAnCzMwKOUGYmVmhmjtJLWkccBbQBzg7Ir5X5ZDMOqXa936YdVVNJQhJfYCf\nAQcCi4E7JF0VEfdXNzKzdUdPJKaiK7uqmZiqlYxLdfZqt87qifqutSamPYCFEfFIRPwLmAmMr3JM\nZmZ1SRG1czgq6cPAuIg4Pr8/BtgzIj5fMs8kYFJ+uyPwYCeLGwg804VweyPXSTHXS2uuk9bWpTrZ\nLiK27mimmmpiAoqOE9fIYBExDZjW5YKkuRHR0NX19Cauk2Kul9ZcJ631xjqptSamxcDQkvdDgCer\nFIuZWV2rtQRxBzBS0nBJGwBHAVdVOSYzs7pUU01MEfG6pM8DfyZd5vqbiLivQsV1uZmqF3KdFHO9\ntOY6aa3X1UlNnaQ2M7PaUWtNTGZmViOcIMzMrFBdJghJ4yQ9KGmhpJOqHU81SPqNpGWS7i0Zt6Wk\n6yUtyH+3qGaMPU3SUEk3SnpA0n2STszj67ZeJG0k6XZJd+U6OT2PHy7ptlwnF+aLSuqKpD6S/i7p\n6vy+19VJ3SWIku48PgjsBEyQtFN1o6qKc4FxLcadBMyKiJHArPy+nrwOfDki3gXsBXwufzbquV5e\nBfaLiF2BUcA4SXsB3wd+lOvkWeC4KsZYLScCD5S873V1UncJAnfnAUBEzAZWtBg9HjgvD58HHN6j\nQVVZRCyJiPl5+AXSl39b6rheIlmV366fXwHsB1ySx9dVnQBIGgL8B3B2fi96YZ3UY4LYFnii5P3i\nPM5gUESrpvFoAAAJDUlEQVQsgfRjCbytyvFUjaRhwHuB26jzeslNKXcCy4DrgYeB5yLi9TxLPX6H\nfgx8FXgzv9+KXlgn9ZggOuzOw+qbpH7ApcAXImJlteOptoh4IyJGkXo22AN4V9FsPRtV9Ug6BFgW\nEfNKRxfMus7XSU3dKNdD3J1H25ZKGhwRSyQNJu0x1hVJ65OSw/SIuCyPrvt6AYiI5yQ1kc7PDJDU\nN+8x19t3aB/gMEkHAxsBm5GOKHpdndTjEYS782jbVcDEPDwRuLKKsfS43I58DvBARPywZFLd1ouk\nrSUNyMMbAweQzs3cCHw4z1ZXdRIRJ0fEkIgYRvr9uCEijqYX1kld3kmdM/+PWd2dx7erHFKPkzQD\naCR1UbwUmAJcAVwE/BvwOPCRiGh5IrvXkrQvcDNwD6vblk8hnYeoy3qR9B7SCdc+pB3KiyLifyW9\ng3SBx5bA34GPR8Sr1Yu0OiQ1ApMj4pDeWCd1mSDMzKxj9djEZGZmZXCCMDOzQk4QZmZWyAnCzMwK\nOUGYmVkhJ4heStLbJc2U9LCk+yX9UdIOedq7Jd0g6aHc8+T/5HsAkHS0pLvz62+Sdq1AbNtIuqTj\nOVstd3YlO1aUNLi5Z85uWNfh7cUq6dOSPlHmugZI+mx3xLWuyZ/bAeXWgaQNJM2WVI83AXc7X+ba\nC+Uf+78B50XEL/O4UUB/YC5wL/CZiLhO0iakO4evjoifSdqbdKPYs5I+CJwWEXtWZ0t6lqQzgFsi\noss3OEk6l1SnrRJhyd225a5rWF7Xzl2NqxLy500R8WaHM3e+jGGUWQeSppA65JxeqXjqRkT41cte\npF4lZ7cx7Tjg/BbjRgBPFMy7BfDPNtZzEHArMB+4GOiXxy8CvpOnzQV2Iz1j/GHg03meYcC9efjd\nwO3AncDdwEhgU+Aa4C5SMjsyz9sENOThCaQb2u4Fvl8S1yrg23nZOaSO9gA+kue9q526eQTYMA8f\nS7px8A/Ao8DngS+RboCaA2xZUnfXAvNIN9m9E9ib1FPuo3m7RuTYvwPcBHwZOI10gxXA9sBfcmzz\ngREt4poJvJzXdQap358z8vbcU1I/jbmcS4B/ANNZvRP4PeD+XMdTC7b9NOB3wA3AAuC/SqZ9hdQD\nwd3A6SX/wweAn+c62a7F+nYn7aTclf+//fMyN+dtnA/sXRL3bODyHOMvgfVKPk8DC+qgH6nr9fm5\nDsaXlL0r8Mdqfw97w6vqAfhVgX8qnEDql75o2g+BEwvGPwts1mLcZODsgnkH5i/0pvn914Bv5OFF\npKMTgB/lH5X+wNakDs6af1yaE8T/AUfn4Q2AjYH/BH5dUt7m+W8T0ABsQ7qjeWtSf2I3AIfneQI4\nNA//ADg1D98DbJuHBxRs03BgXsn7Y4GFJbE/z+oE9yNSR37kH6mReXhPUrcLkJ638eGS9TUBPy95\nfxqrE8RtwBF5eCNgkxaxvVVf+f1/knpV7QMMynUxmPRD+zypH6D1SEl6X9KdvQ+yOlkUbf9ppB/z\njfP/94lczwcB00hJaT3gauD9OaY3gb0K1rUBKdnunt9vlv9PmwAb5XEjgbl5uBF4BXhH3qbrm+uO\n1QmiZR30JX9e8/SFJdvXB3i62t/D3vByO139EW33MvnWeEljSUcb+xbMtxfpYUt/zacuNiD9GDVr\n7tvqHtKRxQvAC5Jeae7Xp8StwNdz//qXRcQCSfcAUyV9n9SscHOLZXYHmiLi6RzrdNKP1hXAv0g/\nYpD26g/Mw38FzpV0EXAZrQ0Gnm4x7saS2J8nHU00b9d7cq+vewMX53oA2LBg3c0ubDlCUn9S4roc\nICJeaWf5ZvsCMyLiDVJHgjeR6mQlcHtELM7rvpP0wzqH9AN8tqRrWF0/LV0ZES8DL0u6kdRz676k\nJPH3PE8/0o/748BjETGnYD07Aksi4o68TStzPJsCP83NnW8AO5Qsc3tEPJLnm5HLbe88lYDvSHo/\nKVFtS0qWT0XEG5L+Jal//v9ZJzlB9E73sbrTsKJp7y8dkfuQWdX8Zcr975wNfDAilhesQ8D1ETGh\njTKa+595s2S4+f0an7mI+L2k20gPX/mzpOMj4gZJo4GDge9Kui4i/rdF+W15LfJuJOlHqG8u59OS\n9szl3ClpVItte5m09160HS23pXk71iM9A2BUO/GUerFgXHvb0pb2limN+Q2gb0S8LmkPYH9S53Kf\nJzVDttRyxyFyWd+NiF+tEUA6J1C0Pc3xFe2EfJHU79eupLorTYZFZbfnaNKR3eiIeE3SItb8/23Y\nYv3WCb6KqXe6AdhQ0n81j5C0u6QPkNql95V0QB6/MfATUnMMkv6NtId9TEQ81Mb65wD7SNo+L7NJ\n8xVSaysnp0ci4iekI4/3SNoGeCkiLgCmks5jlLoN+ICkgUqPkJ1Aattvr5wREXFbRHwDeIY1u3wH\neIi0t122vGf8qKSP5DJUctXXC6TmqXLWsVjS4XkdG+YLB0q1XNds4EilB/lsTUr4t7dVRj7S2Twi\n/gh8gfTo0CLjlZ5BvRWp2ecO0vmjT+V1IGlbSR09MOkfwDaSds/L9M9XFW1OOrJ4EziG1BTUbA+l\nHpbXA44EbumgDjYnNVm+lo92tyvZ3q1ITUyvdRCndcAJohfKe9BHAAfmy1zvI7UxP5mbEMYDp0p6\nkNRccgfw07z4N0hPx/q5pDslzS1Y/9OkNvoZku4mJYx3djLcI4F7c3PIO4HzgV2A2/O4rwPfalH+\nEuBkUvfKdwHzo+Mrj86QdI+ke0k/sHe1WOeLwMPNSW8tHA0cJ+ku0tFZ8+NrZwJfUXqo/YgO1nEM\ncEKuy78Bb28R23JSc969+Uqry0nndu4i7Qx8NSKeamf9/YGr8/pvIu3JF7mddHHAHOCbEfFkRFwH\n/B64NTf9XUIHiS/So3yPBP4v18v1pL37nwMTJc0hNS+VHoHcSjqRfi/p5P7lHdTBdKAhfz6PJiWl\nZmOBP7YXo5XHl7maZZKOIDVZnFrtWHqapNNIzYxTq1B2I7nL7G5a32XAyRHxYHesr575HIRZFhGX\n5+YJW0cpPQTsCieH7uEjCDMzK+RzEGZmVsgJwszMCjlBmJlZIScIMzMr5ARhZmaF/j8l/kJ9o7Na\nSgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f115fd21978>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# let's plot a histogram of the emmissions per capita by country\n",
    "\n",
    "# subplots returns a touple with the figure, axis attributes.\n",
    "fig, ax = plt.subplots()\n",
    "\n",
    "ax.annotate(\"USA\",\n",
    "            xy=(18, 5), xycoords='data',\n",
    "            xytext=(18, 30), textcoords='data',\n",
    "            arrowprops=dict(arrowstyle=\"->\",\n",
    "                            connectionstyle=\"arc3\"),\n",
    "            )\n",
    "\n",
    "plt.hist(co2_2011['Value'], 10, normed=False, facecolor='green')\n",
    "\n",
    "plt.xlabel(stage['IndicatorName'].iloc[0])\n",
    "plt.ylabel('# of Countries')\n",
    "plt.title('Histogram of CO2 Emissions Per Capita')\n",
    "\n",
    "#plt.axis([10, 22, 0, 14])\n",
    "plt.grid(True)\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "So the USA, at ~18 CO2 emissions (metric tons per capital) is quite high among all countries.\n",
    "\n",
    "An interesting next step, which we'll save for you, would be to explore how this relates to other industrialized nations and to look at the outliers with those values in the 40s!"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<p style=\"font-family: Arial; font-size:2.0em;color:blue; font-style:bold\">\n",
    "Matplotlib: Basic Plotting, Part 2</p>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Relationship between GPD and CO2 Emissions in USA"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# select GDP Per capita emissions for the United States\n",
    "hist_indicator = 'GDP per capita \\(constant 2005'\n",
    "hist_country = 'USA'\n",
    "\n",
    "mask1 = data['IndicatorName'].str.contains(hist_indicator) \n",
    "mask2 = data['CountryCode'].str.contains(hist_country)\n",
    "\n",
    "# stage is just those indicators matching the USA for country code and CO2 emissions over time.\n",
    "gdp_stage = data[mask1 & mask2]\n",
    "\n",
    "#plot gdp_stage vs stage"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>CountryName</th>\n",
       "      <th>CountryCode</th>\n",
       "      <th>IndicatorName</th>\n",
       "      <th>IndicatorCode</th>\n",
       "      <th>Year</th>\n",
       "      <th>Value</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>22282</th>\n",
       "      <td>United States</td>\n",
       "      <td>USA</td>\n",
       "      <td>GDP per capita (constant 2005 US$)</td>\n",
       "      <td>NY.GDP.PCAP.KD</td>\n",
       "      <td>1960</td>\n",
       "      <td>15482.707760</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>48759</th>\n",
       "      <td>United States</td>\n",
       "      <td>USA</td>\n",
       "      <td>GDP per capita (constant 2005 US$)</td>\n",
       "      <td>NY.GDP.PCAP.KD</td>\n",
       "      <td>1961</td>\n",
       "      <td>15578.409657</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "         CountryName CountryCode                       IndicatorName  \\\n",
       "22282  United States         USA  GDP per capita (constant 2005 US$)   \n",
       "48759  United States         USA  GDP per capita (constant 2005 US$)   \n",
       "\n",
       "        IndicatorCode  Year         Value  \n",
       "22282  NY.GDP.PCAP.KD  1960  15482.707760  \n",
       "48759  NY.GDP.PCAP.KD  1961  15578.409657  "
      ]
     },
     "execution_count": 30,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "gdp_stage.head(2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>CountryName</th>\n",
       "      <th>CountryCode</th>\n",
       "      <th>IndicatorName</th>\n",
       "      <th>IndicatorCode</th>\n",
       "      <th>Year</th>\n",
       "      <th>Value</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>22232</th>\n",
       "      <td>United States</td>\n",
       "      <td>USA</td>\n",
       "      <td>CO2 emissions (metric tons per capita)</td>\n",
       "      <td>EN.ATM.CO2E.PC</td>\n",
       "      <td>1960</td>\n",
       "      <td>15.999779</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>48708</th>\n",
       "      <td>United States</td>\n",
       "      <td>USA</td>\n",
       "      <td>CO2 emissions (metric tons per capita)</td>\n",
       "      <td>EN.ATM.CO2E.PC</td>\n",
       "      <td>1961</td>\n",
       "      <td>15.681256</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "         CountryName CountryCode                           IndicatorName  \\\n",
       "22232  United States         USA  CO2 emissions (metric tons per capita)   \n",
       "48708  United States         USA  CO2 emissions (metric tons per capita)   \n",
       "\n",
       "        IndicatorCode  Year      Value  \n",
       "22232  EN.ATM.CO2E.PC  1960  15.999779  \n",
       "48708  EN.ATM.CO2E.PC  1961  15.681256  "
      ]
     },
     "execution_count": 31,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "stage.head(2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAAEWCAYAAACufwpNAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd8VFX6x/HPl957ld4FFRFCsSsqdsEOFlBRV8Wy6u7a\ndnVtq+uuddeuKCiK2NG1gQUbLfQqVSG0hN4JSZ7fH/fm54jJZBIymZTn/XrNKzNnzr3zHBLy5Nxz\n7jkyM5xzzrnCUC7RATjnnCs9PKk455wrNJ5UnHPOFRpPKs455wqNJxXnnHOFxpOKc865QuNJxblS\nRNKdkl5KdByu7PKk4oo9SQMlTZa0Q1Jq+Pw6SQrff1VSuqRt4WOupIck1Y44x2WSMiVtl7RV0kxJ\nZ+TyeTHXLWB7JOnGMM4dklIkvS3pkP09t5n9w8yuDD+ntSSTVKGAceZ4fPjv/UD4vJKkR8M2bJe0\nXNLjOZzrG0mbJFUuSCyu5PCk4oo1SbcCTwL/ApoAjYFrgCOBShFVHzGzmkBD4HKgD/CDpOoRdSaa\nWQ2gDvAyMEZSvVw+Oj91c4s9t1/mTwI3ATcC9YCOwAfA6fk5fzFxB5AE9AJqAscDMyIrSGoNHA0Y\ncFbRhueKmicVV2yFPY37gOvM7B0z22aBGWZ2sZnt2fcYM9ttZlMJfnnVJ0gw+9bJAoYDVYG20WLI\nqa6kM8Ley2ZJP0rqGhHzz5JukzQb2JHDX/kdgGHAIDP7ysz2mNlOMxtlZg+HdU6XNCPsJa2U9PeI\n47N7D1dLWi1pTZh4s9//u6TXw5ffhl83h72IwyW1k/SVpA2S1ksaJalOtH+DPPQE3jez1eH35mcz\nG7lPncHAJOBVYMh+fJYrATypuOLscKAy8GF+DzSzbcA4gr+QfyP8RX8lsB1YHO08+9aV1J0gyfyB\nIGk9D4zd57LOIIJeRx0zy9jnlCcAKWY2JcrH7iD4RVwnPM+1kgbsU+d4oAPQD7hd0ok5nOeY8Gsd\nM6thZhMBAQ8BBwCdgRbA36PEkpdJwC3h5chDsi9J7mMwMCp8nCyp8X58nivmPKm44qwBsD7yF3PY\nM9gsaZekY6IcC7Ca4PJStj6SNgNrCX7xn21mW3I5Nre6VwHPm9lkM8s0sxHAHoLLbdmeMrOVZrYr\nh/PWB9ZEC9rMvjGzOWaWZWazgTeBY/epdq+Z7TCzOcArYYx5MrMlZjYu7CGlAY/lcO78eAj4J3Ax\nkAyskvT/vRFJRwGtgDFmNg1YCly0H5/nirkCDeA5V0Q2AA0kVchOLGZ2BICkFPL+o6gZsDHi9SQz\nOyrGz86tbitgiKQbIsoqEfzln21llPNuAJpG+2BJvYGHgYPDc1cG3t6nWuRn/ALENMgvqRHwFEEP\nribBv+GmXKpnJ/OKEc+zX+8FMLNM4GngaUlVgSuA4ZKmmNkCgstdX5jZ+vDYN8Ky3w3mu9LBeyqu\nOJtI0Avon98DJdUATgS+K+SYVgIPmlmdiEc1M3szok60pb+/BJpLSopS5w1gLNDCzGoDzxFctorU\nIuJ5S4Je2b5yiuOhsLyrmdUCLsnh3NnWECSP1vuUtyFIZL/9MLNdZvY0QZLqEiaZC4BjJa2VtBa4\nGThU0qG5fKYr4TypuGLLzDYD9wLPSDpPUg1J5SR1A6rndIykypJ6EMym2kRwaagwvQhcI6l3ODW4\nejiwXjOWg81sMfAM8Kak48IpuVUUTJu+PaxWE9hoZrsl9SLny0V/k1RN0kEEkxHeyqFOGpDFbycj\n1CQYH9osqRnw5yixZgLvAg9Kqi+poqRBQBfgUwBJfwzbUVVShfDSV02CGWADgMywfrfw0Zkg0Q/O\n8x/LlUieVFyxZmaPALcAfwFSgXUEg+O3AT9GVP2LpG0El7tGAtOAI8xsRyHHk0wwrvJfgqS1BLgs\nn6e5MTz+aWAzwTjD2cBH4fvXAfeF7bkbGJPDOSaEn/0l8G8z+yKHWHcCDxJMrd4sqQ9Bku4ObAH+\nB7yXR6zXEfybzib4978eON3M1oXv7wIeJRh7Wk8ws+1cM1tGcJnrFTNbYWZrsx9h2y8u6P0zrniT\nb9LlXMkR3vOxHKiYw8wy5xLOeyrOOecKjScV55xzhcYvfznnnCs03lNxzjlXaMrc7IsGDRpY69at\nEx2Gc86VKNOmTVtvZg3zqlfmkkrr1q1JTk5OdBjOOVeiSPrdDa858ctfzjnnCo0nFeecc4XGk4pz\nzrlC40nFOedcofGk4pxzrtB4UnHOOVdoPKk455wrNJ5UnHOulFuWtp1HPltIUSzL5UnFOedKsR+X\nrGfA0z/w1tSVrN6yO+6f50nFOedKqdFTVjB4+BQa16rCB8OOpFmdqnH/zDK3TItzzpV2mVnGw58u\n4MXvlnNsx4b856LDqFWlYpF8ticV55wrRXbsyeCm0TMYvyCVy45ozV9P70yF8kV3UcqTinPOlRI/\nr9/BNa9PY3Hqdu7vfxCXHt66yGPwpOKcc4UgbdsealapQJWK5RPy+R/OXMVd78+lfDkx/LKeHNsx\nz1Xq4yLufSJJ5SXNkPRx+PpVScslzQwf3cJySXpK0hJJsyV1jzjHEEmLw8eQiPIekuaExzwlSfFu\nj3POZcvMMr6Yt5ZLX55MzwfHc8KjE/h6YWqRxrArPZPb3pnNTaNncmCTmnxy09EJSyhQND2Vm4AF\nQK2Isj+b2Tv71DsV6BA+egPPAr0l1QPuAZIAA6ZJGmtmm8I6VwOTgE+AU4BP49gW55xj/fY9vDV1\nJW9MXsGqzbtoUqsKw45vxxfz1nH5q1M5o2tT7jnzIBrWrBzXOBat28awUdNZkradYce34+YTOxbp\n+ElO4ppUJDUHTgceBG7Jo3p/YKQFd+dMklRHUlPgOGCcmW0MzzkOOEXSN0AtM5sYlo8EBuBJxTkX\nJ1lZxpNfLubZb5aSnpnFke3r87czOnNi58ZUKF+OG0/owPMTlvHfr5bw7aI07jq9MxcktSAeF1HG\nTF3J3WPnUqNyBUZe0YujOySudxIp3j2VJ4C/ADX3KX9Q0t3Al8DtZrYHaAasjKiTEpZFK0/Jofx3\nJF1N0KOhZcuWBW2Lc64M25WeyZ/ensX/5qyhf7cDuKFvB9o3qvGbOpUrlOfGEzpwetem3PHeHG57\ndw7vz1jFi4OTqFlIU3p3783k7g/nMiY5hSPb1+fxC7vRqGaVQjl3YYi5nySpuqSYR6AknQGkmtm0\nfd66AzgQ6AnUA27LPiSH01gByn9faPaCmSWZWVLDhsUjmzvnSo7UrbsZ+MJEPpm7hrtO68wTF3b7\nXUKJ1K5hDUZf1YeHzjmEScs28sw3SwsljhUbdnLusz8yJjmFG/q2Z+QVvYtVQoEoSUVSOUkXSfqf\npFRgIbBG0jxJ/5LUIY9zHwmcJelnYDTQV9LrZrbGAnuAV4BeYf0UoEXE8c2B1XmUN8+h3DnnCs38\n1VsZ8PQPLFq3necv6cFVx7SN6XJWuXJiUK+WDOh2AMO/X87a/Vwi5auF6zjjP9+xcuNOhl+WxK39\nOlG+XPGbmxStp/I10I6gZ9HEzFqYWSPgaIKB8YclXZLbwWZ2h5k1N7PWwEDgKzO7JBwnIZypNQCY\nGx4yFhgczgLrA2wxszXA50A/SXUl1QX6AZ+H722T1Cc812Dgw4L+Qzjn3L7Gz1/Hec/9SJbB29cc\nTr+DmuT7HLf264QZPD5uUYFi2JuZxaNf/MQVrybTvG41Pr7haPoe2LhA5yoK0cZUTjSzvfsWhgPm\n7wLvSirIRcJRkhoSXL6aCVwTln8CnAYsAXYCl2d/nqT7galhvfuyB+2Ba4FXgaoEA/Q+SO+cKxSf\nz1vLNa9P45BmtXlxcBKNaxXsMlOLetW49PBWvPLDcq48ug0dGu87xJyztVt28+aUFbw5ZQWp2/Zw\nQVJz7ut/cMLug4mV8rMUsqRKQEUz2xG/kOIrKSnJkpOTEx2Gc64Ym5OyhQuen0inJjV546reVKu0\nf3OaNu1I55h/fU3vNvV4aUjPXOuZGT8u3cBrE39h3IJ1ZJlxbMeGDDm8Nccf2Gi/YthfkqaZWVJe\n9aL+S0m6CfjMzH6SdCwwBsiU9ICZPVNIsTrnXLGxZssuho6YSr3qlXhxcNJ+JxSAutUrce1x7Xjk\ns5+YsnwjvdrU+12dpWnbGTZqOgvXbqNutYpceVQbLurdklb1q+/35xelvGZ/DSG4HAXwV+AsgpsT\nr41nUM45lwg79mQw9NVkdqZnMvyynoV68+LlR7ShSa0qPPTpgt9tljVhURoDnv6BtG17ePT8Q5l4\nxwnccVrnEpdQIPrsr3uAA4C7JD0AdANOBv4E1JR0t6RjiiZM55yLr8ws46bRM1i4div/vegwOjWJ\nbewjVlUrlefmkzowY8VmPp+3Fggud7303TIuf2UKzetW48Prj+TcHs2L/bhJNLn268zsXkm9wzr1\ngHfN7D5J5YBTzOy+ogrSOefi7aFPFjB+QSr39z+I4zrFZ/zi3O7Neem75Tzy2U8c3aEh94ydxzvT\nUjjloCY8esGhVK9c8tf4zevy11CgMrCB4M54gI7AS/EMyjnnitKoyb/w0vfLueyI1nFdLr5C+XLc\ndsqBLFu/g76PfsM701K46YQOPHNx91KRUCCPgfrwXpDb9ilbSHAjpHPOlXizUzbz97HzOK5TQ/52\nRpe4f94JnRvRu009ZqVs5umLunN616Zx/8yilGtSkXQV8I2ZLQ5vLhwOnAv8DAwxsxlFE6JzzsXH\n1t17uf6NGTSqWYUnLuxWJHeoS+Lly3qyc08GjQp470txFu3y100ECQRgENAVaEOw2vBT8Q3LOefi\ny8y44705rNq8i6cGdaNOtUpF9tk1KlcolQkFoieVjIg76s8gWJZ+g5mNB0rePDfnnIvw5pSV/G/2\nGv7UrxM9Wv3+vhFXMNGSSpakppKqACcA4yPeqxrfsJxzLn4WrNnKvR/N45iODfnDMW0THU6pEm2g\n/m4gGSgPjDWzeQDhnfXLiiA255wrdDvTM7j+jenUqlqRxy44lHLFcKXfkizafSofS2oF1Ay37s2W\nDFwY98iccy4O7v5wHsvW72DU0N40qBHf7X7Lomizv86JeA7BBljrgZlmti3+oTnnXOEa/v1y3pmW\nwo0ndOCI9g0SHU6pFO3y15k5lNUDukoaamZfxSkm55wrVGbGo18s4r9fL+GkLo25sW/7RIdUakW7\n/HV5TuXhJbExQO94BeWcc4UlIzOLu96fy1vJKxnYswUPDDiYCuVj3knd5VO+1wUws18KuDmXc84V\nqV3pmdzw5gzGL1jHDX3bc8tJHWPaCtgVXL6TiqROwJ44xOKcc4Vm8850ho5IZvqKTdzf/6C4runl\nfhVtoP4jgsH5SPWApkCue9M751wi7UzP4OPZa3jum6WkbNrFMxd159RDStf6WsVZtJ7Kv/d5bQSr\nFS82s/T4heScc/ljZsxZtYXRU1cyduZqtu/JoF3D6owc2os+besnOrwyJdpA/YSiDMQ556JJ2bST\nl75bzp6MLMoJykmUU3DLw5TlG5m/ZitVKpbj9EMOYFCvFvRoVdfHTxKgdCzg75wr1XbvzeTqkdNY\nkrqd2tUqYmZkGWSZkZVltKpfnfsHHMxZhx5A7ao+jyiRPKk454q9+z6ez/w1Wxl+WRJ9D2yc6HBc\nFD5Z2zlXrH04cxVvTF7BH45t6wmlBMg1qUg6JeJ5bUkvS5ot6Q1JMX9nJZWXNEPSx+HrNpImS1os\n6S1JlcLyyuHrJeH7rSPOcUdY/pOkkyNjDMuWSLo9f013zhV3S9O2c+d7c0hqVZc/9euU6HBcDKL1\nVP4R8fxRYA3B0i1Tgefz8Rk3AQsiXv8TeNzMOgCbgKFh+VBgk5m1Bx4P6yGpCzAQOAg4BXgmTFTl\ngaeBU4EuwKCwrnOuFNiVnsmwUdOpXLE8/7noMCr6XfAlQqzfpSQz+6uZ/WJmjwOtYzlIUnPgdOCl\n8LWAvsA7YZURwIDwef/wNeH7J4T1+wOjzWyPmS0HlgC9wscSM1sWTnEeHdZ1zpUCfx87j4Vrt/HY\nBYfStLZv4VRSRBuobyTpFkBALUkys+ybIWNNRk8AfwFqhq/rA5vNLCN8nQI0C583A1YCmFmGpC1h\n/WbApIhzRh6zcp/yHNcjk3Q1cDVAy5YtYwzdOZco701P4a3klQw7vh3HdWqU6HBcPkRLDi8SJIMa\nBD2IBgCSmgAz8zqxpDOAVDObFlmcQ1XL4738lv++0OwFM0sys6SGDRtGido5l2grNuzkrx/MpXeb\netx8YsdEh+PyKdrNj/fmUr4WGBzDuY8EzpJ0GlAFqEXQc6kjqULYW2kOrA7rpwAtgBRJFYDawMaI\n8myRx+RW7pwrgcyM296dTTmJxy/s5qsJl0BRv2OSTpb0rKSxkj4Mn58S7ZhsZnaHmTU3s9YEA+1f\nmdnFwNfAeWG1IcCH4fOx4WvC978KL7eNBQaGs8PaAB2AKQQTBjqEs8kqhZ8xNsZ2O+eKoTemrGDi\nsg3ceVpnDqjj4yglUbQFJZ8AOgIjCXoLEPQGbpR0qpndVMDPvA0YLekBYAbwclj+MvCapCUEPZSB\nAGY2T9IYYD6QAQwzs8wwxuuBz4HywHAzm1fAmJxzCbZq8y4e+mQhR7avz6BeLfI+wBVL+nXsfZ83\npEVm9rsLmuGMrEXhlOASJykpyZKTkxMdhnMugplx2StTmbJ8I1/cfAwt6lVLdEhuH5KmmVlSXvWi\nXf7aLalXDuU9gd0Fjsw55/bx7vRVTFiUxm2ndPKEUsJFm1J8GfCspJr8evmrBbA1fM855/Zb6tbd\n3PfRPHq2rstg30irxIs2+2s60DucQtyMYApvSjj7yznn9puZcdcHc9mTkcUj5x1KuXK+VH1JF3WV\nYkm1gWMJkooBqyV9bmabiyI451zpNnbWasbNX8ddp3WmTYPqiQ7HFYJoC0oOBqYDxwHVgOrA8cC0\n8D3nnCuwlE07+dsHczmsZR2uOKpNosNxhSRaT+UuoMe+vRJJdYHJBFONnXMu3/ZmZnHjmzMwgycv\nPIzyftmr1IiWVETOy55kkfMSKc45F5Mnxi9i+orNPDXoMFrW99lepUm0pPIgMF3SF/y6cGNL4CTg\n/ngH5pwrnX5csp5nvlnKhUktOOvQAxIdjitkuY6pmNkIIAmYAOwB0oFvCJbBf7UognPOlS4btu/h\nj2/NpG2D6txzlm9/VBpFnf1lZpsI9ilxzrn9kpVl/OntWWzetZcRV/SiWqWov35cCRVt9lcLSaMl\nfSfpTkkVI977oGjCc86VFsN/WM7XP6Xx19M707lprUSH4+Ik2jItwwkud90ANAUmSKofvtcqznE5\n50qROSlb+OdnC+nXpTGX9vFfH6VZtP5nQzN7Lnx+g6RLgG8lnUUum2E559y+tuzay3VvTKNhjcr8\n89yuBGvSutIqWlKpKKmKme0GMLPXJa0lWGreb311zuXJzPjLO7NYs3k3Y645nLrVKyU6JBdn0S5/\nvcQ+e76b2XjgfGBuPINyzpUOw3/4mc/nreP2Uw+ke8u6iQ7HFYFoC0o+nkv5DIJ7VZxzLlfTV2zi\noU8W0K9LY4b6Mixlhm8A7ZwrdJt3pnPDGzNoWqcK/zrvUB9HKUN8orhzrlBlZRm3jplF2rY9vHPt\n4dSuVjHvg1ypkWdPRdLv+q05lTnnHMAL3y3jy4Wp/PWMznRtXifR4bgiFsvlr3dzKHunsANxzpV8\nH8xYxT8/W8jpXZv6/ShlVK6XvyQdCBwE1JZ0TsRbtYAq8Q7MOVey/G/2Gm4ZM5PD29bn0fN9HKWs\nijam0gk4A6gDnBlRvg24Kp5BOedKlnHz13HT6Bl0b1mXl4YkUaVi+USH5BIk2pTiD4EPJR1uZhOL\nMCbnXAkyYVEaw0ZN56BmtXnl8p6+UGQZF8uYypJwQckXJA3PfuR1kKQqkqZImiVpnqR7w/JXJS2X\nNDN8dAvLJekpSUskzZbUPeJcQyQtDh9DIsp7SJoTHvOUvL/tXJH6cel6rh6ZTPtGNRh5eS9qVvGZ\nXmVdLH9SfAh8B4wHMvNx7j1AXzPbHq5w/L2kT8P3/mxm+w72nwp0CB+9gWeB3pLqAfcQ7O1iwDRJ\nY8Nl+Z8FrgYmAZ8ApwCf4pyLu68XpjLsjem0ql+N14b28qnDDogtqVQzs9vye2IzM2B7+LJi+Ii2\nEGV/YGR43CRJdSQ1BY4DxpnZRgBJ44BTJH0D1Mq+NCdpJDAATyrOxdXKjTt54H/z+XzeOjo0qsHr\nV/amfo3KiQ7LFROxXP76WNJpBTm5pPKSZgKpBIlhcvjWg+ElrsclZf80NuPXbYsBUsKyaOUpOZTn\nFMfVkpIlJaelpRWkKc6Vebv3ZvL4uEWc+NgEvl20nj+f3ImPbjiKRjV9Mqj7VSw9lZuAOyXtAfYC\nIuiI5LnLjpllAt0k1QHel3QwcAewFqgEvADcBtwXnvd3pyhAeU5xvBB+FklJSb5sv3P59Nnctdz/\n8XxWbd7FGV2bcudpnTmgTtVEh+WKoTx7KmZW08zKmVlVM6sVvs7Xtm1mtplgw69TzGyNBfYArwC9\nwmopQIuIw5oDq/Mob55DuXOuED03YSnXvD6NGpUr8OZVffjvRd09obhcxbSgpKS6knpJOib7EcMx\nDcMeCpKqAicCC8NxEsKZWgP4dRn9scDgcBZYH2CLma0h2L+lXxhDXaAf8Hn43jZJfcJzDSaYVOCc\nKySvT/qFhz9dyBldm/K/G4/i8Hb18z7IlWl5Xv6SdCXBJbDmwEygDzAR6JvHoU2BEZLKEySvMWb2\nsaSvJDUkuHw1E7gmrP8JcBqwBNgJXA5gZhsl3Q9MDevdlz1oD1wLvApUJRig90F65wrJ+zNS+NuH\ncznhwEY8fmE3KpT3Rc1d3hRMtopSQZoD9AQmmVm3cPmWe83swqIIsLAlJSVZcnJyosNwrlj7fN5a\nrhs1nV6t6/HK5T39DnmHpGlmlpRXvVj+9NidvaWwpMpmtpBgCRfnXCn03eI0bnhjBoc0q82LvuSK\ny6dYZn+lhGMjHwDjJG3CB8SdK5Wm/ryRq0dOo23D6rx6eU9qVPYlV1z+5PkTY2Znh0//LulroDY+\nduFcqbJq8y6eGLeId6en0Kp+dV4b2ps61SolOixXAsUyUP+amV0KYGYTssuAS+Mcm3MuzjZs38PT\nXy/l9Um/AHDZEW24vm976lX3hOIKJpa+7UGRL8LZXD3iE45zrijsTM/g+QnLeOm7Zezam8n5PVpw\n44kdaOb3n7j9FG2TrjuAO4GqkrZmFwPphHenO+dKnqVp27nmtWksTt3O6Yc05ZZ+HWnXsEaiw3Kl\nRLT9VB4CHpL0kJndUYQxOefyYeXGnXyzKI0JP6VRu2pFru/bnjYNqudY9/N5a7l1zCwqVSjHqCt7\nc2T7BkUcrSvtYrn89bGk6ma2Q9IlQHfgSTP7Jc6xOedyYGZMXLaBrxem8vVPaSxJDRYDb1GvKuu3\npfPBzFVckNScG/p2+P/lVDKzjEe/+IlnvlnKoc1r8+wlPXypFRcXsSSVZ4FDJR0K/AV4GRgJHBvP\nwJxzv5eZZfz57Vm8N2MVlcqXo3fbegzq1ZLjOjWkbYPqrN+eztNfL+GNySt4d/oqLu3TikG9WnLv\nR/P4bvF6BvVqwT1nHuT3nri4ieWO+ulm1l3S3cAqM3s5u6xoQixcfke9K6kyMrO4ecwsPpq1mhtP\n6MA1x7bNdevelE07eerLxbwzLYUsg0oVynF//4O4sGfLIo7alRax3lEfS09lWzhofwlwTDj7y7d4\nc64I7c3M4sY3Z/Dp3LXcfuqBXHNsu6j1m9etxiPnHcrVx7TjjckrGHDYAXRtXqeIonVlWSxJ5ULg\nImComa2V1BL4V3zDcs5l25ORyfVvzGDc/HX89fTOXHl025iPbd+oBnef2SWO0Tn3W7HcUb8WeCzi\n9QqCMRXnXJzt3pvJdaOm89XCVO7rfxCDD2+d6JCciyrPBSUlnSNpsaQtkrZK2hZx34pzLk7SM7L4\nw2vT+GphKv84+xBPKK5EiOXy1yPAmWa2IN7BOOcCZsbt785mwqI0Hj7nEAb28gF2VzLEsvT9Ok8o\nzhWtx8Yt4r0Zq7j1pI6eUFyJEktPJVnSWwRL3+/JLjSz9+IWlXNl2OgpK/jPV0sY2LMF1/dtn+hw\nnMuXWJJKLYLtfftFlBngScW5QvbNT6nc9cFcju3YkPsHHIykRIfkXL7EMvvr8qIIxLmybu6qLQwb\nNZ1OjWvy9MXdqeh7wrsSKJbZX80lvS8pVdI6Se9Kal4UwTlXVqzavIsrXp1K7aoVecV3XHQlWCx/\nCr0CjAUOAJoBH4VlzrlCsCR1Gxc+P5FdezN59YpeNK5VJdEhOVdgsSSVhmb2ipllhI9XgYZxjsu5\nMuHHpes555kf2b03i1FX9qZj45qJDsm5/RJLUlkv6RJJ5cPHJcCGeAfmXGn37rQUhgyfQuNaVfhg\n2BG+NpcrFWJJKlcAFwBrgTXAeWFZVJKqSJoiaZakeZLuDcvbSJoc3qX/lqRKYXnl8PWS8P3WEee6\nIyz/SdLJEeWnhGVLJN2en4Y7lyhmxpPjF3Pr27Po2boe71x7BM3rVkt0WM4Vilhmf60AzirAufcA\nfc1su6SKwPeSPgVuAR43s9GSngOGEuzZMhTYZGbtJQ0E/glcKKkLMBA4iGBcZ7ykjuFnPA2cBKQA\nUyWNNbP5BYjVuSKRnpHFHe/N4d3pKZzbvTkPnXMIlSr4LC9XesQy+2uEpDoRr+tKGp7XcRbYHr6s\nGD4M6Au8E5aPAAaEz/uHrwnfP0HBJP3+wGgz22Nmy4ElQK/wscTMlplZOjA6rOtcsbR5ZzqDh0/m\n3ekp/PHEDvz7/K6eUFypE8tPdFcz25z9wsw2AYfFcvJwDGYmkAqMA5YCm80sI6ySQjCjjPDryvAz\nMoAtQP3I8n2Oya08pziulpQsKTktLS2W0J0rVD+v38E5z/zI9F828/iFh/LHEzv6jY2uVIolqZST\nVDf7haR6xHYnPmaWaWbdgOYEPYvOOVXLPnUu7+W3PKc4XjCzJDNLatjQJ665ojVl+UYGPPMDm3am\n8/qVvTlicUJ+AAAdF0lEQVT7ML/Ny5VesSSHR4EfJb1D8Ev7AuDB/HyImW2W9A3QB6gjqULYG2kO\nrA6rpQAtgBRJFYDawMaI8myRx+RW7lyx8P6MFG57Zw7N61Zl+GU9ad2geqJDci6u8uypmNlI4Fxg\nHZAGnGNmr+V1nKSG2WMxkqoCJwILgK8JZpABDAE+DJ+PDV8Tvv+VmVlYPjCcHdYG6ABMAaYCHcLZ\nZJUIBvPH5t1k5+Ive4bXzW/NonurOrx33RGeUFyZkGtPRVKN7IH2cEbV72ZVRdbJQVNgRLinfTlg\njJl9LGk+MFrSA8AM4OWw/svAa5KWEPRQBoafPU/SmPDzM4BhZpYZfv71wOdAeWC4mc3LX/Odi49n\nvlnK4+MXcU73Zjx8jg/Iu7JDQWcghzekL4GZBD2JaWa2IyxvCxxPcBnsRTN7J8cTFFNJSUmWnJyc\n6DBcKfbapF/42wdzOfuwZjx6/qGUK+cD8q7kkzTNzJLyqpdrT8XMTpB0GvAH4MhwgH4v8BPwP2BI\nuH+9cy704cxV3P3hXE7s3IhHzuvqCcWVOVEH6s3sE+CTIorFuRLtywXruGXMLHq3qcd/L/Kl613Z\n5D/1zhWCiUs3cN2o6Rx0QC1eGtKTKhXLJzok5xLCk4pz+2nuqi1cNTKZlvWq8erlvXwvFFemeVJx\nbj+YGXe+P4calSvw2tDe1KteKdEhOZdQMf9JJakR8P+7B4ULTTpXpn0xfx2zU7bwyHldaVLbN9dy\nLpYFJc+StBhYDkwAfgY+jXNczhV7mVnGo1/8RNuG1TnnsByXnXOuzInl8tf9BMurLDKzNsAJwA9x\njcq5EuCjWatZtG47t5zUkQo+08s5ILakstfMNhAsLFnOzL4GusU5LueKtb2ZWTw2bhGdm9bitIOb\nJjoc54qNWMZUNkuqAXwLjJKUSrBcinNl1tvJKazYuJPhlyX5DY7ORYilp9If2AncDHxGsCfKGfEM\nyrnibPfeTJ76cjHdW9bh+E6NEh2Oc8VKLEnlbjPLMrMMMxthZk8Bt8U7MOeKq9cn/cLarbv508md\nfKMt5/YRS1I5KYeyUws7EOdKgh17Mnj2m6Uc1b4BR7RrkOhwnCt2oi19fy1wHdBW0uyIt2ris79c\nGfXKD8vZsCOdP53cKdGhOFcsRRuof4PgfpSHgNsjyreZ2ca4RuVcMZS6dTfPf7uMEzs3pluLOokO\nx7liKdrlLzOzn4FhwLaIR/Y+9c6VeIvWbePj2avJbV+hbFt27mXw8ClkZhl/OcV7Kc7lJq+eyhnA\nNIK96SNHJA1oG8e4nCuQtG17eHz8IjbvTOfMrgdw/IGNclwxePqKTTzz9VLGL1gHwJeHpfLwuYdQ\nucLv6+5Mz+DyV6ewLG0Hwy/rScfGNePeDudKqmibdJ0Rfm1TdOE4VzAZmVm8NukXHvtiEXsysqhV\ntSKfzFlLzSoVOP2Qpgw4rBm9Wtfj+yXreeabJUxatpHaVSty0wkdkOCJ8YtJ2bST5y9N+s2ikOkZ\nWVzz+nRmrtzMMxd356gOPjjvXDQxLSgp6RzgKIIeyndm9kFco3IuH6Ys38jdH85l4dptHNOxIX8/\nswut6lfnx6XreX/GKsbOWs3oqSupUbkC2/dk0KRWFf56emcG9WpJ9XCZ+nYNa3Dr27M4+5kfeOWy\nnrRtWIPMLOPmMTP5dlEaj5zblVP8znnn8pTrHvX/X0F6BmgPvBkWXQgsNbNhcY4tLnyP+tJjy669\n3PvRPN6bvopmdarytzO6cPJBjX9378jO9AzGzV/HhJ/S6N22HgMOa5bjZa5pv2zi6pHJZGQZz13S\ng7GzVvPmlBXcdVpnrjrGr/a6si3WPepjSSrzgIMtrCipHDDHzA4qlEiLmCeV0iEzy7jslSlMXLqB\na45tx7Dj21O10v7vtrhiw06uGDGVJanbARh2fDv+fPKB+31e50q6WJNKLJe/fgJaAr+Er1sAs3Ov\n7lz8PfL5Qr5bvJ6HzzmEgb1aFtp5W9avxrvXHsEd782mRb1q/Kmfz/RyLj9iSSr1gQWSpoSvewIT\nJY0FMLOz4hWcczn5ePZqnp+wjIt7tyzUhJKtdtWKPHNxj0I/r3NlQSxJ5e6CnFhSC2Ak0ATIAl4w\nsycl/R24CkgLq95pZp+Ex9wBDAUygRvN7POw/BTgSaA88JKZPRyWtwFGA/WA6cClZpZekHhdybBg\nzVb+/PZserSqyz1nlsgrsM6VankmFTObUMBzZwC3mtl0STWBaZLGhe89bmb/jqwsqQswEDgIOAAY\nL6lj+PbTBGuQpQBTJY01s/nAP8NzjZb0HEFCeraA8bpibvPOdK5+LZlaVSvw7MXdqVTBN8ZyrriJ\nZTvhPpKmStouKV1SpqSteR1nZmvMbHr4fBuwAIi252p/YLSZ7TGz5cASoFf4WGJmy8JeyGigv4Ip\nPn2Bd8LjRwAD8orLlUyZWcYNb85g7ZbdPHtJDxrV8v3gnSuOYvlT77/AIGAxUBW4MiyLmaTWwGHA\n5LDoekmzJQ2XVDcsawasjDgsJSzLrbw+sNnMMvYpz+nzr5aULCk5LS0tpyqumPvX5z/x3eL13Nf/\nYLq3rJv3Ac65hIjp+oGZLQHKm1mmmb0CHBfrB4S7Rr4L/NHMthJcnmpHsCXxGuDR7Ko5fXQBynOK\n/wUzSzKzpIYNG8Yauism/vvVYp6bsJSLerdkUBwG5p1zhSeWgfqdkioBMyU9QpAIqsdyckkVCRLK\nKDN7D8DM1kW8/yLwcfgyhWC6crbmwOrweU7l64E6kiqEvZXI+q4UMDMeG7eI/3y1hLMPa8Z9Z/nA\nvHPFXSw9lUvDetcDOwh+wZ+b10HhmMfLwAIzeyyiPHKti7OBueHzscBASZXDWV0dgCnAVKCDpDZh\nchsIjA1vxvwaOC88fgjwYQztcSWAmfGPTxbwn6+WMLBnC/59/qFUKO8D884Vd7H0VNYD6Wa2G7hX\nUnmgcgzHHUmQkOZImhmW3QkMktSN4FLVz8AfAMxsnqQxwHyCmWPDzCwTQNL1wOcEU4qHm9m88Hy3\nAaMlPQDMIEhiroTLyjL+/tE8Rk78hcGHt+LvZx5EuXK+ba9zJUEsy7RMAk40s+3h6xrAF2Z2RBHE\nV+h8mZbiLTPLuOv9OYyeupKrjm7Dnad19n3gnSsGCnOZlirZCQXAzLZLqrZf0TmXA7NfE8oNfdtz\ny0kdPaE4V8LEcpF6h6Tu2S8k9QB2xS8kV1Y99eUSRk9dyfXHt+fWfp08oThXAsXSU/kj8Lak7JlV\nTQmWv3fuN1K37WbBmm0c0a4+FfM5qP7OtBQeH7+Ic7o349Z+HfM+wDlXLMWyTMtUSQcCnQjuDVlo\nZnvjHpkrUbbs2sugFyaxNG0HDWpUYkC3Zpyf1IJOTfLeevf7xeu5/d3ZHNm+Pg+f09V7KM6VYDHt\n/Bgmkbl5VnRl0t7MLK5/YzorNu7k7jO6MGX5RkZM/JmXvl9O1+a1Ob9Hc07vesBvtunNtmDNVq55\nfRrtGtbg2Ut6+HpezpVwec7+Km189lfhMjP++sFcRk1ewSPndeWCpOA+1Q3b9/DBzNW8nbyShWu3\nUU6Q1Loe/bo05qQujWlVvzprtuzi7Kd/xDDev+5IDqhTNcGtcc7lplB2fgxvYGxuZitzrVTCeFIp\nXMO/X859H8/nmmPbcfupv98h0cyYt3orn89by7j561i4dhsAnRrXJD0zi7Rtexjzh8PpckCtog7d\nOZcPhTKl2MxM0geA71jkfuerhet44H/z6delMX85OecdEiVxcLPaHNysNrf268SKDTsZt2Ad4+av\nZd7qrTxzcXdPKM6VIrGMqUyS1NPMpsY9GldiLFy7lRvemEHnprV4YmC3mO94b1m/GkOPasPQo9rE\nOULnXCLEklSOB66R9DPB2l8i6MR0jWdgrvhau2U3Q19NpkaVCrw8pCfVKsU038M5VwbE8tvg1LhH\n4UqM1K27GfTiJLbs2subV/WhSW3fLMs596s852+a2S8EKxP3DZ/vjOU4V/qkbgsSSurW3Yy4oieH\nNK+d6JCcc8VMLNsJ30OwGvAdYVFF4PV4BuWKn/Xb93Dxi5NZvXk3r1zeix6t6iU6JOdcMRRLj+Ns\n4CyC8RTMbDWQ923SrtTYuCOdi1+czMpNOxl+WU96tfGE4pzLWSxJJT3cEMsAJMW066MrHTbtSOfi\nlybz84YdvDykJ4e3q5/okJxzxVgsA/VjJD1PsHXvVcAVwIvxDcslWnpGFh/MWMV/v17C2q27eWlw\nEke2b5DosJxzxVwsC0r+W9JJwFagI3C3mY2Le2QuIXbvzWT0lBW88O0yVm/ZzcHNavHIeb3o09Z7\nKM65vMV6g8EcoCrBJbA58QvHJcq23Xt5fdIKXv5+Geu3p9OzdV3+cc4hHNuxoa8a7JyLWZ5JRdKV\nwN3AVwQ3Pv5H0n1mNjzewbn427B9D6/88DMjJ/7M1t0ZHN2hAdcf357e3jNxzhVALD2VPwOHmdkG\nAEn1gR8BTyol2KrNu3jx22WMnrqCPRlZnNylCdcd346uzeskOjTnXAkWS1JJAbZFvN4GlJpVi8ua\nZWnbefabpbw/YxUA/bs149rj2tK+kc8Sd87tv1iSyipgsqQPCcZU+gNTJN0CYGaPxTE+V0gWrt3K\n018v5X+zV1OxfDku6dOKK49uQ/O61RIdmnOuFIklqSwNH9k+DL/6n7YlwOyUzfz3qyV8MX8d1SuV\n5+pj2jH0qDY0rFk50aE550qhWKYU31uQE0tqAYwEmgBZwAtm9qSkesBbQGvgZ+ACM9sUbgj2JHAa\nwfpil5nZ9PBcQ4C/hqd+wMxGhOU9gFcJZqZ9AtxkZW0ry1zsTM/gtnfn8NGs1dSqUoGbTujA5Ue2\npk6132/p65xzhSWea5ZnALea2XRJNYFpksYBlwFfmtnDkm4HbidYW+xUoEP46A08C/QOk9A9QBLB\n5bdpksaa2aawztXAJIKkcgrwaRzbVCKs3ryLK0cks3DtVm46oQNXHt2GmlUqJjos51wZELekYmZr\ngDXh822SFgDNCMZkjgurjQC+IUgq/YGRYU9jkqQ6kpqGdceZ2UaAMDGdIukboJaZTQzLRwIDKONJ\nZebKzVw1Mpld6Zm8PKQnxx/YKNEhOefKkFhWKd7vtTkktQYOAyYDjcOEk514sn/rNeO3s8pSwrJo\n5Sk5lOf0+VdLSpaUnJaWtr/NKbbGzlrNhc9PpErFcrx33RGeUJxzRS7XpCLpTElpwBxJKZKOKMgH\nSKoBvAv80cy2RquaQ5kVoPz3hWYvmFmSmSU1bNgwr5BLnIzMLB4bt4gb35xB1+a1+eC6I+nY2OdR\nOOeKXrTLXw8CR5vZQkm9gUeAY/NzckkVCRLKKDN7LyxeJ6mpma0JL2+lhuUpBJuBZWsOrA7Lj9un\n/JuwvHkO9cuE7Xsy+HZRGuPnr+Orn1LZvHMv5/VozoNnH0zlCuUTHZ5zroyKllQyzGwhgJlNDgfb\nYxbO5noZWLDPvSxjgSHAw+HXDyPKr5c0mmCgfkuYeD4H/iGpblivH3CHmW2UtE1SH4LLaoOB/+Qn\nxpLGzHh/xio+mLmaSUs3kJ6ZRZ1qFel7YCNOOagJJ3Vp7Ot0OecSKlpSaZR9g2NOr2O46fFI4FKC\ny2czw7I7CZLJGElDgRXA+eF7nxBMJ15CMKX48vBzNkq6H5ga1rsve9AeuJZfpxR/SikfpH9uwjL+\n+dlCWtevxmVHtubEzo3p3rIOFcr77s7OueJBud3WEW4jnKuC3r+SaElJSZacnJzoMPJtTPJK/vLO\nbM469ACeuLAb5cp5j8Q5V3QkTTOzpLzq5dpTKalJozQaP38dd7w3h6M7NODf5x/qCcU5V2xFvW4i\n6XhJ70qaFz7ekXRcEcXmgOSfNzLsjekcdEAtnr2kB5Uq+KUu51zxFW1K8ekEy9t/DFwEXEww7jFc\n0mlFE17Z9tPabVzx6lSa1anKK5f1pEbleC6A4Jxz+y/ab6k/AwPMbFZE2UxJyQSzrD6Ja2RlXMqm\nnQwePpmqlcoz4ope1K/hC0A654q/aNdSmuyTUAAws9lA4/iF5NZs2cXFL01mV3omI67oRYt6vjy9\nc65kiNZT2VHA99x+WLtlN4NemMSG7emMHNqLA5vUSnRIzjkXs2hJpZ2ksTmUC2gbp3jKtLVbdjPo\nxUms357OiCt60b1l3bwPcs65YiRaUukf5b1/F3YgZd26rbu56MVJpG7dzcihvejRyhOKc67kiXaf\nyoSiDKQsS90a9FDWbt3NyCt60aNVvUSH5JxzBRJtSnF/ScMiXk+WtCx8nFc04ZV+a7bsChLKlt2M\nuKIXSa09oTjnSq5os7/+QrDIY7bKQE+CFYOvjWNMZcb0FZs48z8/sG7rHl69vBc9PaE450q4aGMq\nlcwscnOs781sA7BBUvU4x1XqvTsthTvem0OT2lV446revv+Jc65UiJZUfjNSbGbXR7wsfTtdFZHM\nLOORzxby/LfLOLxtfZ65uDt1q1dKdFjOOVcool3+mizpqn0LJf0BmBK/kEqvbbv3cuWIqTz/7TIu\n7dOKkUN7eUJxzpUq0XoqNwMfSLoImB6W9SAYWxkQ78BKEzPjs7lrefizhazatIsHBhzMJX1aJTos\n55wrdNGmFKcCR0jqCxwUFv/PzL4qkshKiYlLN/DwZwuZtXIzHRrV4PUre9Onbf1Eh+Wcc3GR57K3\nYRLxRJJP81dv5Z+fLWTCojSa1q7CI+d15dzuzSnve6E450oxX0u9kO3NzOIfnyzg1R9/plaVitx5\n2oEMPrw1VSqWT3RozjkXd55UClHqtt0MGzWdqT9vYvDhrbi1XydqV62Y6LCcc67IeFIpJNN+2cR1\no6axZddenhzYjf7dmiU6JOecK3KeVPaTmTFq8gru/WgeTWtX5f3retG5qS9X75wrmzyp5FN6Rhap\n23azdstu1mzZzdcLU3lvxiqO7diQJwd2o041v+/EOVd2eVKJ0dBXpzIrZQvrt+/5TXk5wfXHt+fm\nkzr6zC7nXJkXt6QiaThwBpBqZgeHZX8HrgLSwmp3mtkn4Xt3AEOBTOBGM/s8LD8FeBIoD7xkZg+H\n5W2A0UA9gpszLzWz9Hi1p3WD6jSsWZkmtavQpFYVmtSuQtPaVWlapwq1qvhgvHPOAcjM4nNi6Rhg\nOzByn6Sy3cz+vU/dLsCbQC/gAGA80DF8exFwEpACTAUGmdl8SWOA98xstKTngFlm9mxecSUlJVly\ncnJhNNE558oMSdPMLCmvetHW/tovZvYtsDHG6v2B0Wa2x8yWA0sIEkwvYImZLQt7IaOB/pIE9AXe\nCY8fgS8d45xzCRe3pBLF9ZJmSxouKXsl5GZA5DL7KWFZbuX1gc1mlrFPeY4kXS0pWVJyWlpabtWc\nc87tp6JOKs8C7YBuwBrg0bA8pxFuK0B5jszsBTNLMrOkhg191X7nnIuXIp39ZWbrsp9LehH4OHyZ\nArSIqNocWB0+z6l8PVBHUoWwtxJZ3znnXIIUaU9FUtOIl2cDc8PnY4GBkiqHs7o6EOzZMhXoIKmN\npErAQGCsBbMLvgbOC48fAnxYFG1wzjmXu3hOKX6TYD/7BpJSgHuA4yR1I7hU9TPwBwAzmxfO5poP\nZADDzCwzPM/1wOcEU4qHm9m88CNuA0ZLegCYAbwcr7Y455yLTdymFBdXPqXYOefyL+FTip1zzpU9\nZa6nIikN+KWAhzcgmCRQWpX29kHpb6O3r+Qrrm1sZWZ5Tp8tc0llf0hKjqX7V1KV9vZB6W+jt6/k\nK+lt9MtfzjnnCo0nFeecc4XGk0r+vJDoAOKstLcPSn8bvX0lX4luo4+pOOecKzTeU3HOOVdoPKk4\n55wrNGU6qYTL76dKmhtRdqikiZLmSPpIUq2I97qG780L368SlvcIXy+R9FS430uxkJ82SrpY0syI\nR1a4rE6xbWM+21dR0oiwfEG422j2MadI+ils3+2JaEtO8tm+SpJeCctnSTou4phi+f0DkNRC0tfh\n92SepJvC8nqSxklaHH6tG5YrbMOScBuN7hHnGhLWXyxpSKLaFKkA7Tsw/P7ukfSnfc5VLH9Of8PM\nyuwDOAboDsyNKJsKHBs+vwK4P3xeAZgNHBq+rg+UD59PAQ4nWJL/U+DURLetIG3c57hDgGURr4tl\nG/P5PbyIYDM4gGoE68+1JlhXbinQFqgEzAK6JLptBWjfMOCV8HkjYBpQrjh//8LYmgLdw+c1CXZ7\n7QI8Atwelt8O/DN8flrYBgF9gMlheT1gWfi1bvi8bglsXyOgJ/Ag8KeI8xTbn9PIR5nuqVjOu1N2\nAr4Nn48Dzg2f9wNmm9ms8NgNZpapYOXlWmY20YLv/EiK0S6U+WxjpEEEWzxTnNuYz/YZUF1SBaAq\nkA5sJZcdRuMdeyzy2b4uwJfhcanAZiCpOH//AMxsjZlND59vAxYQbLrXn2BXV/jt7q79CbYpNzOb\nRLANRlPgZGCcmW00s00E/zanFGFTcpTf9plZqplNBfbuc6pi+3MaqUwnlVzMBc4Kn5/Pr/u5dARM\n0ueSpkv6S1jejGA/mGxRd6EsJnJrY6QLCZMKJa+NubXvHWAHwQZxK4B/m9lGct9htLjKrX2zCLbb\nrqBgC4ke4Xsl5vsnqTVwGDAZaGxmayD4xUzwFzzkf6fYYiPG9uWm2LcPPKnk5ApgmKRpBF3V9LC8\nAnAUcHH49WxJJ5DPXSiLidzaCICk3sBOM8u+jl/S2phb+3oBmcABQBvgVkltKT3tG07wiyYZeAL4\nkWAriRLRPkk1gHeBP5rZ1mhVcyjL946wRS0f7cv1FDmUFZv2ZSvSnR9LAjNbSHCpC0kdgdPDt1KA\nCWa2PnzvE4Jr3a8T7DyZrdjvQhmljdkG8msvBYK2l5g2RmnfRcBnZrYXSJX0A5BE8NdfbjuPFju5\ntc+CXVBvzq4n6UdgMbCJYv79k1SR4BfuKDN7LyxeJ6mpma0JL2+lhuW57RSbQrCHU2T5N/GMO1b5\nbF9uou2QW2x4T2UfkhqFX8sBfwWeC9/6HOgqqVp4Tf5YYH7Ybd0mqU84o2YwxXwXyihtzC47n+B6\nLfD/XfMS08Yo7VsB9A1nD1UnGORdSC47jBZ95LHJrX3hz2b18PlJQIaZFfuf0TCml4EFZvZYxFtj\nCXZ1hd/u7joWGBx+H/sAW8I2fg70k1Q3nEnVLyxLqAK0Lzcl4+c00TMFEvkg+Gt8DcGAWAowFLiJ\nYHbGIuBhwlUHwvqXAPMIrmk/ElGeFJYtBf4beUyiHwVo43HApBzOUyzbmJ/2ATWAt8Pv4XzgzxHn\nOS2svxS4K9HtKmD7WgM/EQwEjydYqrxYf//C2I4iuIwzG5gZPk4jmGH5JUFv60ugXlhfwNNhW+YA\nSRHnugJYEj4uT3TbCti+JuH3eivBZIsUgokWxfbnNPLhy7Q455wrNH75yznnXKHxpOKcc67QeFJx\nzjlXaDypOOecKzSeVJxzzhUaTyrOxVF4L8X3kk6NKLtA0meJjMu5ePEpxc7FmaSDCe6POYxgpdmZ\nwClmtnQ/zlnBgjvonStWPKk4VwQkPUKwmGV1YJuZ3R/u9zGMYBnzH4HrzSxL0gsESwBVBd4ys/vC\nc6QAzxOsvPuEmb2dgKY4F5Wv/eVc0bgXmE6w+GNS2Hs5GzjCzDLCRDIQeINgj42N4XJAX0t6x8zm\nh+fZYWZHJqIBzsXCk4pzRcDMdkh6C9huZnsknUiwEVNysDQUVfl1WfNBkoYS/P88gGCflOyk8lbR\nRu5c/nhSca7oZIUPCNavGm5mf4usIKkDwdpevcxss6TXgSoRVXYUSaTOFZDP/nIuMcYDF0hqACCp\nvqSWQC1gG7A1YjdD50oM76k4lwBmNkfSvcD4cAn7vcA1BBtszSdYUXgZ8EPionQu/3z2l3POuULj\nl7+cc84VGk8qzjnnCo0nFeecc4XGk4pzzrlC40nFOedcofGk4pxzrtB4UnHOOVdo/g8LRalKlt+3\nOQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f115fd61550>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# switch to a line plot\n",
    "plt.plot(gdp_stage['Year'].values, gdp_stage['Value'].values)\n",
    "\n",
    "# Label the axes\n",
    "plt.xlabel('Year')\n",
    "plt.ylabel(gdp_stage['IndicatorName'].iloc[0])\n",
    "\n",
    "#label the figure\n",
    "plt.title('GDP Per Capita USA')\n",
    "\n",
    "# to make more honest, start they y axis at 0\n",
    "#plt.axis([1959, 2011,0,25])\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "So although we've seen a decline in the CO2 emissions per capita, it does not seem to translate to a decline in GDP per capita"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### ScatterPlot for comparing GDP against CO2 emissions (per capita)\n",
    "\n",
    "First, we'll need to make sure we're looking at the same time frames"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "GDP Min Year =  1960 max:  2014\n",
      "CO2 Min Year =  1960 max:  2011\n"
     ]
    }
   ],
   "source": [
    "print(\"GDP Min Year = \", gdp_stage['Year'].min(), \"max: \", gdp_stage['Year'].max())\n",
    "print(\"CO2 Min Year = \", stage['Year'].min(), \"max: \", stage['Year'].max())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We have 3 extra years of GDP data, so let's trim those off so the scatterplot has equal length arrays to compare (this is actually required by scatterplot)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "52\n",
      "52\n"
     ]
    }
   ],
   "source": [
    "gdp_stage_trunc = gdp_stage[gdp_stage['Year'] < 2012]\n",
    "print(len(gdp_stage_trunc))\n",
    "print(len(stage))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEUCAYAAAAmxTHXAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XnYHEW99vHvHQgQ2cJmhIgGXEDZIfqqQU0AwYMcQFTQ\n1w0XOApHcCEaXEEPLxFUFDkqoMgiAiIYPCJykBBQECEhgQTZBZSAIEIgYiAh+b1/VA3pDLP0s8w8\ns9yf65rr6anu6a6enqd/3VXVVYoIzMysf40a6QyYmdnIciAwM+tzDgRmZn3OgcDMrM85EJiZ9TkH\nAjOzPudA0IMkvUjS+ZLukfQnSb+W9Mo8bxtJMyXdKekuSV+SpDzvvZJuya/rJO1QZ/33SZovaV5+\nnTzA/F03iH36qqQ9Bvq5dpP0vvz93SrpZkk/lDQ2z5sl6Y48/3ZJp1Tm5fnL8/e5QNKFkl5Qte4r\nJL2oKu2Tkj7Qnr0bOEmbSfp5nt5R0t5Nln+xpMuq0n4raYNW5rPvRYRfPfQCBPwB+FghbUfgjcAY\n4B5gz5z+AuAy4PD8/g3ABnn634A/1tnGfcDGI72vnfYC3grMAcbn96sBHwa2yu9nARPz9BrAN4Gr\nC5//Z2H6XODTVes/Cji48H514BZg9WHI+5DXUWIbBwOnlFjuRmBM4f0HgS+M9PHt5deIZ8CvYT6g\nsBtwTZ15HwHOrkp7GfDXGstuACyss56agSCf6E4CrgFuA14DXAzcBfxXYbl/5r+b5mXnAQtysFoN\nODO/nw98Ki97JvDOPL07MDfPPwNYs5CvY4Gb8rytc/qb8zbm5c+tW5XvrwOHFd4fA3ymVv6afPe/\nA6Y0mP9cIMjvVwPuBXYofi95+mPA96o+vw1wfuH9nsCZVev/NnBdzu9rc/ra+Xu6Me//fjn9YOBC\n4H+AmTXy+wFSoLkZOCen/Tvwx7ye3wLjCt/ZOcDMfLwPyekTcl7WAP4C/D1/nwcBr815nZv/VgLm\nicBbq36LC0b6f6uXXyOeAb+G+YDCEcBJdeZ9CziyRvrjwHpVaUcBP6yznvvyibZycq2crGcBX8/T\nRwIP5pPpmsADwEZ5XiUQfIZ8pZdPiusCuwBXFLY1Nv89E3gnsBbwV+CVOf1s4JOFfH0iTx9WyX8+\n0U3K0+tQdfUL7MSqV+Z/Al5SK39NvvvHgPUbzJ9FIRDktBnAQVXfy+rAJcDHa6xjHrBanj62sr+F\n9Z+ep99UOXkC/w94X+X7BO4kBYeD83HZsMZ2tgHuIAf8yjL5pKw8/VHgm3n6GFLAGANsnI/RZuRA\nkJc5mMIdAbBe5VgAewAX5endge9U5eeuyu/Hr+F/rY71EwH1+hR5Ll3SFNLdw64N1jUlIh6tkf7L\n/Hc+cGtEPJTX+Wdgc+AfhWVvBM6QNBqYERHz8nJbSvoucCnwv1Xr3wq4NyLuzO/PAg4nXQlDugOB\nVERzQJ6+FviWpHOBiyPigVV2PGKupBdK2gzYBHg8Iv4i6Xn5a/B9rELSdqQr5HWBz0fEBfUWLUyP\nkVTZxu+AH9VY/gbSlfQfSEH2tqr55+V9ukbSerkOYk9gX0lH5WXWIgU6SEH3sRrb2Q34eeUYF5Z5\nMXCBpE1JV/n3Fj5zSUQsAZZIuirns9F3tj5wlqRXkH5/owv7fkLVso+QAss/sGHnyuLecyvpqrre\nvInFBElbkq5EF+f32wM/JBUfDOaf7pn8d0VhuvJ+lQuPiLiGdOW6EDhH0gci4nFgB9LV7eE5L6tk\nueT2l1e2FxHTSVevY4DrJW1d43M/J91xHAScXy9/TbZ9K7Bz/uz8iNiRVAczptbCklYDtmPlyXxJ\nROyYX5+IiKU1PvYbUv0NwBLSSb2oOtAH6Tt7R2HdL4mIyjafqrMv9S4avku6qt8O+I+q7dfadiNf\nA66KiG1JRU5rAeT9fiD/NivWIu2vtYADQe+ZCawp6ZBKgqTXSHozqQJy10rrG0ljgJPJV1+SXkK6\non5/4Yq7ZSS9FHgkIk4nXf3uLGljYFREXAR8iXxiLbgdmCDp5fn9+4Grm2znZfnE/HVgNlArEJwP\nvJsUDCqtXJ6Xvya7dDzwDUkvLqTVCwKj8/J/jYhbmqy36LekohNIAeTlVfMPyuvfFXgiIp4ALgc+\nUWgdtlOJ7VwJHChpo/yZDXP6+qTACKkSt2g/SWvlz0wm3fEVLSbdIVUU13Vw1bK/IVW+k/P9IlLR\nn7WAA0GPiYgA3g68JTcfvZVUfvtgvm3fD/iipDtIxTc3Aqfkj38Z2Aj4Xm7GOLvBpq4qNB89e5DZ\nnQzMkzQXeAfwHWA8MCsXkZwJHF21f08DHwIulDSfdKfxgybb+WRuknkz6arysuoFIuJW0klqYaU4\nq07+yE1CJ9ZYx69JgfWy3Gz3OtKdyeWFxc6VdAupAnVt0vEoLSKeBJZK2iTvx5uqFnk8b/cHpOI9\nSFfeo4FbJC3I75tt51bgOODq/L19K886hvTd/w6oLhq8gVScdz3wtYh4sGr+VcCr82/mINIFyPGS\nriXVwRRdRg4EpDvc6yPi2Wb5tsGpVPqYWZeQNI10J3GupF8An42IuyTNAo6KiEYBvFV5OoZUxPiN\nYVznHFKT5hOAX0bElcO1bluVK4vNukyu86iYRqo0vmuEstMyEbELgKQFDgKt5TsCM7M+5zoCM7M+\n50BgZtbnuqKOYOONN44JEyaMdDbMzLrKnDlzHo2ITZot1xWBYMKECcye3faGEGZmXU3S/WWWc9GQ\nmVmfcyAwM+tzDgRmZn3OgcDMrM85EJiZ9TkHAjOzPtcVzUdt5MyYu5ATL7+DBxctYbOxY5i611bs\nv9P4kc6WmQ0jBwKra8bchRx98XyWLFsOwMJFSzj64vkADgZmPcRFQ1bXiZff8VwQqFiybDknXn7H\nCOXIzFrBgcDqenBR7ZEB66WbWXdyILC6Nhtbc5TFuulm1p1KBwJJa+fBtq1PTN1rK8aMXvWQjxm9\nGlP32mqEcmRmrVC3sljSKNJg3u8FXgM8QxoU/e/Ar4HTIqLnRkWylSoVwm41ZNbbGrUaugr4LWnw\n8AURsQJA0obAFGC6pF9ExE9an00bKfvvNN4nfrMe1ygQ7BERy6oTI+Ix4CLgIkmjW5YzMzNri7qB\noDoISHohsFZh/l9qBQozM+suTSuLJe0r6S7gXuBq4D7gshbny8zM2qRMq6GvAa8D7oyILYDdgWtb\nmiszM2ubMoFgWUT8AxglaVREXAXs2OJ8mZlZm5Tpa2iRpHWAa4BzJT0CPNvabJmZWbuUuSPYD/gX\n8CngN8A9wD7NPiRpc0lXSbpN0q2SjszpJ0q6XdItkn4haexQdsDMzIamTCD4ckSsiIhnI+KsiDgZ\n+FyJzz0LfCYiXkWqYzhc0quBK4BtI2J74E7ScwrWp2bMXcik6TPZYtqlTJo+kxlzF450lsz6TplA\n8JYaaf/W7EMR8VBE3JSnFwO3AeMj4n8jolK0dD3w4rKZtd5S6eZ64aIlBCu7uXYwMGuvRl1MfBw4\nDNhS0i2FWesywFZDkiYAOwF/rJr1YeCCOp85FDgUYNy4ccyaNWsgm7Qu8PDfFnPY1iuqUp/l4Ttu\nYtYT7r3ErF0UEbVnSOsDGwDHA9MKsxbnp4vLbSBVNF8NHBcRFxfSvwBMBA6IepnIJk6cGLNnzy67\nSesSW0y7lFoHXsC909/W7uyY9RxJcyJiYrPlGhUNRUTcBxwOLC68Kv0NlcnEaFJ3FOdWBYEPkiqc\n39ssCFjvcjfXZp2hUSD4af47B5id/84pvG9IkoAfAbdFxLcK6W8lVTbvGxH/GmS+rQe4m2uzztCo\nr6F98t8tBrnuScD7gfmS5uW0zwMnA2sCV6RYwfUR8bFBbsO6WL1urgEmTZ/prq/N2qRuHcEqC0kH\nALsCAfwuIma0OmNFriMYfjPmLuzIcQYqLYmKYyWPGb0axx+wXUfkz6ybDEcdQWVF3wM+BswHFgAf\nk/TfQ8+ijZRObrZ54uV3rBIEAJYsW86Jl98xQjky631luph4M+kBsACQdBYpKFiXanSyHemr7gcX\nLRlQupkNXZkHyu4AXlJ4vzlwS51lrQt08snWLYnM2q9MINgIuE3SLEmzgD8Bm0j6paRftjR31hKd\nfLJ1SyKz9itTNPTllufC2mrqXlvVrJDthJNtvZZEI11kZdbLmgaCiLi6HRmx9mn1yXaoLZL232m8\nT/xmbdQ0EEh6HfBd4FXAGsBqwFMRsV6L82Yt1OhkO5QTeXXzz0qLpMo2zazzlKkjOAV4D3AXMAb4\naE6zHjTUpqVu/mnWfcoEAiLibmC1iFgeET8GJrc0VzZihnoi7+QWSWZWW5nK4n9JWgOYJ+kE4CFg\n7dZmy0bKUE/km40dw8Iay3ZCiyQzq63MHcH783L/CTxFeo7gHa3MlI2coTYtdfNPs+5TJhA8CiyN\niCcj4lhgKvBga7NlI2WoJ/L9dxrP8Qdsx/ixYxAwfuwY9xNk1uHKFA1dCewB/DO/HwP8L/CGVmXK\nRs5wNC1180+z7lImEKwVEZUgQET8U9ILWpgnG2E+kZv1lzJFQ09J2rnyRtIugJuAmJn1iDJ3BJ8E\nLpRUqRfYFDiodVkyM7N2KtPFxI2Stga2Io0rfntELGt5zqwjdeqANmY2eGXuCMgn/gUtzosNQjtP\nzO4+wqw3lXqy2DpTu0cac/cRZr2pYSBQsnm7MmMD0+4Ts7uPMOtNDQNBHp5yUAPVS9pc0lWSbpN0\nq6Qjc/q78vsVkpoOqmz1tfvE3MkD2pjZ4JUpGrpe0msGse5ngc9ExKuA1wGHS3o1qa7hAOCaQazT\nCuqdgEdJbDHtUiZNnzmsxUTuPsKsN5UJBFNIweAeSbdImi+p6ZjFEfFQRNyUpxcDtwHjI+K2iHCh\n8jCodWIGWB7RkjoDdx9h1puUSn8aLCC9tFZ6RNxfeiPSBNIdwLYR8WROmwUcFRGz63zmUOBQgHHj\nxu1y/vnnl91cX1m0ZBkPP/E0S5evQIgUAla1xmqj2OpF645A7sxsJE2ZMmVORDQtgi/zHMH9knYF\nXhERP5a0CbBO2YxIWge4CPhkJQiUERGnAacBTJw4MSZPnlz2o31ri2mX1ggD6eGPe6dPbnNuzKxb\nNC0akvQV4HPA0TlpNPCTMiuXNJoUBM6NiIsHm0krx5W5ZjYYZeoI3g7sSxqLgIh4EGhaziBJwI+A\n2yLiW0PJpJXTzZW5M+YuZNL0mS2p5Dazxso8Wbw0IkJSAEgqOzrZJNKgNvMlzctpnwfWBL4LbAJc\nKmleROw1wHxbDcPRhfRI8BPLZiOrTCD4maRTgbGSDgE+DJze7EMR8XtS8XQtvyifRRuIbuxCutGD\ncd22L2bdqExl8TckvQV4Engl8OWIuKLlObO+4SeWzUZWqU7ngPmkkckiT5sNGw94bzayyrQa+ihw\nA+lp4HeSHi77cKszZv2jmyu5zXpBmTuCqcBOEfEPAEkbAdcBZ7QyY9Y/urWS26xXlAkEDwCLC+8X\nA39tTXasX3VjJbdZrygTCBYCf5R0CamOYD/gBkmfBvAzAmZm3a1MILgnvyouyX/deY2ZWQ8o03z0\n2HZkxMzMRkbZ5qPW4zwovVn/ciAwd/Fg1uc8eL15UHqzPlfmgbITJK0nabSkKyU9Kul97cictYe7\neDDrb2XuCPbMA8rsQ3qm4JWkh8ysR3gcA7P+ViYQjM5/9wbOi4jHWpgfGwHu4sGsv5WpLP4fSbcD\nS4DD8lCVT7c2W9ZO7uLBeplbxDXXdPB6AEkbAE9GxHJJLwDWi4i/tTx32cSJE2P27Jpj3JuZ1VXd\nIg7S3e7xB2zXF8FA0vAMXp+9Cpggqbj82YPKmZlZm3jQo3KaBgJJ5wAvA+YBlW80cCAwsw7nFnHl\nlLkjmAi8OsqUIZnZsHC59vDwoEfllGk1tAB4UaszYmZJpVx74aIlBCuf9J4xd+FIZ63ruEVcOWXu\nCDYG/iTpBuCZSmJE7NvoQ5I2JxUfvQhYAZwWEd+RtCFwATABuA84MCIeH1TuzbpQs6t9l2sPH7eI\nK6dMIDhmkOt+FvhMRNwkaV1gjqQrgIOBKyNiuqRpwDTgc4PchllXKdOvk8u1h5cHPWquadFQRFwN\n3E4af2Bd4Lac1uxzD0XETXl6MXAbMJ40sM1ZebGzgP0Hl3Wz7lOmXyc/6W3t1vQ5AkkHAicCswAB\nbwSmRsTPS29EmgBcA2wL/CUixhbmPR4RG9T4zKHAoQDjxo3b5fzzzy+7ObOONX/hE3XnbTd+fQAW\nLVnGwseXsKLwvzlKYvwGYxg7ZnS9j5s9z5QpU0o9R1AmENwMvCUiHsnvNwF+GxE7lMmIpHWAq4Hj\nIuJiSYvKBIIiP1C2kluTdLdJ02fWbMUyfuwYrp2223PvfZxtOAznA2WjKkEg+wclu6+WNBq4CDg3\nIi7OyQ9L2jQiHpK0KfBI/TVYkccN6H5T99qq5pOu1a1YXK5t7VTmhP4bSZdLOljSwcClwGXNPiRJ\nwI9IdQrFAe5/CXwwT3+QlWMgWxMeN6D77b/TeI4/YDvGjx2DSHcC/dLdgXWuMmMWT5V0ALArqY7g\ntIj4RYl1TwLeD8yXNC+nfR6YDvxM0keAvwDvGlTO+1A/tSbp5aIRX+0PTS//NkZKmS4mvh4RnwMu\nrpFWV0T8nhQ4atl9QLk0oH+eknQRWHdqxwnav43WKFM09JYaaf823Bmx5vrlKUkXgXWfdj0N3Ym/\njRlzFzJp+ky2mHYpk6bP7MonwOveEUj6OHAYsKWkWwqz1gWubXXG7Pn65SnJfioC6xXtehp6ML+N\n4p3K+mNGI8Gify0blv+fXrlDaVQ09FNSpfDxpKd/KxZ7lLKR0w/ly/1SBNZLBhu8B1qcNNDfRvWJ\netGSZc/NG46T9kACYCfXbdQtGoqIJyLivoh4T0TcX3g5CFhL9UsRWC8ZzNPQgylOGuhvo9aJuqhW\nsdJAinrKBsBO70iw1PMAZu3kJpbdZzDBezDl/QP9bZQpTiwuM9ATdtkA2Il1G0VlRygza6t+KALr\ndAMpyhhM/dVgi5MG8tuoV5RUvUzFQOs6yj4g2On1XmWaj24BPBQRT+f3Y4BxEXFfi/NmZiNkMJWg\nAw3e7agLqnWiLqo+aQ/0hF02AHZ6vVeZoqELSeMJVCzPaWbWo9pRlNGOuqDqoqSxY0azwQtG1y1W\nGkxdx/47jefaabtx7/S3ce203ereOXRyvVeZoqHVI2Jp5U1ELJW0RgvzZGYjrB1FGe1qDj2QO5Wy\nRT2DyQN0btPvMoHg75L2jYhfAkjaD3i0tdkys2bKluEPptliu4oyOq0uqJUn7E7b16IygeBjwLmS\nTiF1GfFX4AMtzZVZF2pnO/GyZfiDfeCpVVfG3aCTT9it0nQ8gucWTOMKKI821lYej8A6XfUJF9KJ\nc6jNXusFl7LjGtRbbuyY0ay95uoNg1YnPwBl5Qx5PAJJ74uIn0j6dFU6AFVdS5v1tVZ0sdDoar5s\nGX695RYtWfbcU7b17hL68cq4XzVqNbR2/rtunZeZZa2oXG0UXMq2bilbpt9JDzdZ0s7O7OreEUTE\nqZJWA56MiJNalgOzHtCKytVGweWkg3YsVYbfrB19me1Z+7W7M7uGzxFExHJg32HfqlmPaUU78UZX\n/WW7Wqi13AYvGD2g7Vn7tbtLijKthq7LLYYuAJ6qJEbETS3JkVkXakWzw1pX8yJdHU6aPpOpe221\nSsVwo7w1akkE/dMiqJMVK+frNeFp1V1bmUDwhvz3q4W0AJr/Am3I3HKjewx35WoxuCxctATBcyeI\noRQVdPrDTf2oVnCupVV3bWUCwUci4s/FBElbtiQ3topeGfTCyqsV+K+dtlvNZqBDaZXkFkGdpVl3\n2dDau7YyfQ39vEaa+xpqg07vutaGV6MukDu990obmkbHsR1dsTd6jmBrYBtgfUkHFGatB6zVbMWS\nzgD2AR6JiG1z2g7AD4B1gPuA90bEk4POfY/zP39/adZctJN7r7ShqXd8qx8QbJVGdwRbkU7kY4F/\nL7x2Bg4pse4zgbdWpf0QmBYR2wG/AKYOML99ZTA9IVpt3TDAeKPA3+m9V9rQjPTxbfQcwSXAJZJe\nHxF/GOiKI+IaSROqkrcCrsnTVwCXA18a6Lr7RT/39zKcuqWupdFVvyt4e9tIH9+mfQ1JeiXwfdJg\nNNtK2h7YNyL+q+nKUyD4VaFo6Drg6xFxSe664tiIqPmUsqRDgUMBxo0bt8v5559ffq96yKIly3j4\niadZunwFa6w2inHrr8XYMbXbgVttd/xtMUuXr3he+hqrjWKrF638+Y30d71oyTIWPr6EFYX/yVES\n4zcY42NugzJlypRSfQ2VCQRXk4pwTo2InXLagsrJvclnJ7BqINgaOBnYCPglcEREbNRsPe50zoZi\ni2mX1myXLeDe6W8DWtdp3EC5ubANpyF3Olfwgoi4odLZXPbsYDIVEbcDe+YMvhJ422DWYzYQZSpa\nW9Fp3GC4WaeNhDLNRx+V9DLysyyS3gk8NJiNSXph/jsK+CKpBZFZS5WpiHMLLetnZe4IDgdOA7aW\ntBC4F3hfsw9JOg+YDGws6QHgK8A6kg7Pi1wM/HgwmTYbiDIVcW6eaf1sIAPTrA2M8sA01ovK1BG4\n/N66zbDVEUgaSxqacgKwemFgmiOGmEezjtHsrqFbmqCaDUaZoqFfA9cD84Hnt8Ez6xGNKmo7pTLZ\nrBXKBIK1IuLTzRezwXBxQ3dwZbL1sjKB4BxJhwC/Ap6pJEbEYy3LVZ/o9+KGMkGwUwKlK5Otl5Vp\nProUOBH4AzAnv1xzOwz6uXfRRj1tDmSZdhnpvmDMWqlMIPg08PKImBARW+SXxyMYBv1c3FAmCHZS\noCw7NKRZNypTNHQr8K9WZ6Qf9XNxQ5kg2GmB0k/9Wit0QvFnmTuC5cA8SadKOrnyanXG+kE/FzeU\n6WLb3XBbr+uU4s8ygWAGcBxwHSvrCOa0MlP9op+LG8oEwX4OlNYfOqX4s2nRUESc1Y6M9Kt+LW4o\n0+3DSPfRbtZqnVL82Wioyv8h9TH0m4hYVjVvS+Bg4L6IOKOlObSeVSYI9mugtP7QKfWEjYqGDgHe\nCNwu6UZJv5Y0U9K9wKnAHAcBM7PB65Tiz0ZDVf4N+Czw2TzAzKbAEuDOiHArIjOzIeqU4s8yzUeJ\niPuA+1qaEzOzPtQJxZ9lWg2ZmVkPcyAwM+tzpYqGKiRtAGweEbe0KD99oROeJDQzqygzMM0sYN+8\n7Dzg75KudtfUg9PvPY6aWecpUzS0fkQ8CRwA/DgidgH2aG22elenPEloZlZRJhCsLmlT4EDSmAQ2\nBJ3yJKGZWUWZQPBV4HLg7oi4MT9VfFezD0k6Q9IjkhYU0naUdL2keZJmS3rt4LPendyRmpl1mqaB\nICIujIjtI+Kw/P7PEfGOEus+E3hrVdoJwLERsSPw5fy+r3TKk4RmZhVlKos3IXU3MaG4fER8uNHn\nIuKa/ETyKsnAenl6feDB8lntDZ3yJKGZWUWZ5qOXAL8Dfksam2AoPglcLukbpLuRN9RbUNKhwKEA\n48aNY9asWUPcdOcYCxz3ulHA2inhibuYNatpaZuZWUsoIhovIM3LRTkDX3m6I/hVRGyb358MXB0R\nF0k6EDg0Ipq2QJo4cWLMnu1hks3MBkLSnIiY2Gy5MpXFv5K09zDkCeCDwMV5+kKg7yqLzcw6TZlA\ncCQpGDwtaXF+PTnI7T0IvDlP70aJ1kdmZtZaZUYoW3cwK5Z0HjAZ2FjSA8BXSJXO35G0OvA0uQ7A\nzMxGTqm+hiTtC7wpv50VEU0fLIuI99SZtUvJvJmZWRs0LRqSNJ1UPPSn/Doyp5mZWQ8oc0ewN7Bj\nRKwAkHQWMBeY1sqMmZlZe5Qdj2BsYXr9VmTEzMxGRpk7guOBuZKuAkSqKzi6pbkyM7O2KdNq6Lw8\nJsFrSIHgc3lgezMz6wF1i4YkbZ3/7gxsCjwA/BXYLKeZmVkPaHRH8GlSO/9v1pgXpAfCrMBDUJpZ\nUbecE+oGgog4NP+d0r7sdC8PQWlmRd10TijzHMG7JK2bp78o6WJJO7U+a93FQ1CaWVE3nRPKNB/9\nUkQslrQrsBdwFvCD1mar+3gISjMr6qZzQplAUAlpbwO+HxGXAGu0LkvdyUNQmllRN50TygSChZJO\nJQ1e/2tJa5b8XF/xEJRmVtRN54QyD5QdSBp7+BsRsUjSpsDU1mar+3gISjMr6qZzQpkRyl4GPBAR\nz0iaDGwPnB0Ri9qQP8AjlJmZDcZwjlB2EbBc0suBHwFbAD8dYv7MzLrejLkLmTR9JltMu5RJ02cy\nY+7Ckc7SoJQpGloREc9KOgD4dkR8V9LcVmfMzKyTddNzAs2UuSNYJuk9wAeAyoA0o1uXJTOzztdN\nzwk0UyYQfAh4PXBcRNwraQvgJ63NlplZZ+um5wSaKdP76J+AIwrv7wU8QpmZ9bXNxo5hYY2Tfic+\nJ9BM3UAg6WcRcaCk+aRO5p6bBUREbN9oxZLOAPYBHomIbXPaBUClEe1YYFFE7DiUHTAzGwlT99pq\nlToCaPycQCd3QNfojuDI/HefQa77TOAU4OxKQkQcVJmW9E3giUGuu+06+SCaWfsN5DmBTq9YbtT7\n6EP57/0AktZrtHyNz18jaUKteZJEelCtK7qy7vSDaGYjY/+dxpc6BzSqWO6Ec0iZ3kf/Q9LDwC3A\nnPwa6tNdbwQejoi7hrietuil1gFm1n6dXrFc5gr/KGCbiHh0GLf7HuC8RgtIOpQ0MA7jxo1j1qxZ\nw7j5gXn35oth81pzFo9ovsysO0zbcQVLl694Xvoaq43qiHNImUBwD/Cv4dqgpNWBA4BdGi0XEacB\np0HqYmLy5MnDlYUB+8L0mTVbB4wfO4ZPvHdy+zNkZl1lUVXxMqSK5eMP2I7JHVA0VCYQHA1cJ+mP\nwDOVxIg4ov5HGtoDuD0iHhjk59tuoK0DzMyKOr0DujKB4FRgJjAfeP69TR2SzgMmAxtLegD4SkT8\nCHg3TYqFOk2nH0Qz63xlK5ZHQpneR6+LiDe0KT81ufdRM7OBG87eR6+SdKikTSVtWHkNQx7NzKwD\nlCka+r+PVV3bAAAQzElEQVT579GFtAC2HP7smJlZu5Xpa2iLdmTEzMxGRpkHyl4g6YuSTsvvXyFp\nsN1OmJlZhylTR/BjYClQqTB+APivluXIzMzaqkwgeFlEnAAsA4iIJaQeSM3MrAeUqSxeKmkMuSvq\nPJj9M40/0p3cw6iZ9aMygeArwG+AzSWdC0wCDm5lpkaCexg1s37VtGgoIq4g9Q10MOmJ4IkRMau1\n2Wo/9zBqZv2q1PgCEfEP4NIW52VEdXo3sWZmrVKmsrgv1BtntBvHHzUzGwgHgmzqXlsxZvRqq6S5\nh1Ez6welh57sde5h1Mz6Vd1AIGk74HRgPHAZ8LmIeDzPuyEiXtueLLZPJ3cTa2bWKo2Khr4PHANs\nB9wJ/D4/QwAwusX5MjOzNmlUNLRORPwmT39D0hzgN5LeT364zMzMul+jQCBJ60fEEwARcZWkdwAX\nAR6PwMysRzQqGvo68KpiQkTcAuwOXNzKTJmZWfvUvSOIiJ9WpiWtk5LiqYj4C3BIOzJnZmat1/A5\nAkkfl/QX4H7gr5Lul3RYe7JmZmbt0Kj56BdJYxBMjog/57Qtge9I2jAiGo5JIOkMYB/gkYjYtpD+\nCeA/gWeBSyPis0Pfjedr1pOoexo1M0saVRa/H9ghIp6uJETEnyUdCNxM88FpzgROAc6uJEiaAuwH\nbB8Rz0h64WAz3kiznkTd06iZ2UoNi4aKQaCQtgRY0WzFEXEN8FhV8seB6RHxTF7mkfJZLa9ZT6Lu\nadTMbKVGdwQPSNo9Iq4sJkraDXhokNt7JfBGSccBTwNHRcSNtRaUdChwKMC4ceOYNWtW6Y28e/PF\nsHmtOYuZNWtW0/lmZv1EEbWfDZO0DXAJ8HtgDukhsteQBqbZLyJubbpyaQLwq0odgaQFwEzgyLyu\nC4Ato14msokTJ8bs2bPL7REwafpMFtboPnr82DFcO223pvPNzHqBpDkRMbHZcnWLhvKJflvgGmAC\nsGWe3rZMEKjjAeDiSG4gFTFtPMh11dWsJ1H3NGpmtlKjVkMvB8ZFxBlV6W+U9GBE3DOI7c0AdgNm\nSXolsAbw6CDW01CznkTd06iZ2UqNioZ+BXw+P01cTJ8IfCUi/r3hiqXzgMmkK/6HSWMfnwOcAewI\nLCXVEcxslsmBFg2ZmVn5oqFGlcUTqoMAQETMzmX/DUXEe+rMel+zz5qZWfs0aj66VoN5Hr/RzKxH\nNLojuFHSIRFxejFR0kdIrYg6np8eNjNrrlEg+CTwC0nvZeWJfyKpgvftrc7YUPnpYTOzcho1H304\nIt4AHAvcl1/HRsTrI+Jv7cne4PnpYTOzcpoOXh8RVwFXtSEvw+rBGg+MNUo3M+tXDfsa6mabja1d\nn10v3cysX/VsIPDTw2Zm5TQtGupWfnrYzKycng0EkIKBT/xmZo31bNGQmZmV40BgZtbnHAjMzPqc\nA4GZWZ9zIDAz63N1xyPoJJL+DtxfYtGNacFANyPE+9KZvC+dp1f2A4Z/X14aEZs0W6grAkFZkmaX\nGYShG3hfOpP3pfP0yn7AyO2Li4bMzPqcA4GZWZ/rtUBw2khnYBh5XzqT96Xz9Mp+wAjtS0/VEZiZ\n2cD12h2BmZkNkAOBmVmf68hAIOkMSY9IWlBIO0bSQknz8mvvwryjJd0t6Q5JexXS35rT7pY0rZC+\nhaQ/SrpL0gWS1mjRfmwu6SpJt0m6VdKROX1DSVfk7V8haYOcLkkn5/zeImnnwro+mJe/S9IHC+m7\nSJqfP3OyJLV5X7rxuKwl6QZJN+d9ObbR9iWtmd/fnedPGOw+tnFfzpR0b+G47JjTO/Y3lre1mqS5\nkn6V33fdMWmwL517TCKi417Am4CdgQWFtGOAo2os+2rgZmBNYAvgHmC1/LoH2BJYIy/z6vyZnwHv\nztM/AD7eov3YFNg5T68L3JnzewIwLadPA76ep/cGLgMEvA74Y07fEPhz/rtBnt4gz7sBeH3+zGXA\nv7V5X7rxuAhYJ0+PBv6Yv++a2wcOA36Qp98NXDDYfWzjvpwJvLPG8h37G8vb+jTwU+BXjX4TnXxM\nGuxLxx6TjrwjiIhrgMdKLr4fcH5EPBMR9wJ3A6/Nr7sj4s8RsRQ4H9gvR87dgJ/nz58F7D+sO5BF\nxEMRcVOeXgzcBozPeT6rxvb3A86O5HpgrKRNgb2AKyLisYh4HLgCeGuet15E/CHSr+PsEdiXejr5\nuERE/DO/HZ1f0WD7xeP1c2D3nN8B7WOb96Wejv2NSXox8Dbgh/l9o99Exx6TWvvSxIgfk44MBA38\nZ751OkO5OIV0MvprYZkHclq99I2ARRHxbFV6S+Vb151IV2zjIuIhSCdY4IV5sYHuy/g8XZ3eUlX7\nAl14XPJt+zzgEdI/2D0Ntv9cnvP8J3J+B7qPLVG9LxFROS7H5eNykqQ1q/elZJ7b+Rv7NvBZYEV+\n3+g30dHHhOfvS0VHHpNuCgTfB14G7Ag8BHwzp9cqG4tBpLeMpHWAi4BPRsSTjRatkdbp+9KVxyUi\nlkfEjsCLSVeLr2qw/a7aF0nbAkcDWwOvIRUtfC4v3pH7Imkf4JGImFNMbrDtjtwPqLsv0MHHpGsC\nQUQ8nH/wK4DTSf+8kKLh5oVFXww82CD9UdKt1+pV6S0haTTpxHluRFyckx/Ot3fkv4/k9IHuywN5\nujq9JWrtS7cel4qIWATMIpXN1tv+c3nO89cnFV0OdB9bqrAvb81FeRERzwA/ZvDHpV2/sUnAvpLu\nIxXb7Ea6qu7GY/K8fZH0k44+JkOpYGjlC5jAqpXFmxamP0UqBwTYhlUrh/5MqhhaPU9vwcrKoW3y\nZy5k1Qqow1q0DyKV3327Kv1EVq0sPiFPv41VK41uiJWVRveSKow2yNMb5nk35mUrlUZ7t3lfuvG4\nbAKMzdNjgN8B+9TbPnA4q1ZM/myw+9jGfdm0cNy+DUzv9N9YYZ8ms7KCteuOSYN96dhj0tIvYQhf\n3nmkYoZlpOj3EeAcYD5wC/BLVj0BfYFUxnsHhdpzUm38nXneFwrpW5Jq3e/OP7Q1W7Qfu5Ju2W4B\n5uXX3qSyzCuBu/LfysEV8N85v/OBiYV1fTjn927gQ4X0icCC/JlTyE+Lt3FfuvG4bA/MzXleAHy5\n0faBtfL7u/P8LQe7j23cl5n5uCwAfsLKlkUd+xsrbG8yK0+eXXdMGuxLxx4TdzFhZtbnuqaOwMzM\nWsOBwMyszzkQmJn1OQcCM7M+50BgZtbnHAh6mKRxkn4q6c+S5kj6g6S353mTJT2Re0e8Q9I1+YnI\nymeLvYoukLTvyO3JwEj6taSx+XXYID4vSTMlrdeK/FVt6/ND/Pz+kl5dZ96nJf0pd2lwpaSXFuYN\nqFdLNehltvDZyco9bRbSzpT0zjy9T/693Zzz9R9Vyx5T9X6N/LtcHWspB4Ielf+BZwDXRMSWEbEL\n6cGb4hOJv4uInSJiK+AI4BRJuxfmnxSp64J3AWdIGrbfSz7ZtuT3FxF7R3rKdiypl8qB2hu4ORp3\nBzJchhQISJ2N1QwEpOcLJkbE9qSO2U6A1A068BXg/5Cebv1KoY+o7wOHAq/Ir7cW1ndSROyYX78e\nSCbzU+mnAf8eETuQ+qqaledtJuky4LAcZD4FEKlzuCuBgwayLRs4B4LetRuwNCJ+UEmIiPsj4ru1\nFo6IecBXgf+sMe824Flg42J6vko8J1893yXpkMK8qZJuzFejlT7yJyiNZ/A94CZWfXweSa+RdF2+\nYrxB0rr5M7+TdFN+vSEvOzlfLf4iX13+oBJYJN0naWNgOvCyfHI5UdI6+cr4pnzVW6/3yfcClxTy\n9YG8HzdLOienvTSvq3K1/ZKcfma+kr4u34lVroY3zfmt3GG9UdJ0YExOOzcvN0Pp7u1WSYcW8vBP\nScflPFyvdLf3BmBf4MS8jpdVHberIuJf+e31rLwIGImeRtclPd37j5y3ZyLijjzvk8Ac4Hukfnh+\nU/jcDNLxsBZyIOhd25BOtgNxE6lTrFVI+j+kXhT/XuMz25MekX898OV8dbcn6WrytaTO6HaR9Ka8\n/FakLnd3ioj7C9tYA7gAODJfMe4BLCH1w/SWiNiZdGV4cmHbrwU+A2xH6vjugKq8TQPuyVewU4Gn\ngbfndU0Bvlkp+qgyiXRiQtI2pCdVd8v5OjIvc0rej+2Bc6vytSnpSex9SMEI4P8Cl+c7rB2AeREx\nDViS81c52X04371NBI6QtFFOXxu4PufhGuCQiLiO9DT31LyOe2rsS8VHSF0RwOB7tazVy2wpEfFY\nzuv9ks6T9N7CHeFS0tP2oyNiWb7wqFhACg7WQg4EfULSf+eryRsbLVb1/lNK3Rt/Azgoaj+GfklE\nLImIR4GrSCfnPfNrLiuDyyvy8vdH6nO92lbAQxFxI0BEPBmpe+HRwOmS5pO6FCgWg9wQqX/55aRu\nSXZtsG+V/ft/km4Bfks6yY2rsdyGkcZcgNwfft6/ygkNUuD7aZ4+p2rbMyJiRUT8qbD+G4EP5XLw\n7Qrrr3aEpJtJV/Cbs/J7WwpUyt/nkPriKkXS+0iB5cRKUo3FmvVqWa+X2VrL1kyPiI8Cu5O6hDgK\nOCPPP5F0LvqPfHc1+bkPpmO7VNK6ddZtw8CVML3rVuAdlTcRcXguLpnd4DM7kQacqTgpIr7RZDvV\n//yVE8rxEXFqcYbSOAZP1VmPaqwLUkd2D5OuokeRruobbbuR95I6adslIpYp9Q65Vo3lnpU0KlKP\nqvXyVa24zDOFaUEabCnfFb0NOEfSiRFxdnEF+QS4B/D6iPiXpFmF/C0rBOLllPzflbQH6Y7mzZF6\nvYR0pT+5sNiLSeX1dXu1jIiHC+s8nZVBqegfpM7RijYk9SxLXs98YH4uYrsXODginiAFgYeAy4FL\nJL0kIirHek1WPe42zHxH0LtmAmtJ+ngh7QX1Fpa0PfAlUudXA7Gf0ri5G5FOLjeS/pk/rDR2AZLG\nS3phg3UA3A5sJuk1+TPramX3wg/lk/L7ST1JVrxWaUzbUaRio99XrXMxqWy6Yn1SP/HLJE0BXkpt\nd5A6O4NUWXlgpYhGqaIV4DpS5TukAFO97VUotdh5JCJOB35EGooVYJlSRWolf4/nILA1qXfJZqr3\nsbjNnYBTgX0j4pHCrMuBPSVtkIt49iQVWz0ELJb0ulxk9gFyXUmuP6h4O6nIptpdpGP4qsI+7wDM\ny/UzkwvL7gjcn5d7VaGYaD6pGHJ0nrcR8PeIWFbiu7BB8h1Bj4qIkLQ/cJKkz5LK959i5WAYAG+U\nNJcUIB4BjoiIKwe4qRuAS4GXAF+LiAeBB/PJ4A+5CP6fwPtIV7L18rtU0kHAdyWNIdUP7EGqQLxI\n0rtIRU/FO4o/kMrgtyOVm/+iap3/kHStpAWk8vGvA/8jaTap99Tb62TnUlJQuzsibpV0HHC1pOWk\n4q6DSa2szpA0lfTdfqjJ9zQZmCppWf4+PpDTTwNukXQTqafJj+WiqztIxUPNnE8qOjuCNB5usZ7g\nRGAd4MJ8HP4SEftGxGOSvkYK2gBfLRR5fZw0tu4Y0ndWqVc4QWmw9QDuA1Zp+gmpAjgXQ/1Y0lqk\n3oM/GhFP5KKdz0o6lXRsnyJ9j5DqZH5Mqlt5F3BcoehsCjCgFko2cO591AYtl3f/s0TxUSu2PRk4\nKiL2abbsINa9Kaki+C3DvW6rT9IxEXFMVdrFwNGFFkbWAi4aMquSi0hOVxseKLNVzCq+yS3JZjgI\ntJ7vCMzM+pzvCMzM+pwDgZlZn3MgMDPrcw4EZmZ9zoHAzKzP/X8LglGbXzq0OAAAAABJRU5ErkJg\ngg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f11610b7f98>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%matplotlib inline\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "fig, axis = plt.subplots()\n",
    "# Grid lines, Xticks, Xlabel, Ylabel\n",
    "\n",
    "axis.yaxis.grid(True)\n",
    "axis.set_title('CO2 Emissions vs. GDP \\(per capita\\)',fontsize=10)\n",
    "axis.set_xlabel(gdp_stage_trunc['IndicatorName'].iloc[0],fontsize=10)\n",
    "axis.set_ylabel(stage['IndicatorName'].iloc[0],fontsize=10)\n",
    "\n",
    "X = gdp_stage_trunc['Value']\n",
    "Y = stage['Value']\n",
    "\n",
    "axis.scatter(X, Y)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This doesn't look like a strong relationship.  We can test this by looking at correlation."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([[ 1.        ,  0.07676005],\n",
       "       [ 0.07676005,  1.        ]])"
      ]
     },
     "execution_count": 36,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "np.corrcoef(gdp_stage_trunc['Value'],stage['Value'])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "A correlation of 0.07 is pretty weak, but you'll learn more about correlation in the next course.\n",
    "\n",
    "You could continue to explore this to see if other countries have a closer relationship between CO2 emissions and GDP.  Perhaps it is stronger for developing countries?"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": true
   },
   "source": [
    "## Want more ? \n",
    "\n",
    "### Matplotlib Examples Library"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "http://matplotlib.org/examples/index.html"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "%%javascript\n",
    "IPython.OutputArea.auto_scroll_threshold = 9999;"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.5.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
